Informe
Informe
Informe
1.1. Generalidades
Ejemplo:
1
Realizando un análisis de la transmisión de calor, permite predecir la
velocidad de la transferencia térmica del agua a la barra y de esta
información se puede calcular la temperatura de la barra, así como la
temperatura del agua en función del tiempo.
Q X K A
T
x
BTu , w …………….…….…………..
h (1,1)
T
= gradiente de temperatura
x
3
Existen dos tipos de convección:
4
Figura N° 1.2 Distribución de la temperatura y velocidad de un fluido
sobre una placa plana en convección forzada
5
g
Gr 2
T L3 ………………………. (1,6)
V
Qr A1 T1 T2
4 4
………………………… (1, 9)
Qr = A Eb = A T4 (1-10)
7
1.5. Ecuación Fundamental de la Transmisión de Calor por Conducción
Para el flujo térmico de la dirección (x), la ley de Fourier viene dada por:
T Qx T
Qx kA q k (1.12)
x x
x
A
8
Energía que atraviesa Variación de la energía
Energía generada en
por conducción el + interna dentro del elemento
el elemento de control
elemento de control de control
(1,12)
Qx = qx ∆y ∆z
Qz Q q
Q Q
z z z z x y z (1-15)
z z z
z z
3
q 0 Energía generado por unidad de volumen (W/m ), (BTU/h.m3)
9
La variación U de la energía interna de dt, para el caso de sólidos y
t t t
x y z
0
t
Ó en notación simbólica:
T
(.k T) q0 C P ……………….……….. (1,23)
t
10
Nota 1: El operador Laplaciano en coordenadas cartesianas:
2 2
2
2 ……………………….. (1,25)
x 2 y 2 z
k
Nota 2: , difusividad térmica. (1.26)
Cp
x 2 y 2 z 2 t
11
Fig. Nº 1.5 Conducción tridimensional del calor a través de un elemento
Reemplazando valores
12
T T
E U mc c rz r (1.33)
p p .
t t
0
r.z r r.r
p
rrz z t
r
r
rz r
rz
r
r r r r
r 0 rz
r . r r .r z rz z
z
z0 z z z z
8. Reemplazando en 6, se tiene
13
1.5.3 Ecuación de conducción de calor en coordenadas esféricas
r = radial
= Polar, cenital o colatitud
V r sen
2
= azimutal o longitud
2
Ar r sen
1. Balance de energia :
Velocidad de Velocidad de Velocidad de Velocidad de
generación de
conducción de conducción de cambio de
(1.35)
calor de entrada calor de salida calor en el interior energía
2. Remplazando:
E
Qr
Q Q Q
r r
Q
Q
Q generado
t
(1.36)
14
3. El contenido de energía en dicho elemento y la velocidad de generación de
calor dentro del mismo se pueden expresar como:
E U T T
mCp VCp 1.37)
t t t t
4. Reemplazando se tiene:
¶T
Q - Q - Q - Q - Q - Q + q V = r Cp V
( r + Dr r) ( f + Df f ) (
q+ D q q ) 0
¶t
T T
q kA k r sen
r r
r
2
r
1 2 T
Se tiene : kr
r 2 r r
T T
q kA k rr
1 T
Se tiene : k
r 2 sen 2
15
1 q q 1 q
Lim
rsenr r r senr
r 0 3
T T
q kA k rsenr
T
ksen
1
Se tiene :
r 2 sen
16
a) Temperatura especificada constante (condición de Dirichlet)
T ()x, T0S x
dT
k qS " , x 0
dx
dT
h T T k , x0
dx
17
(b) Ambiente radiactivo
T 4 T 4 k
dT
dx
, x0
18
En ella se incorpora R"t ,c que es precisamente la resistencia térmica de
LA LB
R"t ,c
KA KB
La presión de contacto
19
1.7 Problema Resueltos
Problema N° 1
SOLUCIÓN:
1. Diagrama de flujo
20
Figura N° 1.13-b Cilindro con fuentes de calor
2. Balance de energía
Donde:
σ=
C.F1
CF2
Dado que,
22
14. Para una emisividad de ε = 0.85
16. Luego:
Problema N° 2
23
Figura N° 1.14 Tubo de acero sometido a fluidos interior y exterior
Fuente: Elaboraciòn propia, Ing Alberto Emilio Panana Girio
1. Circuito térmico
Obteniendo Ts = 502K
Luego:
24
Problema N° 3
b) El calor disipado por unidad de tiempo y por metro de longitud del tubo.
Solución.-
1. Diagrama de flujo
hC 20w / m 2 º C
4. Despejando (hr)
hr 13,88 w / m 2 K
5. Por lo tanto:
h 20 13,88 33,88 w / m 2 K
Q d e LhTtubo Text
Problema N° 4
Solución.-
1. Diagrama de flujo:
26
Figura N° 1.16 Tubo de acero sometido a fluidos interior (vapor) y exterior
(aire)
Fuente: Elaboración propia, Ing. Alberto Emilio Panana Girio
Q = hc A(T1-T2 ) + hr A(T1-T2 )
Q = (hc + hr )A(T1 -T2 )
h = hc + hr
27
7. Realizando el reemplazo de los valores de las variables se tiene
hr = 6.221125W/m2k
h = hc + hr
h = 25 + 6.221125 = 31.22W/m2k
qr/A = 0.87KW/m2
%Q = (0.87/4.37)x100 = 19.9%
%Q = (3.5/4.37)x100= 80.1%
Problema N° 5
Un tanque Cilíndrico de oxigeno liquido (LOX, por sus siglas en ingles), tienes
un diámetro de 4 pies, una longitud de 20 pies y extremos hemisféricos. La
temperatura de ebullición del LOX es de -29oF. Se busca un aislante que
reduzca la razón de evaporación en estado estacionario a no más de m = 25
lb/h. El calor de vaporización del LOX es de ΔHv = 92 Btu/lb. Si la temperatura
exterior del aislante es de 68oF Y el espesor de este no debe ser mayor a 3
pulg, ¿Qué valor debería tener su conductividad térmica?
28
Solución.-
1. Diagrama de flujo
2. Condiciones
Según el problema, hay transferencia de calor por conducción.
Asumiremos calor conductivo unidimensional en el tanque cilíndrico.
3. Datos:
r1= 2pies ε = 3pulg = 0.25pies; ε = r2-r1 r2 = 2.25pies
QK = -0.63889 Btu/s
1
2kL(T
)L n T(r2 /r ) Q
Q k 2 1
29
Problema N° 6
Solución.-
1. Diagrama de flujo:
Qk Qc Qr
30
2 4 4
8
8, 7266 10 (630 530 )
4 4
Qr A(T
)02 , 8 T0al,r1714 10
Qr 9, 408Btu / h
En el Sistema Internacional:
Q 41,1207 W ; qt 5158, 2290 W / m 2
Problema N° 7
1. Diagrama de flujo
Qk Qc Qr
r 2 alr 2
3. El calor transferido
522 497, 94
Qk 1,1 825
0, 8
Qk 27293, 06Btu / h
P1. Dos superficies perfectamente negras están dispuestas de tal manera que
toda la energía radiante que sale de una de ellas, que se encuentra a
800°C, es interceptada por la otra. La temperatura de esta última superficie
32
se mantiene a 250 °C. Calculese la transferencia de calor entre las
superficies por hora y por unidad de área de la superficie que se mantiene
a 800 °C.
P2. Una placa de metal está perfectamente aislada por una de sus caras y por
la otra absorbe el flujo radiante del sol de 700 W/m2. El coeficiente de
transferencia de calor por convección en la placa es 11 W/m2.°C y la
temperatura del ambiente 30 °C. Calcúlese la temperatura de la placa en
condiciones de equilibrio.
P4. Se deja una plancha de 1000 W sobre una tabla de planchar con su base
expuesta al aire a 20 °C. El coeficiente de transferencia de calor por
convección entre la superficie de la base y el aire circundante es 35
W/m2.°C. Si la base tiene una emisividad de 0,6 y un área superficial de
0,02 m2, determine la temperatura de la base de la plancha.
Para una pared plana simple cualquiera sometida a convección por una superficie
(Izquierda) y a ( convección + radiación ) por la otra (derecha), tal como se muestra
en la figura se tiene,
34
Figura N° 1.22 Circuito térmico de un Superficie sólida sometido a
convección interior y convección y radiación exterior
35
CAPITULO II TRANSFERENCIA DE CALOR POR CONDUCCION
Donde:
36
d 0 dT
x 0 ………………….. (2.4)
dx dx
Figura N° 2.1 Pared Plana, para evaluar el perfil de temperatura y el flujo de calor
Fuente: Elaboración propia Ing° Alberto Emilio Panana Gario
d 2T
1). 0 …………….. (2.5)
dx 2
d dT
0
dx dx
dT
C1 ……………………. (2.8)
dx
T( x ) C1 x C 2 ………………(2.9)
C 2 TO ……………………………….(2.10)
Para xL T TL
37
TL C1 X T0
T T0
C1 L ………………….. (2.11)
L
æTL - To ÷ö
Lo que da: Q = - KA (2.14)
ç
x ç L ÷
è ø
KA
La cantidad:, ), para una placa o
L es la conductancia térmica (Ct
pared plana.
L
Resistencia térmica = = Rt
KA
K K 0 1 T ,
T2
Si K Ko 1 T
T
K 1
; (2.15)
K ()T dT
m T2 dT
T1
T T T 2
T 2
T2
K 0 1 T dt 2 1
2 2 1
K K
T1
m T2 0
T2 T1
T1
dT
K m K o 1 T2 T1
2
K 0 A
Qx T T T 2 T 2 (2.16)
1 2 1 2
L 2
T1 T2
Q Km A
L
Donde:
KO = conductividad térmica en T = 0
T1 T2
2
39
2.2.2 Cilindros Huecos
dr dr
dT C1
2. Separando las variables e integrando se obtiene:
dr r
CF: 2 T( r 0) T0 , r r0
Ti C1 Lnri C2 T0 C1 Ln r0 C2
Ti TO
C1 C2 Ti Ci Ln ri (2.20)
r
Ln O
ri
40
Ti T0
C T Ln r
2 i Ln r i
/r
o i
Ti To T T
T Ln r T i 0 Ln r
Ln r / r i i
L n r /r
o i o i
T T
Ti T0 Ln
r
(2.21)
i
r0 ri
Ln
ri
dT C1 Ti T0
; Donde: C1 (2.22)
dr r r
0
Ln
ri
Ln r0 / ri
2KL
Qr T T
i 0
Ln r0 / ri
41
T
T
AK
Q T
i 0
L L
AK
Ti T0
Q
r
Ln 0 r
i
2 KL
ri
Como A 2rL y ro o A , A se puede expresar como:
A
ri i
Ao Ai
A (2.24)
Ao
Ln Ai
Ti T0 Ti To
Q (2.25)
ro ri ro ri
K A K
Ao Ai
AO
Ln Ai
A0
Nota 2: Para valores de Ai < 2, el área media aritmética
42
un 4%, por lo cual la primera puede usarse satisfactoriamente, para
paredes de mayor espesor, esta aproximación no es aceptable.
43
Ti T0 Ti T0
C → C ……. (2.31)
1 1
ro ri
1 1 1 1
ri ro
Ti To
1 1
C1 r ro
C2 Ti C2 i Ti …….. (2.32)
ri ri
dT C1 dT Ti To
r2 (2.34)
dr r2 dr 1 1
ri ro
Sustituyendo Qr queda:
4K
Qr T T
1 1 i 0
ri r0
T T ……………
4Kri r0
Q . (2.35)
r0 ri
r i 0
45
K = conductividad del aislante.
a. Circuito térmico
Qt R2
R1
L n r T2
r0
T∞
TC i 1
R1 R 2
2KL 2r0 Lh
Ln 0 r
r 1
i R2
R 2 r Lh
2 KL
1 0
2LT1 T
Qr
Ln r0 r
i 1
K ro ho
Figura N° 2.5 Pared plana compuesta con supeficies convectivas para evaluar el
coeficiente de convección global
1
U= (2.42)
1 D x23 D x34 1
+ + + ... +
h12 K 23 K 34 h(2n)+(3 ) + n+
b. Para un cilindro
47
Q U 4 A4 (T1 T5 ) ; A4 2r4 L
Basado en el radio r4
Figura N° 2.6 Cilindro hueco compuesto con supeficies convectivas para evaluar
el coeficiente de convección global
48
Suponga que dos sólidos A y B se colocan uno junto al otro en
paralelo y que la dirección del flujo de calor es perpendicular al
plano de la superficie expuesta de cada sólido.
Circuito térmico
Resistencia A = R A L
K A A1
Resistencia B = RB L
K B A2
O sea:
1 1 RA RB
Q T T T T
t 1 2 1 2
RA RB RA RB
T1 - T2
Qt = (2.45)
RA RB
RA + RB
Problema N° 1
Solución:
1. Diagrama de flujo
2. Datos:
ε = 0,6
50
CF:
CF:
51
Reemplazamos estos datos en la ecuación (3):
Despejando T1:
CF1:
CF2:
PROBLEMA Nº2
Una barra de oro está en contacto térmico con una barra de plata, una a
continuación de la otra, ambas de la misma longitud y área transversal. Un
extremo de la barra compuesta se mantiene a T1 = 80º C y el extremo opuesto
a T2 = 30º C. Calcular la temperatura de la unión cuando el flujo de calor
alcanza el estado estacionario.
52
Figura N° 2.9 Barras planas colocados en serie, para evaluar la temperatura de
contacto
Solución:
PROBLEMA N° 3
Diagrama de flujo.-
Datos.-
KA = 20w/mk KC = 50w/mk KB =?
LA = 0,3 m LB = 0, 15 m LC = 0, 15 m
Solución.-
Qconv = Qa = QB = Qc
t F ts ,o
qx
1 L L L
A B C
h K A K B KC
54
800º C 20º C 780
q ………………… ()
x
1 0.3 0.15 0.15 0.15
0.058
25 20 K B 50 KB
Problema N° 4
Una pared plana grande, tiene un espesor de 0.35m; una de sus superficies se
mantiene a una temperatura de 35ºC mientras que la otra superficie esta a
115ºC. Únicamente se dispone de dos valores de la conductividad térmica del
material de que esta hecha la pared: así se sabe que a T = 0ºC, K = 26w/mk y
a T = 100ºC, K = 32w/mk. Determine el flujo térmico que atraviesa la pared,
suponiendo que la conductividad térmica varía linealmente con la temperatura.
Solución.-
1. Diagrama de flujo
55
i. La temperatura media de la pared.
0.35 115 35
x m 0.175m , Tm 75º C
2 2
0 35 m
x m
A e
q x 6971,4285 w / m 2
3. Segunda forma de solución:
De (1) K 0 26
32
1 / 100 2.30769 x10 3
26
K (T ) 26 1 2.30769 x10 3 T K K 1
1
, 2
T T
m 0
2
56
Km
K (T )dT ; K
26 26 x2.30769 x10
T dT
dT
m
3
115
115 35
35
26 x2.30769x10
3 2 2
26(115 35) 115 35
Km 2 30.499 w / mº C
80
PROBLEMA N°5
Datos:
: AA = AD = 36 cm2 = 36 x10-4m2
Solución:
57
R B RC
R2
RB RC
T max imo
Q 3
1
Rn
1
Rn R1 R2 R3
XA 0.02m
R1 0.079365k / w
K A AA (70w / mk ) (36x10 4 m 2 )
R B RC
R2
R B RC
XB 0.025m
RB 0.23148 k / w
K B AB
60
w
mk
18x10 4 m 2
XC 0.025m
RC 0.34722 k / w
K C AC
40
w
mk
18x10 4 m 2
0.23148x0.34722 0.080374
R2 0.13888 k / w
0.23148 0.34722 0.3787
XD 0.04m
R3 0.55555 k / w
K D AD 20w / mkx36x10 4 m 2
58
1
3
w
Rn 0.079365 0.13888 0.55555 0.773795 k
1
5. Reemplazando:
150º C
Qx 193.849 w
ºC
0.773795
w
PROBLEMA N° 6
Se supone que la interfaz entre los dos materiales tiene una resistencia de
contacto infinita y que toda la superficie externa está expuesta al aire, para
el cual y
a) Dibuje el circuito térmico del sistema. Usando los símbolos precedentes,
marque todos los nodos y resistencias pertinentes.
b) Para las condiciones que se establecen. ¿Cuál es la pérdida total del
calor del tubo? ¿Cuáles son las temperaturas de la superficie externa
Ts,2(A) y Ts,2(B)?
Solución.-
Problema N°7
Vapor con calidad del 98% fluye a una presión de 1,37. 105 N/m2, a una
velocidad de 1 m/s por un tubo de acero de 2,7 cm. de diámetro exterior y 2,1
cm. de diámetro interior. El coeficiente de transferencia de calor en la superficie
interna donde ocurre condensación es de 567 w/m2k. Una película de grasa en
la superficie interna añade una resistencia térmica unitaria de 0,18 m2k/w.
Estime la razón de pérdida de calor por metro de longitud de tubo, si:
60
1. Diagrama de flujo
V g 1,2687 m 3 / kg h g 2689,3471 KJ / kg
h1 455,3593 2189,3081
h1 2644,6673 KJ / kg
V1 0,0011 1,2422
V1 1,2433m 3 / kg
(0,0105m) 2
1m 1kg
m 3
s 1,2433m
61
m 2,78581.10 4 kg / s
KJ
H 1 0,7368 736,8 watt
s
K C 43,27 w / mk
Q TV Ta
L 1 Ln(2,7 / 2,1) 1
0,18
(2,1.10 )567 2 (43,27) (2,7.10 2 )(11)
2
Q w
38,7822
L m
Si L = 3m → Q 206,3467 watt
H 2 H1 Q → H 2 736,8 206,3467
H 2 530,4533 w
H2 530,4533 w
h h
→
2
m
2
2,78581.10 4 kg / s
62
KJ
h 2 1904,1258
Kg
h2 hF 1904,1258 455,3593
x → x
h g hF 2689,3471 455,3593
x 64,8551%
Q TV Ta
1 Ln(r2 / r1 ) Ln(r3 / r2 )
L
RT
2r1 hV 2K C 2Ka
e 10cm r3 r2 e 2, 7 5 7, 7cm
Q (108, 6 21)ºC
1 Ln(2, 7 / 2,1)(7,L7n/ 2, 7)1
L 0,18
(2,1.10 2 )567 2 (43, 21)2 (0
, 07788)(7, 7.10 )(11 2
)
Q w
32,14
L m
H 2 640, 3625w
H2 640, 3625w
h2 → h2 2298, 658199Kj / kg
m 2, 78581.104 Kg / s
x 81, 76%
63
2.4 PROBLEMAS PROPUESTOS
PROBLEMA 1
k 3,35x10 4 T 0,1179
2
Caso T1 T2 dT / dC (K/m) qx (W/m )
1 400K 300K
2 100ºC -250
3 80ºC 200
4 -5ºC 4000
5 30ºC -3000
PROBLEMA 2
64
b) Evalúe cualitativamente la influencia de T, ext y hext sobre las
temperaturas.
Solución:
a) Tint 7,7º C y Text 4,9º C
PROBLEMA 3
Solución:
a) P"elec 1,27 kW / m 2
b) Text 11,1º C
c) P"elec aumenta al aumentar hext y disminuye al aumentar T, ext
PROBLEMA 4
Una casa tiene una pared compuesta de madera, aislante de fibra y tablero de
yeso, como se indica en el esquema. En un día frío de invierno los coeficientes
de transferencia de calor por convección son hext 60W / m 2 .K y
300 16 0,046
300 28 0,038
300 40 0,035
Solución: b. 4.214 W; c. 0, 45 %
PROBLEMA N° 5
Una hielera cuyas dimensiones exteriores son: 30cm x 40cm x 40cm está
hecha de espuma de estireno ( ). Inicialmente la hielera está llena
con 40 Kg de hielo a 0 ºC y la temperatura de la superficie interior se puede
tomar como 0 ºC en todo momento, el calor de fusión del hielo a 0 ºC es
y el aire ambiente circundante esta a 30°C. Descartando toda
66
transferencia de calor desde la base de 40cmx40cm de la hielera, determine
cuanto tiempo transcurrirá para que el hielo que está dentro de ella se funda
por completo, si las superficies exteriores de la misma están a 8°C
PROBLEMA N° 6
67
PROBLEMA N° 7
1. Diagrama de flujo
PROBLEMA N° 8
68
CAPITULO III CONDUCCION CON FUENTES DE CALOR (Generación
Interna de calor)
3.1 Generalidades
eléctrica o través º de la placa, provocando en esta una generación de
calor uniforme q por unidad de tiempo y volumen.
º
d 2T qº
4. Simplificando: 2
º
0 Ecuación diferencial de 2do orden (3.4)
dx K
5. Condiciones de Frontera:
dT
CF: 1 x=0 0 (3.5)
dx
xL º
dx 2
KA
dT
1
q 2LA
º
x L º
dx 2
71
71
Reemplazando esta condición en en la distribución de temperatura se
qº
obtiene: Tc Tw º L2 Tmax (3.11)
2K
qx x L K (/ qo)L K
qx x L qo L (3.13)
dT dT d dT
2. K 2 rL qºº 2 rLdr K 2 rL K 2 L r dr
dr dr dr dr
d dT q0 r
r 0 (3.16)
72
72
dr dr K
73
73
3. Las condiciones de frontera para resolver la ecuación (16), son :
dT qºº r C1
(3.19)
dr 2K r
q ºr 2
T º C1 Ln r C2 (3.20)
4K
qºº r02
Tc = T = Tw + (3.24)
max
4K
Reemplazando se tiene
74
74
q r
q0 r0
q K 0 0 q (3.26)
r ro r r ro
2K 2
75
75
3.4 ESFERA SÓLIDA
r
dr
q 4 r dr
2
incremento de energía de volumen Q
gen o
r 0 0
2
: (3.30)
dr dr K
76
76
4. Debido a que ( q0 ) es uniforme a través de la esfera y Tw es constante
sobre toda la superficie exterior de la esfera, es de esperar que la
distribución de temperatura sea simétrica con respecto al centro de la
esfera, por tanto la otra condición de frontera es:
dT
CF: 2 r=0 0 (3.32)
dr
77
77
5. Separando variable (30) e integrando, se tiene:
dT qº r 3
r2 º C1
dr 3K
dT qºr C
º 1 (3.33)
dr 3K r 2
6. Separando variables de nuevo o integrando la relación anterior:
qº r 2 C
T º 1 C2 (3.34)
6K r
q ºº r 2
O sea T C2
6K
T Tw º r r
qº 2 2
(3.37)
0
6K
10. Se puede determinar la temperatura Tc en el centro de la esfera (r=o)
2
qºº ro
Tc Tw T max (3.38)
6K
3.4.2 Flujo de calor
dT dT q 0 r0
qr K r ro (3.39)
dx dx 3K
qr q0 r0
q K 0 0 (3.40)
r
3K 3
NOTA 1.- Puede suceder que algunos problemas no se conozca Tw, pero en
cambio, qºº , h y Tf (temperatura de fluido), son conocidas. En estado
78
78
estacionario, todo el calor generado en el sólido se debe transmitir
por convección hacia fuera al fluido que lo rodea (si no fuera así, se
tendría un crecimiento de la energía en el sólido, que daría por
resultado un incremento de la energía interna del material que
posteriormente requeriría un cambio de temperatura con respecto al
79
79
tiempo), se puede determinar Tw para las tres geometrías en la
forma siguiente, si el coeficiente de transferencia de calor (h) es
uniforme:
Qgen = qo°V = h.Ssuperficie (TW - Tf ) (3.41)
Donde:
V= Volumen de todo el cuerpo
Qgen q ºº r0 L
2
q ºº r 0
T w T f (3.43)
2h
4
80
80
Energía total generado en la esfera, Qgen q ºV q º r 3 (3.44)
º º 0
3
Calor que transfiere por convección la esfera al fluido que lo rodea
81
81
Qc h.4 r02 . Tw T f
qºr
Tw º 0 T f (3.46)
3h
PROBLEMA Nº 01
Solución:
82
82
3. El calor generado por unidad de longitud y tiempo es:
2
x x 0.20m
o 2
q' E' q D 24,000 W 754.0 W
1
gen
4 m3 4 m
2 25
m 2 .Kxp x0.400m
83
83
T0 T0 T1 71.8º C 192º C
4K r 4x0.5 W
m.K
T(0) = 192 °C
84
84
PROBLEMA 2
b) Calcular los valores de los flujos de calor en las dos superficies del muro,
y .
Solución:
1. Diagrama de flujo
d 2T
k + qo = 0
dx 2
85
85
4. Al comparar de (3) con la ecuación del enunciado: T(x) = a + bx + cx2
Se tiene que:
°C/m2
5. Despejando: qo=2x105W/m3
6. Aplicando la ley de Fourier en los dos extremos de la pared:
PROBLEMA 3
1. Diagrama de flujo
86
86
Figura 3.5 Alambre cilíndrico aislado sólida con generación interna de Calor
Uniforme
2. Sea:
Tc = temperatura en el centro del alambre
Ts = Temperatura en la interfase alambre – plástico
Tp = temperatura superficial del plástico
dT qr C
o + 1
dr 2k r
q r2
T´= o C1 Lnr C2
4k
4 Cálculos en el plástico
4,1 Ecuación diferencial de conducción de calor a través del plástico
d dT
r 0
dr dr
h T´T f
dT ´
CF : 4 r r2 -k
dr
Tf Ts r
T´´ Ts Ln (ii)
Ln r2 k2 r1
r1 hr2
81
81
6. El calor generado en el alambre es igual al que se conduce a la
superficie del mismo, igual al que se conduce a través del plástico y es
igual al transferido por convección al fluido.
82
82
qo r1 dT ´´ dT ´´ T f Ts 1
k p ;
2 dr dr Ln r2 kp r1
r1 hr2
7. Reemplazando y despejando, Ts
q or12 r2 kp
Ts Ln T f
2k p r1 hr2
1, 5 106
0, 003
2
Tc 97,1
4 18
Tc 97, 28 97, 3C
25 97,1 0, 7
T 97,1 Ln
0, 3
p
0, 7 1, 8
Ln
0, 3 14 0.007
Tp 93, 92C
PROBLEMA Nº 04
Una placa plana cuyo espesor es 10 cm. genera calor a razón de 30000 w ,
m3
cuando se hace pasar una corriente eléctrica a través de ella, una de las caras
de la pared esta aislada y la otra esta expuesta al aire con temperatura de
83
83
25ºC. Si el coeficiente convectivo de transferencia de calor entre el aire y la
superficie de la placa es h = 50 w 2 , y la conductividad térmica del material
m K
K= 3 w mK . Determine:
84
84
Solución:
1. Diagrama De flujo
Figura 3.6 Pared con superficies una aislada y otra expuesto a un fluido con
generación interna de Calor Uniforme
dT qº q ºº 2
= - º x + C1 ; T= - x + C1x + C2 ;
dx K 2K
dT
CF: 2 x=L K x L h Tx T
dx
dT qo qo
L T L2 C
; ()x2
dx xL k xL
2k
85
85
Reemplazando las relaciones en la condición de frontera y
simplificando
86
86
q q
K o L h o L2 C hT
2
k 2k
qL q L2
q 2k
C T L T
o o o
1 2
2
h 2k 2k hL
qL2 x 2k
2
T T
1
2k L hL
dT
6. La Tmax estará cuando 0
dx
dT
5000(2 x)0 , estará en, x = 0
dx
7. Tmax 135 5000(0) Tmax 135
PROBLEMA 5
Dos Grandes placas de acero a 90ºC y 70ºC están separadas por una barra de
acero de 0.3 m de largo y 2.5 cm de diámetro. La barra está soldada en cada
placa. El espacio entre las placas se rellena de aislante que también aisla la
circunferencia de la varilla. Debido al diferencial de voltaje entre ambas, fluye
corriente a través de la barra, y se disipa energía eléctrica a razón de 12 W.
Calcule:
88
88
Figura 3.7 Barra de acero conectada a placas en sus extremos, con su
superficie aislada y con generación interna de Calor
Uniforme
Solución:
q 0 81487.33086w / m³
4. Condición de frontera
q 0 L²
T2 C1 L T1
2K
89
89
q0 .L²
(T - T )-
2 1
2K T2 - T1 q0 L
C1 = C1 = -
L L 2K
Reemplazando C1 y C2 en ()
q0 T2 T1 q0 .L
T .x² x T1
2K L 2K
q0 T2 T1
T x² Lx x T1
2K L
T 2829.42121x² 915.493029 x 70
x 0.1617
Tmáx 144.05º C
6. Flujo de calor
Qx dT
q x K
A dx
Para X=0.3
dT
q x | x0.3 K ……………… (β)
dx x0.3
7. Remplazando en (β)
91
91
w ºC
q x | x 0.3 14.4 x 782.1596
mº C m
q x | x 0.3 11263.0996w / m²
Qx | x 0.3 q x | x 0.3 xA
w
Qx | x 0.3 11263.0996 x x0.025² m²
m² 4
Qx | x 0.3 5.52876w
9. Para X=0
dT
q x | x0 K ……………… (γ)
dx x0
dT q0 x T2 T1 q0 .L
dx x0 K L 2K
dT ºC
915.493020
dx x0 m
q x | x 0 13183.099w / m²
Qx | x0 q x | x0 .A
w
Q x | x 0 13183.099 . .0.025²
m² 4
92
92
Qx | x 0 6.47123w
93
93
Qx | Qx | x0.3 | | Qx | x0 |
Qx | 5.52876w | | 6.47123w |
Qx 11.999w
13. Comparando el flujo de calor total con el flujo de calor por generación
Q x Qgen
11.999w 12w
PROBLEMA Nº 06
a) La temperatura máxima.
b) La generación interna de calor.
c) Los flujos de calor en x 0 y x 3m .
1. Diagrama de flujo
Figura 3.8 Cilindro sólida con generación interna de Calor Uniforme, con
superficie aislada
94
94
Fuente: Elaboración propia, Alberto Emilio Panana Girio
95
95
2. Ecuación diferencial gobernante Sistema Unidimensional con
generación de calor
d 2T qºº
0
dx 2 K
CF: 2 x1=3 T2 0º C
C2 400º C
qoº 9 qo º
0 x2 C x C 3C 400
2 100
1 2 1
2K
9 q ºo
C1 133, 333 0, 015qo 133, 3333
3 200
7. Reemplazando C1 y C2 en T Ecuación ( 3 )
x 2 0.015q o 133.333 x 400
qo
T
2(100)
Perfil de Temperatura
8. La temperatura máxima (Tmax), se encuentra cuando x 0.8m , por la
dT
condición de máximo: |x 0.8 0
dx
dT
8.1 Hallamos: de la ecuación 6. se tiene
dx
96
96
0
qº
0.8 0015qº133.333
100
91
91
qº 19047.57 w
m3
8.2 Tmax x=0.8
19047.57 2
T 0.8 [0.015(19047, 57)133, 333]0.8
400
2 100
max
Tmax 460.95º C
9. Calculo del flujo de calor en x 0 y x 3
dT
9.1 Q x 0 KA x0
dx
dT qº
x 0.015qº 133.333
dx 100
D
2
qº
Q K x 0.015qº 133.333 ()
100 4
x
100(19047.57)(0.02)
2
Q x 3 13.16
10. Comprobación
qoV q x 0 q x 3 17.947 qoV 19047.57 0.02 3 17.95
4
90
90
Problema N°7
Solucion.-
2
T q0
0
r 2 k
dT qr C
0 1
dx 2k r
q0 r 2
T C 1Lnr C 2
4k
91
91
CF :1 r r0 T TS
dT
CF : 2 r0 0
dr
q0
T T r2
C S 0
4k
q0 r0 q0 r0
h(T
) T h
) S Tf
S f
2 2(T
q0 75030012.01 W
m3
75030012.01
TC 240 (0.001) 2
4
TC 240.0467768C
92
92
75030012.010.001
h 170.4865 W
2(240 20)
m2 C
93
93
3.7 PROBLEMAS PROPUESTOS
94
94
To respectivamente. Obténgase también una expresión para la
temperatura máxima de la pared. Y el flujo de calor en la superficie.
95
95
Figura 3.11 Sistema planos en serie con generación interna de Uniforme
96
96
IV. SUPERFICIES EXTENDIDAS (ALETAS)
4.1 INTRODUCCION
4.2.1 Aletas
Las aletas pueden ser de sección transversal rectangular, como tiras que se
anexan a lo largo de un tubo, se les llama aletas longitudinales, o bien discos
anulares concéntricos alrededor de un tubo, son las aletas circunferenciales.
El espesor o el área de la sección transversal de una superficie extendida
pueden ser uniformes o variables
97
97
Figura N° 4.1 Superficie extendida de área variable
Q KA
dT
K
d
A
dT
dx (4.4)
x dx
98
98
( x)
dx dx
( x)
dx
99
99
d) La transferencia de calor por convección dQc, se expresa como
dQc hdS ( x ) (T T) (4.5)
d 2T 1 dA()x() dT 1 h dS
x
T T 0
(4.6)
2
dx A()
x() dx dx Ax K dx
10
01
00
La aleta se une a una superficie base de temperatura T To , en x=o y
10
11
01
de una aleta es uniforme, son constantes el área A( x ) A , y el perímetro
En consecuencia:
dA dS ( x )
0 P
dx dx
d dT
dx dx
Donde: m 2 hP (4.10)
KA
b. ()x3 C cosh()msex4nh(C
) mx (4.12)
L x 6 C senh m L x
c. ()x5 C cosh m (4.13)
10
21
02
Las funciones hiperbólicas se definen:
a) Distribución de Temperaturas
d 2 ()x
2
m 2 ()x 0
dx
100
100
100
1. Condiciones de Frontera
CF1: ()x To T f o en x = o
CF 2 : ()x 0 , cuando x
Ph
2. Se tiene que: m2 (4.15)
kA
3. Solución de la Ecuación Diferencial
( x ) C1e mx C 2 e mx
C2 o .
hP
Q o hPe0 dx hP e
mx mx
o , m
m
0 KA
Q hPAK 0 (4.18)
101
101
101
b.2 Cálculo de la transferencia de calor por conducción en la base de la
aleta.
102
102
102
Q KA
d
x o KA m 0 e
mx
x o
dx
Q= hPAK .Q 0
a) Distribución de Temperatura
Ecuación Diferencial
d 2 ()x
2
m 2 ()x 0 para: 0 x L
dx
Condiciones de Frontera
CF1: x To T f o x=o
dQ ()x
CF 2 : = 0 x=L
dx
m 2 hP
KA
103
103
103
xL 0 1 2
dx
o C1 C 2 C 2 o C1 (ii)
104
104
104
Resolviendo en forma simultanea (i) y (ii)
oe mL o e mL
C1 mL y C o
e e mL e mL e mL
2
e mx e mx
( x) o
1 e 1 e 2 mL
2 mL
Cosh [m()L - x ]
Q ()x = Q o (4.19)
cosh (mL)
De la ecuación de Fourier
d ( x )
Q AK x o
dx
d ( x )
Sustituyendo, de la ecuación (14) en la ecuación de flujo de calor
dx
y simplificando;
Q AK o mTanhmL
Q o PhKA Tanh mL (4.20)
Ecuación Diferencial
d ( x )
m 2 ( x ) 0
dx
106
106
106
Condiciones de Frontera
CF1: 0 x=o
CF2: Q = QL x=L
Distribución de Temperatura
b) Flujo de calor
cosh (mL)- Q L / Q 0
Qo = hPKA ´ Q 0 ´ (4.22)
senh (mL)
Ecuación Diferencial
d2 ( x )
2
m 2 ( x) o
dx
Condiciones en la Frontera
CF1: ()x0 T T f o en x=o
d ( x )
CF2: K h ( x ) O en x=L
dx
107
107
107
Las constantes de integración C1 y C2 de la solución (b) se determinan
aplicando las condiciones de frontera donde se obtienen
respectivamente:
o C1 cosh mL C 2 senhmL (i)
108
108
108
KC 2 m hC1 o (ii)
senh (mL)+ h
( mK )cosh(mL)
Qo = hPKA ´ Q 0 ´ (4.24)
cosh (mL)+ (h mK )senh (mL)
109
109
109
Figura N° 4.3 Perfil de temperaturas para aletas de área uniforme, con
diferentes condiciones de frontera en el extremo
110
110
110
Encontrar:
111
111
111
a. La distribución de temperatura de las varillas si el extremo se
encuentra aislado.
b. El flujo de calor desde las varillas
Solución.-
T 20º C h 9
Q
T 12 D 1, 2
x
L 30
d 2T d 2
hP
T T o m2 0
dx 2 KP dx 2
CF1: o x=0
d
CF2: o x=L
dx
3. Con:
hP
T T y m2
KA
112
112
112
6. Para determinando el flujo de calor desde la aleta para este caso
113
113
113
Qo K A m 0 Tanh mL
7. Calculus
a. Determinación del perfil de temperatura de las 3 aletas
ph
Aleta de vidrio, m2
KA
P h
m Donde: P= Perímetro = D
KA
A= Área Transversal = D 2 / 4
Por tanto:
4h 4 x9
m m 51.40m 1
DK 1.09x1.25x10 2
m = 0.514cm-1
Para la aleta de acero: m 0.07108cm 1
T1 22.11º C
x = 15 cm T2 = 20,044ºC
x = 22,5 cm T3 = 20,0094ºC
x = 30 cm T4 = 20,00004ºC
En la siguiente tabla se presenta la variación de temperatura
para diferentes distancias para los tres tipos de aletas
114
114
114
c. Calculo del flujo de calor desde las aletas
Qo = KA´ m´ Q o ´ Tanh (mL)
D2
Ka 57w / mº C A
4
o Arleta de vidrio
1.25cm
2
m2
Q 1.09w / mº C 120 20 º C 51, 4m Tanh 0.514 30
4 100cm
2
Q 0.688 W
Calcular.
Solución:
110
110
110
2. Se tiene el caso de aletas rectas de sección transversal uniforme.
Considerando el cuarto caso de aletas en donde el extremo, la
conducción de calor es igual a transferencia por convección.
3. Ecuación diferencial gobernante
d 2 hP
m2 0 m2
dx 2 KA
5. Condiciones de frontera
CF1: o x=0
d
CF2: K h x=L
dx
6. Perfil de temperatura.
h
L x
T T cosh m senh m L x
()x
x
Km
0 T0 T cosh mL h Km senh mL
senhmL h coshmL
Km
Qo hPKA b
cosh mL Km senhmL
h
Datos
h 7w / m 2 º C To 30º C K 202w / mº C
T 20º C H 1000mm TL ?
L 30mm Qt ? 3mm
111
111
111
8.1 o To T 30 20 10º C
8.2 L TL T TL 20
112
112
112
hP
m
KA
Donde: A H 0.0031 0.003 m2
Reemplazando
7 x2.003
m 4.81369m 1
0.003x202
Ó sea:
cosh mL L cosh mo 1
senh L L senh(0)0
Reemplazando
TL 20 1
10 1.01044 7 / 202 X 4.81369 X 0.14491
TL 29.8º C
8.5 Cálculo del flujo de calor por aleta Qa
senhmL h coshmL
Qa hKPA o Km
cosh mL
h
Km
senhmL
Se tiene:
113
113
113
hKPA 7 x2.006x0.003x202 2.917096
114
114
114
o 10º C
senh(mL) 0.14491
cosh(mL) 1.01044
( h Km) 7 202x4.81369 0.0071989
Reemplazando:
0.14491 0.0071989x1.01044
Qa 2.917096x10
1.01044 0.0071989x0.14491
Qa 4.3980w
8.5 Calor transferido por los 4 lados con 40 aletas por lado (Qta )
Qta 4.3980x40x4 702.242w
ASA = AtP Aa
2
AtP = área de la pared = b x H = 1 x 0.8 = 0.8m
Aa = área de aletas= H N = 1 x 0.003 x 40 = 0.12m (área
2
Reemplazando
Qt 892.6w
Respuesta:
a. TL 29.8º C
b. Qt 892.6W (calor transferido con aletas)
PROBLEMA Nº 4.3
Solución.-
1) Ecuación Diferencial
d 2T
2
hP
T T 0
dx KA
2) Si T T y m 2 Lp / KA
2
d
Entonces m 2 0
2
dx
3) Condiciones de frontera
CF : 1 q = q0 = T0 - T¥ x0
CF : 2 q = qL = TL - T¥ xL
117
117
117
con CF : 1 0 C1 C 2
mL
CF : 2 L C1e C2 e mL
C 2 0 C1
L 0 e mL C1 e mL e mL
Se determina, C1 y C2
q - q e- mL éq - q e- mL ù
C = L 0
C = q - êL 0 ú
1 2 0
mL
e - e - mL êe mL - e- mL ú
ë û
mL mL
L 0 e mx L 0 e e mx
e 0
e mL e mL e mL e mL
mL mL
L 0 e mx L 0 e e mx
e mx
e
0 e mL e mL
e mL e mL
mx 0 e mx
mL
L 0 e mx
mL mL e e
e
e
æ öæe - e ö÷
L
e- mL ÷ e- mx
mx - mx
q q
= çç - ÷çç ÷+
- mL
q çq mL
÷çe - e ÷
0 è 0 øè ø
118
118
118
Perfil de temperatura
e mx
hiperbólicas, para su determinación se realiza de la siguiente forma
mL e
mx
mx
L 0e 0e
e mL e mL
119
119
119
L e mx e mx 0 e mLmx e mLmx
e e
mL mL
qLsenh mx + q0senh m (L - x )
q=
senh mL
…….. Perfil de temperaturas
Diferenciando el perfil
d L cosh mx. m 0 cosh m L x m
dx senh mL
dq m q cosh mx - q0 cosh m (L - x )ù
é
= ëL û
dx senhmL
dq m [qL - q0 cosh mL ]
=
dx x= 0 senh mL
æq cosh mL - q Lö÷
Qa = KAm ç 0
senh mL ø÷
Diferenciando el perfil
d cosh mx . m 0 cosh m L x m
L
dx x L senh mL
120
120
120
d m L cosh mL 0
senh mL
dx x L
PROBLEMA Nº 4.4
Solución.-
1. Diagrama de la aleta
d 2
m2 0 ………………….. …(1)
dx 2
C.F . 2
2.5 Despejando:
emL emL
c 2 1 ó c1 2 1
e mL e mL
1
2 sen h mL
e mL e mL
c 2 1 ó c2 1 2 1
2 1 mL mL
2 senhmL e e
x
e mL
e mL
123
123
123
e mx e mx
2
1 e mL x e mL x
x
e mL e mL e mL e mL
sen hm L x
1
sen hmx
x 2 ……………………(5)
sen hmL sen hmL
2
x . senhmx senh m L x
1
1 senhmL
d d
Q KA KA …………( )
neto
dx x0 dx xL
Reemplazando en (α)
KAm 0 cosh mL L
0
Qx0
… ……(6)
senhmL
Reemplazando
cosh mL
Q KAm 0 L ………….(7)
xL senhmL
2.6 Cálculos:
D2
124
124
124
x 0.05
2
A , A 1, 963595x103 m 2
4 4
0 204 25 179 , L 93 25 68
125
125
125
Ph 4h
m , m 5, 497471 m1
KA KD
68
45 1, 9634 10 3179 5, 4974 2, 6976
80, 4545W
179
Qx0
2, 5054
179 68 2, 6976
Q 45 1, 9639 103 5, 4974 0, 86041W
xL 2, 5054
Qneto 79, 57 W
Problema N° 4.5
Una varilla larga pasa a través de la abertura en un horno que tiene una
temperatura del aire de 400ºC y se prensa firmemente en la superficie de un
lingote. Termopares empotrados en la varilla a 25 y 120 mm del lingote
registran temperaturas de 325 y 375ºC, respectivamente ¿Cuál es la
temperatura del lingote?
m = 11,56
: Tb = 300ºC
Problema 4.6
Solución:
1. Se desea hallar:
a. TX1 cuando X=X1
b. TX2 cuando X=X2
c. TL cuando X=L
2. Para el desarrollo del problema se hará el siguiente análisis:
121
121
121
La ecuación diferencial para una aleta de área constante es :
;
valor de k (conductividad térmica) de la varilla de estaño lo ubicamos en
el apéndice A (Tabla A1) del Incropera, para lo cual debemos conocer la
temperatura promedio:
Condiciones de frontera:
CF1: x=0 , = b =
F2: x= L ,
Evaluando las CF
obtenemos el perfil de
temperaturas:
122
122
122
Realizando los cálculos:
1.
2.
3. b =
6. Reemplazando se obtiene:
123
123
123
4.5. Superficies Extendidas de área variable
d 2 1 d
2
m2 0
dx r dr
2. Las condiciones de contorno más sencillas que nos permiten determinar las
dos constantes son, suponer conocida la temperatura en la base y despreciable
el calor convectado en el extremo de la aleta, estas son:
126
126
126
CF:1 r r1 0 (T
) 0 Tf (4.33)
d
CF:2 r r2 0 (4.34)
dr
0 C1 I I mr2 C2 K1 mr2
Operando en forma simultanea, las constantes C1 y C 2 son:
0 K1 mr2
C1 (4.36)
I 0 mr1 K 1 mr2 I 1 mr2 K 0 mr1
0 I 1 mr2
C (4.37)
I mr K mr I mr K mr
2
0 1 1 2 1 2 0 1
El flujo de calor disipado por la aleta puede calcularse mediante el cálculo del
calor que fluye por su base
dT
Q a KA r r A 2r1t
dr 1
127
127
127
d
Q r tK (4.39)
2 1
dr r r1
a
C1 I 0 mr C 2 K 0 mr
Donde:
AS 2 r22 r12
2 2
Q 2h. r2 r1 0 (4.44)
Por tanto:
a
2h. r 2 r 2 128
128
128
2 1 0
129
129
129
2.r1
K1 mr1I 1 mr2 I mr1 K1 mr2
m
(4.45)
a
r 2
r 2
K mr I mr I mr K mr
2 1 0 1 1 2 0 1 1 2
Donde:
131
131
131
temperatura del extremo de la aleta original. En consecuencia se puede
escribir la siguiente equivalencia:
Q hPL TL T f hA TL T f ,
A t
Simplificando, L*
P 2
tanh m L A
P
(4.48)
m L AP
L T= espesor aleta
TL Aislado
Tf . h
L*
P = 2 H + 2 t si H t H=1 P=2
A t
L*
P 2
132
132
132
perfil rectangular, triangular y hiperbólico) y Nº 4.12 (caso de aletas
anulares).
Q Qmax
130
130
130
1/2
L h
3/2
= C
kAp
131
131
131
Figura Nº 4.12. Para determinar la eficiencia en aletas anulares
At: área total incluyendo el área de la base sin aletas y toda la superficie
aleteada, At = Ab + NAf.
132
132
132
0 1
NA f
At
1
f
h At NA f N A f Tb T h A NA f 1 f T
b T
t
hA
NA f
1 T
f
T hA T T
f b 0 t b
t 1
At
133
133
133
Fuente: Alan Chapman, Fundamentos de transferencia de calor, Tercera Edición
, Editorial McGraW Hill, México
134
134
134
Figura N° 4.14 Gráfica para la evaluación de la eficiencia de aletas anulares o
circunferenciales
135
135
135
Tabla N° FUNCIONES DE BESSEL MODIFICADAS DE PRIMERA Y
SEGUNDA ESPECIE
x e x I0 x e x I1 x ex K0 x ex K1 x
0.0 1.0000 0.0000
136
136
136
2.0 0.3085 0.2153 0.8416 1.0335
137
137
137
5.4 0.1762 0.1589 0.5279 0.5749
5.6 0.1728 0.1565 0.5188 0.5633
5.8 0.1696 0.1542 0.5101 0.5525
6.0 0.1666 0.1520 0.2019 0.5422
6.4 0.1611 0.1479 0.4865 0.5232
6.8 0.1561 0.1441 0.4724 0.5060
7.2 0.1515 0.1405 0.4595 0.4905
7.6 0.1473 0.1372 0.4476 0.4762
8.0 0.1434 0.1341 0.4366 0.4631
8.4 0.1398 0.1312 0.4264 0.4511
8.8 0.1365 0.1285 0.4168 0.4399
9.2 0.1334 0.1260 0.4079 0.4295
9.6 0.1305 0.1235 0.3995 0.4198
10.0 0.1278 0.1213 0.3916 0.4108
137
137
137
4.6 Problemas resueltos
Problema N° 4.7
Una aleta anular de aluminio de perfil rectangular se une a un tubo circular que
tienen un diámetro externo de 25mm y a una temperatura superficial de 250°C.
La aleta es de 1mm de espesor y 10mm de longitud y la temperatura y el
coeficiente de convección asociados con el fluido adyacente son 25°C y 25
W/m2K, respectivamente.
SOLUCION:
1. Diagrama de la aleta
t
r2c r2 12.5mm 10mm 0.5mm 23mm 0.023m
2
r2 c
1.84 (i)
r1
138
138
138
t
Lc L 10.5mm 0.0105m
2
Ap Lct 0.0105m x 0.001m 1.05x105 m2
1
1 æ 25 W ö2
3 æ ö
÷
2
3 2 ÷
ç
L 2ç h = (0.0105m) 2 ç m K = 0.15 (ii)
kA ÷ ÷
c
p
÷
ø ç 240 W x1.05x10- 5 m 2 ÷
mK ø÷
0.97
Q ' NQ h 1 Nt 2r1 b
Q ' 2.91KW / m
PROBLEMA Nº 4.8
el calor indicado.
140
140
140
Solución
1) Diagrama de flujo
2T f 1 T f 2 285 24º C
Q
L Ln 33
Ln r2 1 1
1 r1 1 30
0.03 x 1000 50 0.033 x 8
r1 h1 K r2 h2
Q QC2
L QC1 Q K
QC1 h1 2 r1 L T f 1 T1
QK 2 LK T1 T2 / Ln r2 / r1
140
140
140
QC2 h2 x 2 r2 L T2 T f 2
141
141
141
Q
L 100.25 Kcal / h. w
Si son necesarias las aletas, ya que el tubo sin aleta no puede aportar la
energía calorífica necesaria.
Q T f 1 T1 T1 T0 T f 1 T0
L 1
Ln r2 / r1
r
Ln 2
2r1 h1 2K 1
r1
2 r1 h1 2K
Separando T0
Q 1 Lnr2 / r1
T f 1 T0
L 2r h 2K
1 1
Kcal 1 Ln 33 30
T T 460 2,58
h. m 2 x 0,03 x 1000 2 x 50
f1 0
T0 = 82, 42º C
2h 2 x8
4.2 Cálculo: r3 0,066 0.68164 0,682
Kt 50 x 0,003
142
142
142
4.3 Por gráfico la eficiencia (aleta anular)
0.95
143
143
143
4.4 Cálculo de Q h a T T ; a 2 r 2 r 2
max 2 2 0 f 2 3 2
Qmax 8 x 2 0,066 2 0,033 2 82,42 24
Qmax 9.59359
L Na t . N a m
L La 1 tN a
Qt Qta QLa
Qt 460 Kcal / h.m x 1m 460 Kcal / h
Qta 9.11 x N a
QLa 100.25 Kcal / h.m 1 N a .t
5.4 Reemplazando
1 41 x 0,003
S 0,02139 m
41
S 21,39 mm
144
144
144
Problema N° 4.8
SOLUCION:
145
145
145
t
r2c = r2 + = (40mm + 2mm)= 0.042m
2
r2c 0.042m
= = 1.68 ()i
r1 0.025m
t
Lc = L + = 15mm + 2mm = 0.017m
2
- 5 2
A = L t = 0.017m x 0.004m = 6.8x10 m
p c
1
3 æ ö
÷
2
3/2 é ù1/ 2 - 4
L 2ç h = (0,017) ê40 - 5ú = 1.216x10
kA ÷
c
pø 240´ 6,8´ 10
0.97
Q xQmax 2h r2,c r1 b
2 2
Q 2x0.97 x40 W
m K 2 2 2
x 0.042m 0.025m x180C 50W
Q 50W
11.05
hAc,b b 40W 2 x2 0.025m x0.004m x180C
m K
Q ' NQ h 1 Nt 2r1 b
Q ' 125x50 W / m 40 W
m2 K
1 125x0.004x2x0.025mx180C
Q ' 6.82KW / m
146
146
146
Peoblema N° 4.9
147
147
147
Si la pared del tubo esta a 100ºC y el flujo contiguo esta a 25ºC, con un
coeficiente de conveccion de 75 w/m2K. ¿Cuál es la transferencia de calor de
una sola aleta? ¿Cual será la transferencia de calor si la resistencia de
contacto pudiera eliminarse?
Solución.-
1. Diagrama de flujo
Tw T
qf
Rt ,c R f
Rt ,c 2 x10 4
Donde: Rt ,c 1.06 K / W
Ab 2 0.015 m 0.002 m
b b b 1
Además: R f
n f qmax n f h Af b 2 h nf (r 2 2,c r1 )
2
qf
Lc L t / 2 0.016m
r2,c / r1 0.031/0.015=2.07
Ap Lc t 3.2x10 5 m 2
148
148
148
(h / kAp)1/ 2 0.0163 / 2 x(75 / 240 3.2 105 )1/ 2 0.20
3/ 2
Lc
149
149
149
Nf = 0.94
1
Rf 3.07 K / W
2 (75W / m .K )(0.94)((0.031m) 2 (0.015m) 2 )
2
PROBLEMA Nº 4.10
El calor es transferido a través de una pared plana de espesor x 1.25 , de
Se decide que el intercambio de calor entre ambos fluidos usar aletas (rectas)
de perfil rectangular de longitud L = 2.5 cm de espesor w = 0.16 cm y
espaciadas 1.25cm sobre sus centros, la conductividad es la misma del
material de la pared asumiendo que el flujo es unidimensional, encontrar el flujo
de transferencia de calor con unidad de área de la superficie primaria de la
pared, si esta a) no tiene aleta. b) si se adicionan aletas al lado derecho.
Solución.-
1. Diagrama de Flujo
Figura Nº 4.19 Corte transversal de una pared con aletas y sin aletas
150
150
150
Fuente: Elaboración propia, (Alberto Emilio Panana Girio)
151
151
151
2. Para la pared sus aletas por área primaria Ap . El flujo de calor es.
q T f 1 T f 4
H 1 (profundidad unitaria)
Ap 1 x 1
h k h2
q 120 20
Ap 1 0.0125 1
2364.9 W / m2
450 200 25
3
152
152
152
3.5 Cálculo de relación área total At a área sin aleta Ap
At 2L
5
Ap
q 120 20
Ap 1 0.0125 1
459 20 5 x 0.9742 x 25
2 x 450
mL 0.025 1.3258
200 x 0.0016
Af 2L
0.8256
At 2L
Af 2L
5
At
153
153
153
Af
2¡ 1
A
1 1 0.82561 0.6549
t
0.7151
2¡
PROBLEMA N° 4.11
Una aleta circular con una sección transversal rectangular de 3,7 cm. De
diámetro externo y 0,3 cm. de espesor, rodea un tubo de 2,5 cm. de diámetro.
La aleta está construida de acero inoxidable de K = 14,4 W/m.K. El aire que
sopla la aleta produce un coeficiente de Transferencia de Calor de 28,4 w/m2 K.
Si las temperaturas de la base de la aleta y el aire están a 260ºC y 38ºC,
respectivamente.
Calcule:
154
154
154
Fuente: Elaboración propia, (Alberto Emilio Panana Girio)
151
151
151
Solución.-
d 2 1d
m 2 0
dr 2 1dr
2. Condiciones de frontera
C.F : 1 0 ó T T0 ; r r1
d dT
r r2
C.F : 2 0 ó 0 ;
dr dr
dI 0 ()m()r dK 0 mr
C mI ()mr , C m(K
) mr
1 1 2 1
dr dr
150
150
150
I1 ()m()r02 C1 K1 mr2 C2
K1 ()mr2 0
C1
I 0 ()m()r1()()K1 mr2 I1 mr2 K0 mr1
151
151
151
I1 ()mr2 0
C2
I 0 ()m()r1()()K1 mr2 I1 mr2 K0 mr1
K ()m()r(2)() I 0 mr I1 mr2 K0 mr
1
0 I 0 ()m()r1()() K1 mr2 I1 mr2 K0 mr1
7. Cálculos:
2h 2 x28, 4
m 36, 260m 1
Ke 14, 4x0.003
t t 0, 3
r 1, 25cm 0, 0125m ; L L ()r0,6r 0, 75
0, 0075cm m
1 c 2 1
2 2 2
Por Tanto:
mr1 36, 260375 0.0125 0, 45325
(260 38)1.0
1126659 (I)0.m
38r4451 () K 0 mr
T T f
1.45835
T - Tf = 153.9419l 0 ()m
5r8.52+37 () K 0 mr
152
152
152
10. La razón de transferencia de calor de la aleta
Qaleta Qmax ; gráficamente (1)
153
153
153
h
1/ 2
r2 c 3/ 2
...(1)......(2) L
c
r1 KAp
t
r R L L L L r r
2c 1 c c 2 1
2
4
Lc 0, 75cm 75x10 m r2c 2cm 2 x10 2 m
r2 c 2 x102
1,6
r1 1,25x102
1/ 2 1/ 2
h = (75x104 )03 /,219229 28, 4
L 3/ 2
KAp
c
4
0, 225x10 x14, 4
De la grafica: η = 0,985
Qmax h Asup (T
) 0 Tf
A 2 ()r 2
r2
sup 2c 1
Qaleta Qmax
Problema N° 4.12
155
155
155
a) El calor disipado en cada metro de longitud de tubería sin aletas.
b) El calor disipado en cada metro de longitud de tubería con aletas.
c) La temperatura en el extremo aislado de la aleta.
d) El aumento en porcentaje de calor disipado por el hecho de colocar las
aletas.
Solución.-
1. Diagrama de flujo
H a 1 5x103
Na 182 aletas
w a 0,5x103 5x10
156
156
156
3.1 Calculando n` (eficiencia global)
157
157
157
Af Na
1 (1 ) . ……….(5)
At
K ()m()r(2)() I 0 mr I1 mr2 K0 mr
1
0 I 0 ()m()r1()()K1 mr2 I1 mr2 K0 mr1
158
158
158
2 h 60 2 1
m 197 0, 005 34, 90378m
k t
159
159
159
4.4 Cálculo de las funciones de Bessel , los cuales se dan en la siguiente tabla:
Con: mr1 = 0,69807 y mr2 = 1,745189
L I1
1
K ()m()r(2)() I 0 mr2 mr2 K0 mr2
L (T) L T f cuando r r2
0 (T)50 0 T f C r r1
4.6 Calculo del aumento de calor disipado con el tubo con aletas
PROBLEMAS Nº 13
161
161
161
b) Material con el que se debe construir la aleta triangular, tomando como
referencia su conductividad térmica para que en las condiciones anteriores
tenga la misma efectividad que la encontrada para la alta rectangular
Solución
1. Diagrama de flujo
mL KA
P 2 H t H t P 2H A H xt
162
162
162
m 8.2462 m1
163
163
163
3. Reemplazando : Tanh 1.644924
tanh 8.2462m1 0.2m 0.92875
1. 0.563
8.2462m1 0.2 m 1.64924
Ph
Nota: 1 Cálculo de m
KA
A 1 x 0.0125 0.0125m 2
2.025x17
m 68.85 8.29759
0.0125x40
tanh 8.29759m 1 0.2 m 0.93015223
0.56049
8.29759 0.2 1.659518
56 %
Reemplazando
17 x 2.025x9.9125x40
senh 8.29759 x 0.2
17 8.29759 x 40 cosh 8.29759 x0.2
cos h 8.29759 x 0.2
17 8.29759
x 40 s en h 8.29759 x0.2
17 x 0.4
165
165
165
2f
Parámetro: L 0.2 x 8.2462 1.649 1.65
Kt
t 0.0125
1) Longitud corregida, Lc L 0.2 Lc 0.20625 m
2 2
2) Area del perfil Ap Lc.t 0.20625 m 0.0125 2.578 103 m 2
17 / 40 x 2.578x103
1/ 2
0.20625
3/ 2
1.20265
1/ 2
t 2
Área de Perfil: Af 2H L2
2
A x AI Htx
AI Ht, entonces A x
X L L
Ht
dA x
L
166
166
166
2) Área superficial de la aleta : S x
167
167
167
2
S 2bH , b x2 C*
x
*
t/2 C t / 2X
Donde C * se determina: , donde C *
L x L
t / s X
2 2
t
O sea b x 2 X 1
L 2L
2
t
Por tanto: S x 2Hx 1
2L
2
t
Si f 1
2L
dS x
S x 2Hx f y 2Hf
dx
Htx d 2 Ht d
K . h 2Hf 0
L dx 2 L dx
2hLf
Si n además si L b f 1
Kb
dS x
S x 2 XH 2 H=1
dx
169
169
169
dA x dS x
Por tanto si se reemplaza A x
dx dx
Htx d 2 Ht d
K . h 2 0
L dx 2 L dx
2hL
Donde : n m L
kt
I 0 2n L
161
161
161
6.4) El calor disipado al exterior por la aleta longitudinal de
anchura unitaria será igual al que penetra por
conducción por su base, por lo que:
d
Qa KA
dx xL
A tH H 1 A t
d I 0 ax
aI ax
1
dx
d k0 ax
ak ax
1
dx
d 0
x I1 2n1/ 2 2nx1/ 2
dx xL I 0 2n L xL
I0 0
x
2x
2n
I
1/ 2 1 2nx1/ 2 xL
2n L
x
2n
I1 2n L
I0 0 2 L
2n L
Reemplazando en Qa
Q Kt
2n
I1 2n L
a 0
2 L I 0 2n L
2hL 2h
n m L m
kt kt
Cálculo de la Eficiencia
QQ
0 Qmax h2 H x d T0 Tx
Qmax
162
162
162
2
t
d L2 como L t H =1
2
dL
Reemplazamos
163
163
163
I 2n L
kt
2n x 1
0
a
0
2 L I 2n L
h2L0
2hL
donde n
kt
1
x
I1 2n L ó n 2
2hL
a
n L
I 0 2n L kt
Si : n m L n L m L. L n L mL
I1 2n L I1 3.2985 5.195
I 0 2n L I 0 3.2985 0.258
Reemplazando
1 5.195
0 x 0.503
1.64925 6.258
164
164
164
i. Se observa que el rendimiento de las aletas rectangulares es superior al de
la aleta triangular.
j.
a. La Determinación de la Eficiencia mediante gráfico Fig. 217
(Aleta triangular de perfil rectangular
165
165
165
*) Cálculo del perímetro
2 x17
L 2h 0.2 1.649
kt 40 x 0.0125
n L
I 0 2n L
1 I1 2 N
a 0.58
N I 0 2N
166
166
166
4.7 Problemas Propuestos
Problema 4.7.1
Problema 4.7.2
Problema 4.7.3
Problema 4.7.4
167
167
167
La cara superior de una viga de forma (I) de 12 pulgadas se
mantiene a una temperatura de 500 °F, mientras que la inferior está
a 200 °F. El espesor de la estructura es de 1/2 pulgada y a lo largo de
la viga sopla aire a 500 °F se tal suerte que h = 7 Btu/ h.pie2.°F.
Suponiendo que la conductividad térmica del acero es constante e
igual a 25 Btu/h.pie.°F, determine:
Problema 4.7.5
Problema 4.7.6
168
168
168
temperatura de 25 °C con un coeficiente de transferencia de calor h 0 =
200W/m2
Problema 4.7.7
(b) ¿Cuál es el calor disipado por unidad de ancho por una sola aleta?
169
169
169
V. CONDUCCIÓN DE CALOR EN RÉGIMEN ESTACIONARIO EN FUNCIÓN
DE DOS O MÁS VARIABLES
Figura N° 5.1 Placa plana en dos dimensiones, con temperaturas en los borde
iguales
170
170
170
T = X(x) Y (y) ( 5.2)
Como cada miembro de esta ecuación depende solo de una variable, los
dos miembros tienen que ser iguales a una constante ( 2 )
1 2 X 1 2Y
2 (5.4)
X x 2
Y y 2
T 1 Senh y 2Cosh y 3 Sen x 4Cos x (5.7)
Se tiene que (λ) y las (β), son constantes que hay que determinar
mediante las condiciones de contorno.
171
171
171
Entonces la ecuación general (5.7) se reduce a:
172
172
172
T= 1
3 (5.10)
Además : β = β1 β 3 (5.11)
Para que esta ecuación se cumpla para todos los valores de (y), es
necesario que:
- Sen
- Que se satisface 2
λ=0, ,
a a
n
- Y en general por: , siendo n 0,1, 2....
n
a
Para cada valor de (n) se obtiene un valor de (λ) que proporciona una
solución diferente de la ecuación
obteniéndose
T= nSenh n y n (5.15)
n 1
173
173
173
πn
Con: n x
a
Sen i dx
j con i j
0
T=f(x)= n n n 1 1 1 n enh n n
n=0
a a a
0 0 0
Por definición de ortogonalidad se hacen cero (0), todas las integrales del
segundo miembro, menos la correspondiente al coeficiente βn por lo que:
a
f(x)Sen n dx a
2
β Senh λ b 0
a 0
f(x)Sen(λ x)dx (5.17)
n n a n
Sen dx
2
n
0
ny (5.18)
Senh a
171
171
171
5.1.2 Determinación de la temperatura para un sistema bidimensional (placa
rectangular o cuadrada), cuando una frontera se encuentra a una
temperatura uniforme diferente a las temperaturas de las otras fronteras
nx nx
1 1
a a
T0 a
f(x) Sen dx T0 Sen dx
n
(5.19)
0
a 0
a n
a T0
nb n
T0 n=1
Senh n=1,3,..
nb n
Senh
a a
(5.20)
172
172
172
Figura N° 5.3 Placa plana en dos dimensiones, mostrando las isotermas,con
borde inferior caliente
T0
4
n=1,3,..
nb n
(5.21)
Senh
a
nx ny
Senh Sen
T b b
4 (5.22)
T0 n=1,3,..
na n
Senh
b
T0
4
n=1,3,..
na n
(5.23)
Senh
b
174
174
174
b2) Para el caso particular del calor transmitido a través de la superficie x =
0, por unidad de longitud perpendicular al plano (x, y) se tiene:
b
T(,x )y
Q x 0 K x 0 dy
y 0
X
ny
b
2K Senh
a n a nx
Q f(x)Sen dx dy
X 0 a
nb a a
n=1
y 0 Senh 0 (5.24)
a
nb
Cosh 1
a
nx
2K
Q X 0
a n=1
a
nb f(x) Sen
0 a
dx
Senh
a
Problema N° 5.1
(1, 0,5)
T2 -T1
n=1,3,5.. n nW L
Senh
L
Diagrama de flujo:
175
175
175
Figura N° 5.4 Placa plana rectangular (bi-dimensional), con tres bordes a
temperaruras iguales, y la superior a otra temperatura
Solución.-
(1)0, 5
Senh
1
1 1
Sen
(1)(1) 2
(1)(1)
1 2 Senh
T(1, 0,5)1-T
2 2
(1, 0,5)
T -T (3)0, 5 (5)0, 5
2 1 3 5
1 1 (3)(1)( Senh 1 1 Senh
5)(1) 2 2
Sen Sen
3 2 (3)(1)(5)(1) 5 2
Senh Senh
2 2
T(1, 0,5)1-T
(1, 0,5) 0, 4457757
T2 -T1
T 0, 4457757(150 50)50 94, 5757 C
(3)0, 5
(3)(1) Senh 2
Sen
2 (3)(1)
Senh
2
T(1, 0,5)1-T
(1, 0,5) 0, 440742
T2 -T1
T 0, 440742(150 50)50 94, 0742 C
nx
Sen , el valor de = 180 °
L
ny nW
Senh , y Senh el valor de = 3,1416
L L
177
177
177
Figura N° 5.5 Conducción unidimensional (Método numérico)
Q f a Q1a Qgen 0
T1 Ta x
hA T T kA q A 0
x
f a o
2
i 4
Q
i 1
i 0 0
Q
i 1
i 0 0 T1 +T2 +T3 +T4 4T0 0
A todos los nodos que se encuentran situados en toda la periferia del sólido,
hay que hacerles un balance de energía por separado.
181
181
181
Si el sólido esta en contacto con un fluido a Tf, Fig. N° 5.7, con un
coeficiente de transmisión de calor por convección (h), se asigna a cada
nodo de este tipo la mitad de la superficie que a cualquier otro nodo
interior. El nodo (0) puede intercambiar calor por conducción con tres
nodos continuos, y transferir calor por convección al fluido
T1 To T2 T0 T3 T0 T4 T0
k Δy Δx Δy o
Δx
5.32)
180
180
180
2- Si la red es cuadrada: Δx = y
181
181
181
x 2
T1 T2 +T3 T4 4T0 qo 0
K
x 2 (5.33)
T1 T2 +T3 T4 qo
T0 K
4
Considerar una cámara cúbica de paredes gruesas y aristas externas 2a, 2b,
2c, la octava parte se representa en la figura N° 5.8
182
182
182
Figura Nº 5.9 Nodo Interno ( i, j, k ) de una cámara cúbica (con generación
interna de calor q0 W / m
3
, en estado estacionario).
Fuente: Elaboración propia, Ing Alberto Emilio Panana Girio
Nodo 1 i 1, j, k Nodo 5 i, j, k 1
Nodo 3 i, j 1, k
T T T T T T
k yz 1 0 k yz 2 0 k xz 3 0
x x y
T4 T0 T5 T0 T6 T0
k xz k xy k xy (5.34)
y z z
q0 xyz 0
6. Simplificando k
x 2
T1 T2 T3 T4 T5 T6 6T0 q0 0
183
183
183
(5.35)
184
184
184
7. Si en la cámara no hay generación interna de calor
T T T T T T
T 1 2 3 4 5 6 (5.36)
0
6
PROBLEMA Nº 1
Solución.-
185
185
185
2Q21 Qc Q51 0
2K
Y
L
T2 T1
hXL T T K XL
T 5 T1
0
A F 1 A
2 X Y
Donde: X Y
K A T2 T1 hX TF T1 K A T5 T1 0
2K A hX T1 K AT2 K AT5 hXTF ....... 1
Ecuación para el Nodo 2:
186
186
186
Q 34 Qc Q84 0
Y T3 T4 h X Y X T8 T4 0
K L L T T h L T T K L
A F 4 F 4 A
2 X 2 2 2 Y
Donde: X Y , T8 T3
187
187
187
KA K
T3 T4 hX TF T4 A T3 T4 0
2 2
K AT3 K A hX T4 hXTF ....... 4
B A A B
X
2 2 X Y Y
Donde: X Y
2K B T6 T5 K A T6 T5 K A T1 T5 K B T9 T5 0
K AT1 2K A K B T5 K A K B T6 K BT9 ....... 5
Ecuación para el Nodo 6:
B A A A
2 X 2 X Y 2 X
Y T7 T6 K T10 T6
K L XL 0
B B
2 X Y
Donde: X Y
K K K K
T T T T K T T T T T T K T T
B
5 6 A
5 6 A 2 6 A
7 6 B
7 6 B 10 6 0
2 2 2 2
K K K K
K AT2 A A B 6 A B T7 K BT10
B T 2 K K T
5 ....... 6
2 2
188
188
188
B A A A
2 X 2 X Y X
K
X
L
T11 T7 K X
L
T11 T7 0
A B
2 Y 2 Y
Donde: X Y , T11 T6
189
189
189
K K K K
T T T T K T T K T T T T T T
B
6 7 A
6 7 A 3 7 A 8 7 A
6 7 B
6 7 0
2 2 2 2
2K AT3 K A K B T6 3 K A K B T7 0 ....... 7
PROBLEMA Nº 2
190
190
190
Figura Nº 5.11 Sistema bidimensional para evaluar la ecuación nodal de
temperaturas para los nodos: (m,n), (m+1, n), (m-1, n-1),
(sistema bidimensional en estado estacionario).
Solución.-
Qm1,m,n Q Qgen Qc 0
x y
y Tm 1 Tm,n
q 0
x 2 2
K q 1 hL T Tm,n 0
2 x 2
1.2 Cálculo de L:
2 2 2 2
x y x y
L Si x y
2 2 2 2
x 2 1
L 2 x
4 2
2
k
T 1 T , qx q0 x
hx 1 / 2 T T , 0
m mn mn
2 2 4
191
191
191
1.3 Simplificando y multiplicando por (2/k)
192
192
192
x q 0x 2
T 1 T ,
q
h
2 x 1 / 2 T T , 0
m mn mn
k 4k k
2h
qx q 0x 2 2h
T 1 x 1 / 2 T x 1 / 2 T 0
m 1 m,n
k k 4k k
y Tm,n Tm1,n
y Tm 2,n Tm1,n Tm1,n1 Tm1,n y
k k k x qx1 q0 x 0
2 x 2 x y 2
2 x 2
Tm,n 4Tm 1,n Tm 2,n 2Tm1,n 1 q x q0 0
k k
3) Para el nodo (m+1,n-1)
Qm1,m1,n1 Qm 2,n1m1,n1 Qc Qgen 0
H x 2 y 2 x 2
Multiplicando 2/k
q x 2
2 Tm1,n Tm1,n1 2 Tm 2,n1 Tm1,n1 2h
x 2 T Tm1,n1 0 0
k k
2h q x 2hx 2 2
2T 2T 4 x 2 T 0 T 0
m 1,n m 2,n 1 m 1,n 1
k k k
193
193
193
PROBLEMA Nº 3
1) Balance de Energía ; d = 1
Qm1, n m,n Qm, n1 Qc1 Qc2 Qgen 0
k y .1 Tm 1 , n Tm,n x Tm, n 1 Tm,n hy
T T h x T T q x y
. .1
f m,n f m,n 0
2 x 2 y 2 2 2 2
2
Multiplicando por
k
x x x 2
Tm1 Tm,n1 2Tm,n 2h T f 2h Tm,n q0 0
k k 4
2
2hx .T f x
q0 x
T T h T
m,n1 m1,n 2 1 m,n 0
k k 4k
194
194
194
PROBLEMA Nº 4
Hallar la ecuación nodal para nodo situado en una esquina interna con
convección y generación, en un sistema bidimensional., para una placa
rectangular y placa cuadrada
Solución
1) Balance de energía d = 1
k y
Tm 1,n Tm,n k y Tm 1,n Tm,n k x Tm,n 1 Tm,n k x Tm,n 1 Tm,n
x 2 x y 2 y
2 2 2 2 2
2
Si el sistema es una placa cuadrada x y , multiplicando por
k
x x 2 x 2
2Tm1,n Tm1,n 2Tm,n 1 Tm,n 1 6Tm,n 2h Tm,n 2h x T f q0
k k
2 4
x 1,
2 T T T T 2h T m
n
m
190
190
190
m 1,n m,n 1 f
hx x 2
2 3 Tm,n q0 .3
k k 4
191
191
191
Problema N° 5
Solución.-
3. Ahora : si x y z ,
192
192
192
x x x x x 2
k T1 T0 k T2 T0 k T3 T0 k T4 T0 h T f T0
2 4 4 2 2
2 3
x
h. T f T0 q0 x 0
2 4
4
4. Si multiplicamos todo por obtenemos
k.x
h h x 2
2T1 2T0 T2 T0 T3 T0 2T4 2T0 4 xT f 4 xT0 q0 0
k k k
4h 4hx q0 2
2T T T 2T T x T 6 x 0
1 2 3 4 f 0
k k k
PROBLEMA Nº 6
Tm 1 T Tm1 0
k 2kA k 2kA
193
193
193
Donde: (m) denota el nodo que se encuentra en la punta de la aleta.
kA
El flujo de calor en la base es: q x 0 Tm 1 Tm
x
194
194
194
Donde Tm en la temperatura de base y Tm1 es la temperatura en el primer
incremento.
Solución.-
1. Diagrama
195
195
195
Q(1m)(1) m Q m m Q
convec
T T T ,T
kA m m 1 kA m 1 m hPx T T
m
x x
hPx 2 hPx 2
Tm 2 T Tm 1 Tm 1 0
kA kA
T1 300º C
5 6
x
2
x
Balance de Energía
Q56 Q Qcon 6 0
con 6
196
196
196
T T hPx
kA 5 6 T T6 hA T T6 0
x 2
5. Ecuaciones Nodales
Nodo 1: T1 300º C
D2 P D 4
Nodo 2: con P D A
4 A 2 D
D
4
h4x 2 h4x 2
T2 2 T T1 T3 0
kD kD
Reemplazando valores
Nodo 3
4hx 2 4hx 2
T3 2 T T2 T4 0
kD kD
T3 2 T2 T4 0
228 x 0.025 228 x 0.025
197
197
197
20107 T3 0.408 T2 T4 0
Nodo 4
198
198
198
4hx 2 4hx 2
T2 2 T T3 T5 0
kD kD
4 17 0.03 4 17 0.03 x 38
2 2
T4 2 T3 T5 0
228 x 0.025 228 x 0.025
20107 T4 0.408 T3 T5 0
Nodo 5
4hx 2 4 h x 2
T5 2 T T 4 T6 0
kD k D
4 1 7 0 .0 3 4 1 7 0 .0 3 x 3 8
2 2
T5 2 T 4 T6 0
2 2 8 x 0 .0 2 5 2 2 8 x 0 .0 2 5
20107 T5 0.408 T4 T6 0
Nodo 6
2hx 2 hx 2 x 2h x 2
T6 1 h T T5 0
kD k k kD
17 x 2 x0.032 17 x 0.03 17 x 0.03 2 x17 x 0.032
T6 1 x 38 T5 0
228x 0.025 228 228 228 x 0.025
1.007605 T6 0.289 T5 0
T2 228.049 T3 278.779 T4 272.083
Rpta:
T5 267.89 T6 266.157
PROBLEMA 7
199
199
199
Considere la transferencia de calor bidimensional en estado estacionario en
una barra sólida larga de sección transversal cuadrada en la cual se genera
calor de manera uniforme con una velocidad de q0 0.19 x105 BTU / h. pie 2 .
Determine:
200
200
200
a) Las temperaturas en los nueve nodos
b) La velocidad de la pérdida de calor desde la barra a través de una
sección de 1 pie de largo.
Solución.-
1.- Datos
BTU BTU
q0 0.19x105 2
h 16 T 70º F
h. pie h pieº F
BTU
h 7, 9 x y 0.2 pies
h pie º F
q x 2
T4 T5 T6 T5 T2 T5 T8 T5 0 0
k
0.19x105 x 0.22
T2 T4 4T5 T6 T8 0
16
201
201
201
Para el nodo 2:
Q12 Q32 Q52 Q62 Qgen 0
Para el nodo 8:
Para el nodo 6:
Para el nodo 4:
Para el nodo 1
Para el nodo 3:
Para el nodo 7:
202
202
202
T4 2.1975T7 T8 37.575 0 ................... 8
203
203
203
Para el nodo 9:
Resolviendo
T5 328, 0394º F
T1 304, 8482º F
T6 316,1644º F
T2 316,1644º F
T7 304, 8482º F
T3 304, 8482º F
T8 316,1644º F
T4 316,1644º F
T9 304, 8482º F
Q hyT
y y
Q h T3 T y T6 T T9 T8
2 2
1
304, 8482 70 316,1644 70 304, 8482 70
1
Q hy
2 2
BTU
Q 7, 9 0, 2 pie 481, 0126 º F
h pieº F
BTU
Q 759, 9999
h pie
PROBLEMA N° 8
Solución:
1. Balance de energía
204
204
204
2. Reemplazando la ecuación de Fourier en cada una de las cantidades de
calor por conducción, el calor convectivo y generado y simplificando se
tiene,
2hD xT q D x2 2hD x
f o
T +T +T +T +2T + + - (6 + )T
0 =
1 2 3 4 5 o
k k k
Problema N°1
200
200
200
Problema N° 2
Determinar la ecuación nodal, para un nodo identificado por (i, j, k), en tres
dimensiones, en estado transitorio, con generación de calor interno por unidad de
volumen y por unidad de tiempo (qo), situado en un lado de una cámara cúbica o de
un horno, de conductividad térmica (k), expuesto a un medio ambiente de
temperatura (Tf) y coeficiente de transferencia de calor (h). Las celdillas para un
sistema tridimensional x, y, z, son de amplitudes Δx = Δy = Δz, el material tiene una
densidad (ρ) , conductividad térmica (k), y difusividad térmica (α).
Problema N° 3
Problema N° 4
201
201
201
b) La superficie frontera está sometida a un flujo de calor entrante de valor
constante,
202
202
202
vecinos, y de las características térmicas y geométricas del nodo. Se ha de
tener en cuenta la energía almacenada en cada nodo de la celdilla, en el
tiempo considerado, la cual se puede expresar como la variación de la energía
interna del mismo.
U 0 T T t t T0t
m.Cp 0 .A.x.Cp 0
t t t
luego : (6.3)
U T Tt t
T T
t t
T t t
T t
Q Q 0
kA 1 0
+kA 2 0
.A.x.Cp 0 0
t x x t
1o 2o
Despejando, T0t t
t t k t t k t
t t t t
T0 T0 T T 2T F
x 2 .Cpx 2
.Cp.x 2 1 2 0 0
&.4)
F0 T1 T2 T0 1 2F0
t t t t t
T0 203
203
203
Figura N° 6.1 Nodo interno (0) en una placa plana, sistema unidimensional
transitorio
x t
F 0
2
t
F0 2
x
2F0 T1 BiTF 1 2F0 2F0 Bi T0
t t t t
T0 ; (6.7)
Bi h.x
k
204
204
204
Figura N° 6.2 Nodo en frontera convectiva de un plana, sistema unidimensional
transitorio
h.L
Bi 0,1
k
205
205
205
Realizando un balance de energía al sistema mostrado en la figura N° 6.3, que
se encuentra a T=T(t) en el instante (t), se tiene que la variación de su energía
interna en ese instante es igual a la energía que es transferida al fluido que le
rodea en dicho instante:
T h A T (t)- T f
h A T (t)- T f
T
Q - V Cp - (6.8)
t t V Cp
t V Cp t V Cp V Cp
(6.9)
hA
t hL t
(t) e V Cp
e Bi Fo
; Bi ; Fo
o o
k L2
f o
q t dt h eA T T
t t
Bi Fo
Q(t)() dt
o f t 0
t 0
Bi Fo
V Cp
=h A T T 1-e = h A T t
Bi Fo 1-e
T (6.10)
o f o f
hA Bi Fo
Como:
Q(t)
Q V Cp T T 1- e Bi Fo (6.11)
0 o f
Q0
Q(t)
Fracción de pérdida de energía
Q0
208
208
208
C T T f TC T f T T f
(6.15)
0 0 C TC T f T0 T f T0 T f
209
209
209
6.3.2 GRAFICA PARA CALCULAR EL CALOR TRANSFERIDO
Q0 – Q(t).
Bi h.L
Q()t k
(6.20)
h
2
t
Qmax Bi F
2
o
k2
210
210
210
6.5 CONDUCCIÓN TRANSITORIA EN UN CILINDRO INFINITO CON
CONDICIÓN DE CONTORNO DE CONVECCIÓN
211
211
211
Este problema se resuelve análogamente al caso anterior. Si el radio exterior
del cilindro es (R), que se calienta o enfría y que inicialmente para t = 0 en el
intervalo (0 r R) tiene una distribución de temperatura conocida de la forma
f ()r ó () To T f , la ecuación en coordenadas cilíndricas es:
2 1 1
, con T T f , se puede resolver gráficamente
r 2 r r t
1
T T Bi
f
(6.21)
.t
T T
0 0 f Fo 2
R
r
T T R
(,r )t(. ) rt f
(6.22)
1 k
T T
C C f
Bi h.R
212
212
212
La relación de calor a dimensional se determina el la gráfica de la figura N°
4.14 c, conocida como gráfica de Gröber, interceptando los parámetros:
k2
213
213
213
Bi h.R
Q()t k
(6.24)
2
Qmax h t
Bi 2 Fo
k2
2 2 1
, con T T f (6.25)
r 2 r r t
1
T T Bi
f
(6.26)
.t
T T
0 0 f Fo 2
R
r
T R T
(,r )t(. )
rt f
(6.27)
1 k
T T
C C f
Bi h.R
214
214
214
c. La expresión del flujo de calor a dimensional, Q/Qo se conoce como
fracción de energía pérdida, y es la pérdida real de energía en el tiempo (t),
dividida entre la pérdida total necesaria para alcanzar la temperatura del
medio ambiente.
215
215
215
La pérdida total de calor necesaria para alcanzar la temperatura del medio
ambiente se halla mediante la relación. También es conocida como cantidad
de:
Qmax m Cp (T Ti (
) V ) Cp T Ti (6.28)
k2
Bi h.R
Q()t k
(6.30)
h
2
t
Qmax Bi 2 Fo
2
k
Q = (,x)t f
= Ae
- l1t
(l x ) t > (6.33)
i
1
f cos .
(,x)t1 pared
T - T L
0,2
T T 2
r
Q = = Ae J çl t > (6.36)
(,r )t1cilindro 0 - l1t
(,r )t - f æ ö÷ 0,2
1
T- T çè . r ÷
ø
i f 0
217
217
217
- Donde las constantes A1 y λ1 son funciones sólo del número de Bi y
en la Tabla (6.1) se de una lista de sus valores con respecto al
número de Bi.
218
218
218
-
CALCULO DE TEMPERATURAS EN EL CENTRO, PARA UN
CILINDRO
- De la ecuación anterior (para determinar la temperatura local), se tiene
que: J0(0) = 1, la temperatura en el centro se halla, mediante la
ecuación
T (0,t )- T f 2
Q
= = A e- l 1 t
(,O)1
t cilindro (6.37)
T- T
i f
219
219
219
- En el límite de (sen x)/x, es igual a 1, de la ecuación anterior (para
determinar la temperatura local), la temperatura en el centro del
cilindro, se halla mediante la ecuación:
220
220
220
T (0,t )- T f 2
Q
= = A e- l 1 t
(,O)1
t esfera (6.40)
T- T
i f
- Pared Plana
æ Qt ö÷ senl 1
1
ç ÷ = - q (6.41)
0,pared
çèQ ø÷ l
max pared 1
- Cilindro sólido
æ Qt ö
÷ J1()l 1
1 2
ç ÷ = - q (6.42)
0,cilindro
èçQ ø÷
max cilindro
l 1
- Esfera sólida
æ Qt ö
÷ Senl 1 - l 1Cosl 1
1 3
ç ÷ = - q (6.43)
0,esfera 3
èçQ ø÷
max esfera
l
1
222
222
222
Los problemas de conducción transitoria estudiados se limitan a
configuraciones especiales como son la placa, el cilindro y la esfera, con
diversas situaciones de contorno. Estas formas se han escogido para
asegurarnos de que la temperatura del sólido depende solo de la coordenada
espacial en el tiempo. En ciertas aplicaciones se tiene que considerar la
conducción transitoria en función de más de una dimensión espacial.
¶ 2T ¶ 2T 1 ¶T
+ = (6.45)
2 2
¶x ¶y a ¶t
æ ö
dT hcT ÷
çen, x = 0,en, x = A, = ± ÷
dx
k ÷
Para, t > 0, ç ÷ (6.46)
ç
dT hT
çen, y = 0,en, y = B, = ± c ÷
èç dy k ø÷
223
223
223
Para t = 0, T = f (x y) = f1 (x) f2 (y)
224
224
224
Y si esto es posible, la solución de la ecuación (1), con las condiciones
indicadas, se puede expresar como el producto de dos soluciones
transitorias unidimensionales.
T = Tx (x, t) Ty (y, t)
1. Estas ecuaciones para placa finita son idénticas a las que regulan la
conducción transitoria de calor en una placa plana infinita. Por tanto, la
solución al problema de conducción transitoria del calor en la barra
rectangular se obtiene como el producto de las soluciones para dos placas
infinitas cuya intersección forma la barra en cuestión.
4. Para el cilindro circular como el producto de las soluciones para una placa
infinita y para un cilindro circular de longitud infinita.
225
225
225
Si por ejemplo, se desea determinar la temperatura en el punto P del cilindro de
longitud finita que se muestra en la Fig, 6.4, dicho punto vendrá localizado por
dos coordenadas (x, r), siendo x una coordenada axial medida desde el centro
del cilindro y r su posición radial.
En las graficas que se representan en las Fig. 6.6 y Fig 6.7, se hace un
resumen de las soluciones para sistemas bidimensionales y tridimensionales,
mediante gráficos de Heisler ( para sistemas unidimensionales), utilizando las
siguientes relaciones:
F (x,t )
S (x ) = (solido semi ¥ )
F0
F (x,t )
P (x ) = (placa ¥ ) (6.49)
F0
F (r ,t )
C (x ) = (cilindro ¥ )
F0
Para hallar el calor total, se puede utilizar una expresión debida a Langston, de
la forma:
220
220
220
Q = Q r Cp V (T0 - TF ) (6.50)
Qt
En la que es la fracción de energía disipada, , que se puede aplicar
Q0
en la forma:
( ) ( )( )
Q = QPLACA (1) + QPLACA (2) 1- QPLACA (1) + QPLACA(3) 1- QPLACA (1) . 1- QPLACA (2) (6.52)
221
221
221
Figura 6.6 Soluciones en forma de producto a los problemas de
222
222
222
Figura 6.7 Soluciones en forma de producto a los problemas de
223
223
223
6.9 Número de Biot (Bi)
semi-infinito
224
224
224
La ecuación diferencial para la distribución de temperatura T(x,t),
cuando las propiedades son constantes, es:
2T 1 T
(6.55)
x 2 α t
t
e (6.60)
x π x 2 αt παt
2Q -x 2 Q x x
T T 0
exp 0
= 1-erf (6.63)
i 0
kA 4 kA 2
226
226
226
6.13 PROBLEMAS RESUELTOS
SOLUCIÓN:
1. Diagrama de flujo
3) Para una placa de ancho (a), altura (H) y espesor (e) y despreciando el
efecto de los bordes, se tiene:
227
227
227
hL V a l H l
Bi L
K A 2 a H 2
h. l 400 w .0, 005m
Bi 2 m 2K 0.04 0.1
K w
50
mK
6) Calculo de la Temperatura:
Tt TF
e Bi.F o
Ti TF
Determine:
228
228
228
Figura 6.11 Cilindro, para evaluar la temperatura en el eje, así como a
una distancia radial en función del tiempo
Solución:
h
1)
500
22.2 condicion de contorno con convección
K 22.5
2.22 0.1
hR 500 x 0.1
Bi
K 22.5
Bi 1 0.45
229
229
229
Por tanto el tiempo que transcurrirá hasta que el eje del cilindro alcance una
temperatura de 250°C.
r 0.05
0.5 Bi 1 0.45
R 0.1
T R ,t 50 0.85250 50 220C
T R ,t 220C
También:
q(R,t ) q(0,t ) q(R,t ) q(0,t ) 250 - 50 q(R,t )
= x ; = = 0.21 ; = 0.85
R ,t T R ,t TF
0.21x0.85 0.1785
i Ti TF
Por tanto:
T R ,t 220C
230
230
230
6.8.3.- Una lámina de cierto material de propiedades térmicas
0.0002 m h
2
K 2 Kcal tiene un espesor de 3cm y se encuentra a
hmC,
231
231
231
una temperatura de 24°C. En estas condiciones, se introduce en un horno a
Determine:
Condición de criterios:
232
232
232
h 60
30 por convección
k 2
Se obtiene:
hR 60x0.015
Bi = = = 0.45 ; Bi - 1 = 2.22
K 2
i Ti TF 24 400
, t
Interceptando en la grafica con Bi 1 se obtiene :
i
2
t 2 x10 4 m xt
F0 3.56 h t 4horas
L2 0.015m2
x 1
Bi 1 2.22 0.66 se tiene que
L 1.5
T R ,t 400
0.92 T x 1,t 0.92(100) 400 345.6C
300 400
T x 1,t 345.6C T1
233
233
233
c) Cálculo del tiempo necesario, si el calentamiento continúa para que el plano
situado a 1 cm del plano central se alcancen 350°C.
234
234
234
æq ö- 1 T(x,t ) - TF 300 - 400
ç (x,t ) ÷
÷ = = mediante los parametros
ç
q(0,t ) ÷
ø T(0,t ) - TF T(R,t ) - 400 = 0.92
üï
x 1
= = 0.66 , Bi = 0.45ý en la Fig 4.13b
L 1.5 ïþ
T(0,t ) = 345.6°C
2
t 2 x10 4 m xt
Po lo tanto: F0 5,1 h t 5.7horas
2
L 0.015 m 2
2
235
235
235
x ,t T x ,t TF 200 400
se calcula con los parametros :
0,t T0,t TF T0,t 400
hL* 60 x0.03
Bi 0.9
k 2
Bi 1 1.11
T0,t 181.4C
se obtiene :
2
2 x10 4 m x t
t
Fo 3.83 *2 h t 17,335horas
L 3m 2
Determinar:
236
236
236
a) La temperatura máxima y mínima del cilindro treinta minutos después de
haberlo introducido en el horno.
237
237
237
b) El calor absorbido en ese tiempo.
SOLUCIÓN:
b) Qabsorbido t=30minutos
238
238
238
3) Cálculo del # Bi:
hL 30 x0.08
Bi 4.8
K 0.5
w
h 30 m 2 k
4) 60 1 (condición de criterios de convección)
k 0.5 w m
mk
p x , r ,t
P x C R
i
x ,t
P x para una placa inf inita, mediante Fig 2.37a y b
i
R ,t
C R para un cilindro l arg o(inf inito), Fig 2.38a y b
i
hL 30 x0.08
Bi 4.8 Bi 1 0.208
K 0.5
*Mediante la Fig. 4.13a, con Bi 1 y F0 se tiene:
239
239
239
0, L T0,t TF
P0 0.9
T T
i i F
x
Por tanto:
PL 0.27x0.9 0.249
Los parámetros:
2
5x10- 7 m x1800s
a t s
F0 = = 2
= 0.36
R2 (0.05m)
hR 30x0.05
Bi = = = 3 ® Bi - 1 = 0.333
K 0.5
R ,t
0.36 i
Fig.2.38b Bi 1 0.33
240
240
240
R/R0 1
Por tanto:
241
241
241
R ,t
CR 0.47 x0.36 0.1692
i
Tmax 480C
Q C V (Ti TF )
Se tiene que
Qp Qcil (1 Qp )
0.52 0.7(1 0.52) 0.856
242
242
242
Por tanto:
Q = F r C V (T
),i - TF
K k
a= ® rC = D = 0.1m H = 0.16m
r Cp a
k
Q V (T T )
i F
0.5 w mk
Q 0.856x 2
7
x R 2 H m3 20 500k
5x10 m
seg
0.5 w
Q 0.856x
m 2
7
mk x 0.052 x0.16H m3 20 500k
5x10
seg
1kcal
Q 516327Wseg 516327 Joule 516.327KJ
4.186J
Q 123.34Kcal
h m2 C
Determinar:
243
243
243
1. Diagrama de flujo
241
241
241
Figura 6.13 Cilindro, para evaluar la temperatura en el eje, así como el
calor absorbido
2) Cálculo del tiempo necesario para que la temperatura del centro sea de
138°C, ya que de los demás puntos será superior.
hL 20 x0.025
Bi 0.125 Bi 1 (0.125) 1 8
K 4
a t 0.04 m 2 h x t
F0 = = = 16 ® t = 0.25horas t = 15 min
L2 (0.025m)
2
i 0,t o
L ,t TL TF x
0.93 con 1
0,t T0 TF L
0,t Tc TF
P0 0.1517
0 T0 TF
PL 0.93x0.1517 0.1411
-Cilindro Infinito:
q(0,t ) Tc - TF
Bi - 1 = 1.333 ; C = = = 0.58 eje
(0) q(0) T0 - TF
q(EJE ) q(R ,t )
= 0.58 = 0.62
q(0) q(EJE )
Por lo tanto:
CR 0.58x0.62 0.3596
241
241
241
X TX TF
0.0818
0 T0 TF
TX 148.14C
TL 139.5C
2 36 üï
0.04 m x ï
a t h 60
F0 = = = 38.4°Cïï Q = 0.97
L2 (0.025m)
2 ýï
Q
0
ï
Bi 2F0 = 0.6 , Bi = 0.125 þï
Q = 0.97x51.24 = 49.70kcal
Problema N° 1
242
242
242
En un proceso de fabricación, de cilindros de acero inoxidable, que se
243
243
243
temperatura de 300 K con h = 500 W/m2.K. la longitud del cilindro es de 60 mm
enfriamiento y determine:
valores Tm = 450 K
Problema N° 2
Una esfera de acero con diámetro de 7.6 cm se tiene que endurecer calentándola
Calcule:
244
244
244
Una esfera de acero con diámetro de 7.6 cm se tiene que endurecer calentándola
Calcule:
Problema N° 4
inicial la temperatura del cilindro era igual en todas las partes de éste:
y la temperatura en la superficie.
245
245
245
TEMA VII TRANSFERENCIA DE CALOR POR CONVECCIÓN.
7.1 DEFINICIÓN.
Convección Natural
Convección Forzada
246
246
246
Pueden ocurrir situaciones donde ambas formas de convección actúen
simultáneamente, siendo la convección forzada de mayor interés practico que
la convección natural, debido a las aplicaciones industriales.
Donde:
Qc = calor transferido, w
A = área de transferencia, m2
5. Análisis numérico.
247
247
247
En este capítulo se incidirá en el estudio de convección forzada por el interior y
exterior tubos.
248
248
248
K= conductividad térmica del fluido (w/m. °C), evaluado en y = 0
(interface en la frontera sólido – fluido)
T
= gradiente de temperatura en el fluido en y = 0. La coordenada (y),
y y 0
se mide a lo largo de la normal a la superficie.
T - T f
Ts T f
-
h.Lc
Nu= (7.5)
k 0
249
249
249
Es la temperatura promedio del fluido en una sección transversal dada del
tubo, llamada también temperatura media de mezclado. Esta temperatura varia
de una sección transversal del tubo a otra, se expresa:
250
250
250
energía térmica total a travesando una sección
Tb
del tubo en una unidad de tiempo
(7.6)
capacidad calorífica del fluido a travesando la
misma sección en una unidad de tiempo
La energía térmica total que cruza una sección transversal del tubo en
una unidad de tiempo es:
R
C p T v 2 r d r
0
T
0
C pT v 2 rdr
(7.7)
b
R
C p v 2 rdr
0
251
251
251
Figura Nº 7.2 Tubería de sección circular, donde se muestra el volumen de
control para el balance de energía
251
251
251
Fuente: Elaboración propia, Ing Alberto Emilio Panana Girio
T 1 T r 2
r 2vo 1 2 r (7.10)
r r x R
250
250
250
5. Las constantes de integración se calculan teniendo en cuenta las
siguientes condiciones:
251
251
251
v T r 2 r 4
r 2
0
T T c v 2 v 0 1 (7.12)
α x 4
; 2
16R 2 R
T
k
h T T x r R
T (7.13)
k ; h
c p f c
x r R Tp T f
R R
Tf ρCp v2πrdr = TρCp v2πrdr
0 0
R r 2 2v T r 2 r 4
R
2v 1- T
0
rdr
T
Tvrdr
0
0 0
R 2
C
α x 4 16R 2
(7.14)
r2
f R
R
0
vrdr
2v0 1- 2 rdr
0
R
7 v0 R 2 T
T f TC
48 x
Tp T f r 3 v0 R T T 7 v0 R T 11 R
r=R 2 2
T
x x
C C
8
48 k
hc ; siendo D = diámetro= 2 R
11 D
253
253
253
k
hc 4, 3636 ; Nu = 4,3636 (7.17)
D
Esto puede ocurrir por ejemplo cuando fluye vapor condensado sobre la
superficie exterior de la tubería. A una distancia suficiente del punto en el
que empieza el calentamiento corriente abajo, el flujo se vuelve totalmente
desarrollado térmicamente, la forma del perfil de temperatura no cambia, y
el número de Russel tiene un valor constante dada por la ecuación:
a Tp Cte Nu 3.656 (7.18)
254
254
254
En la ecuación se tiene, (L) la longitud de la tubería y (d) el diámetro; el
parámetro (ηc), se utiliza para expresar el efecto de la diferencia de
255
255
255
temperaturas del fluido Tf y de la pared Tp, sobre las propiedades del fluido.
Se aplica en aquellos casos en que la viscosidad del fluido cambie
marcadamente con la temperatura, η =η (T); en muchos casos (ηc) se considera
la unidad, siendo de interés en los fluido muy viscosos.
Eecuación de Hausen
a. Con temperatura de pared constante, es:
0, 0668Gz
Nu 3, 66 (7.10)
1 0, 04Gz 2/3 C
d
0, 065 Re Pr
Nud 3, 66 L d ; Red < 2300 (7.22)
2
d 3
1 0, 04 Re d Pr
L
256
256
256
7.5.2 Ecuación de Sieder y Tate
η f
0,14
0, 8 1/3 Re >10000 ; L >60
Nu = 0,027 Re Pr con: d (7.24)
η p 0,7<Pr<16500
104<Re< 5.106; 0,5<Pr<200 ; error < 5 + 6%
cuyo campo de validez es: 104<Re< 5.106; 0,5<Pr<2000 ; error 10%
η
0< F <40
ηp
El valor del coeficiente de rozamiento viene dado para Pr > 0,5 por:
257
257
257
Tomándose las propiedades del fluido a la temperatura media TF, excepto
p que lo es a la temperatura de la pared Tp.
ν kf
258
258
258
d. Para una conducción formada por dos tubos concéntricos, el diámetro
hidráulico se determina mediante la relación siguiente:
Donde:
259
259
259
d, diámetro de los tubos.
Nu = C Re n Pr 1/3 (7.29)
Las propiedades del fluido se calcular a una temperatura media, entre la del
flujo TF y la de la pared exterior Tp.
260
260
260
proyecto de la inmensa mayoría de intercambiadores de calor., en la Fig. N° 7.5
el flujo forzado a través de un haz de tubos en batería.
a. Primer Método
Se utiliza una ecuación parecida a la de un solo tubo, en la que los
valores de C y n dependen de las distancias entre tubos adyacentes.
Estos parámetros varían si los tubos están alineados (disposición
regular), o están al tresbolillo o en quincunce, ambas disposiciones
triangulares, ver la fig. Nº 7.5
261
261
261
Para el caso de un flujo turbulento sobre baterías de 10 ó más tubos en
la dirección del flujo, se utiliza la siguiente ecuación y los valores C y n
se determinan de la tabla N°7.2
2000<Remáx <40000
Ren Pr1/ 3 ;
Nu C máx (7.31)
Pr>0,7
Tabla N°7.2 Para evaluar los valores de las constantes C y N
Sn / d
1.25 1.5 2 3
Sp / d C n C n C n C n
En línea
Al tresbolillo
1 - - 0.552 0.558 - - - -
Disposición triangular.
260
260
260
ST d v f ST
2 2 2
=> v máx = (7.34)
T S 2 d
S Paso mínimo
2 Ly
Donde: ST = espaciamiento transversal, desde los centro de tubos
2
S
S D T S L 2y (7.35)
2
b. Segundo método
Cuando el número (Nt) de tubo por fila sea superior a 20, se recomienda
utilizar la ecuación de Zukauskas, más moderna que la anterior, de la forma:
PrTF
Para gases, Nud = C* Re mmáx Pr 0,36 (7.36)
PrTp
m PrTF
Para líquidos, Nud = C* Re Pr 0,36 4 (7.37)
máx PrTp
0,7<Pr<500
6
10<Re med <10
C* y m están tabulados, tabla N° 4
261
261
261
Para líquidos, las propiedades se toman a TF, excepto los números de Pr
de la raíz, que lo son a las temperaturas respectivas.
Para haces con menos de 20 tubos por fila, N<20, el número de Nud
obtenido con la ecuación de Zukauskas se corrige mediante un factor de
corrección x que se determina a partir de la Fig. N° 7.6 o de la tabla N°
7.5, de tal forma que el factor de corrección es::
262
262
262
Figura N° 7.6 Gráfica para determinar el factor de corrección para el
número de filas en un banco de tubos, a utilizar en la
ecuación de Zhukauskas
263
263
263
1 S
P f N U 2
,
max
2
C2
f C1 Re m ,
Dm (7.39)
n
ST / D 1
Donde:
f = factor de fricción
N = número de filas de tubos paralelo a dirección del flujo
Problema Nº 1
264
264
264
a. Determinar el coeficiente de transferencia de calor superficial.
265
265
265
b. Encontrar la caída de presión del aire para este arreglo del banco de tubos.
Solución.-
1. Diagrama de flujo
A la temperatura Ts = 5 a.C.:
S S 2 T
3, 5
2 3, 808 cm
D L
2 2
4. El área frontal de flujo, de los tubos es: (ST – D) = 3,0 – 1,59 = 1,41 cm.
266
266
266
(ST - D) < 2 (SD – D); entonces la velocidad máxima ocurre en el
espaciamiento: (ST – D), y se determina, mediante la expresión:
267
267
267
ST 3, 0
V Vf 5 10, 638m / s
ST D 3, 0 1, 59
max
7.1 Para este número de Reynolds, mediante la tabla Nº 6,5 las constantes
C y n para la ecuación 6,54 a usar, con, ST/SL = 3,0/3,5 = 0,857, son:
C 0, 35 ST / S L 0, 35(0, 857)0, 339
1/ 5 1/ 5
n 0, 6
Pr f
1/ 4 1/ 4
0, 712
C Re Pr 0, 339 1,138 10 0, 714
n 0,36 4 0,36
N
u20
Pr 0, 717
s
Nu20 81, 44
k 0, 02545
h Nu 71, 66
D 0, 0159
h 114, 7w / m 2 .C (20, 2 Btu / h pie 2 F
11. Para evaluar la caída de presión a través del banco de tubos, se determina
primeramente el factor de fricción, para banco de tubos alternados:
1,08 1,138 10
0,472 4 -0,16
f = 1,0+ 0, 345
3,0/1,59-1
268
268
268
17,45
0,14
1
ΔP = 0,345 4 1,1644 10, 6382
2 18,65
ΔP = 90,1N/m2 (0, 013 l b f / pul 2 )
Problema Nº 2
Un haz de tubos utiliza una disposición en línea con Sn = SP = 1,9 cm. y 6,33
mm de diámetro de tubos. Se emplean 6 filas de tubos que constan de una pila
de 50 tubos de altura. La temperatura de la superficie de los tubos se mantiene
constante a 90 a.C. y transversalmente a ellos circula aire atmosférico a 20
a.C. siendo 4,5 m/s la velocidad antes que la corriente entre al haz de tubos.
Calcúlese el calor total por unidad de longitud transferido en el haz de tubos.
Estímese la caída de presión para esta disposición.
Solución.-
1. Diagrama de flujo
269
269
269
2. Las constantes para utilizar en la ecuación, pueden obtenerse en la tabla Nº
6,4, empleando
n
h.D V .D 1/ 3
Nu C f Pr f
k v f
Sp 1, 9 S
3, 00157 n
D 0, 633 D
C = 0,317; n = 0,608
Sn D 1, 9cm 0, 633cm
max
s
m
Vmax 6, 748
s
270
270
270
0, 317 2, 33104
h.D 0,608
Nu (0, 70751/ 3 127, 7619
k
127, 7619 0, 02836
h 572, 4056w / m2 .C
0, 00633
p p f2 f1
2
20 T f 2
Q 0, 94 572, 4056 0,119317 90 4, 6037 1007, 5 T 20
f2
2
Tf2 = 20,96 °C
271
271
271
Problema N° 3
271
271
271
En la calefacción de locales, normalmente se suele usar un haz de tubos en el
que el agua se hace circular en ellos y aire en flujo cruzado sobre ellos (por el
exterior). Considere un arreglo escalonado, el diámetro exterior de los tubos es
16,4 mm y los espaciamientos longitudinal y transversal: Sl = 34,3 mm y ST =
31,3 mm, existen 7 hileras de tubos en la dirección del flujo de aire y 8 tubos
por línea. En condiciones de operación típicos la temperatura superficial del
cilindro es de 70 a.C., mientras que la temperatura del flujo de aire cruzado y la
velocidad son 15 a.C. y 6 m/s respectivamente.
Calcular:
Solución.-
Ts = 70 °C : Pr = 0,701
2.Diagrama de flujo
270
270
270
Figura N° 7.9 Banco de tubos de arreglo alternado (triangular)
S S2 T
3, 43
2 37, 7 cm
D L
2 2
4. El área frontal de flujo, de los tubos es: (ST – D) = 31,3 – 16,4 = 14,9 mm.
ST D 31, 3 16, 4
max
7.1 Para este número de Reynolds, mediante la tabla Nº 6,5 las constantes
C y n para la ecuación 6,54 a usar, con, ST/SL < 2, por tanto
C 0, 35 ST / S L 0, 35(0, 91)0, 34 0n, 6
1/ 5 1/ 5
Pr f
1/ 4 1/ 4
0, 710
C Re n Pr 0,36 0, 34 13943 0, 71
0,6 0,36
N
u20
Pr 0, 701
s
Nu20 87, 9
271
271
271
9. El número de Nusselt calculado es para 20 filas de tubos, en la dirección
del flujo, como el banco de tubos es de ocho filas, se ha de hacer la
corrección respectiva, utilizando la figura Nº 6,7, obteniéndose:
272
272
272
Nu7
0, 98 Nu 7 0, 98 92, 5263 87, 9
Nu20
h 135, 6w / m 2 .K
Mediante la relación:
DNh
T T T T f ,i e
V f Nt ST Cp
s f ,s s
12 Reemplazando valores
DNh
V f Nt ST Cp
T T T
s f ,s s
T
f ,i
e
(,o 1o64).56.135,6 W/ m 2. K
T T 55 e
C
1217 x 6 X 8 X 0,0313
44, 5C
s f ,s
Q N h Tml 56135, 6x0, 0164x49, 6
Q 19, 4KW / m
Problema Nº 4
274
274
274
b. La caída de presión a través del banco de tubos.
Solución
Datos:
N1 144 tubos
D 1, 5 cm ; L 1 ; ST S L 2 cm
TS 50 k ; T f 300 k ; v 6m / s
1/ 2 1/ 2
1) Cálculo S D ST2 S L2 22 22 2, 828cm
ST D 2 1.5 0.5
2 S D D 2 2.828 1.5 2.65 ; Como ST D 2 SD D
2) La velocidad máxima
ST
v v
máx f
ST D
max 6 2
24 m / s
2 1.5
3) Cálculo Rmax
max D 24 m / s 0.015
R 13899, 6135
275
275
275
……. Flujo turbulento
max
25, 9x106 m2 / s
276
276
276
4) Para lograr que el Nº Re se determina las constantes C y por la ecuación
de Zukauskas para banco de tubos
1/ 4
n0,36 Pr f
Nu CR P
e r
Pr s
Como: ST / S L 1 0, 7
C 0, 27 ; n 0, 63
1/ 4
0, 708
Nu20 0, 27 13899, 6138 0, 689 0,36
0,63
0, 68
Nu20 97,18416 f c 0.98
h 213, 6561 W / m2 º C
Tf 2 Tf 1
Q mCp T f 2 T f 1 hA Ts
2
m 1, 695436 Kg / s
6).Cálculo de Área
7).Reemplazando en (5)
277
277
277
Tf 2 300
1, 695456x1014 J / Kg º C T f 2 300 213, 6561 x 6, 7858 500
2
T f 2 418, 65 K
máx
f 2
f
C2
f C1 R m
c
ST / D 1n
C1 0.176
2
C2 0, 34 SL / D 0, 34 x 0, 4533
1, 5
1,13x1.5
n 0, 43 1, 277.5 ; m 0.15
2
f 0.176
0.4533 1.3899, 61380.15 ; f 0, 483098
1.2775
2
1.5 1
278
278
278
Problema Nº 5
279
279
279
Fluye aire a razón de 5pies/s, sobre un intercambiador de calor de flujo
cruzado, que consta de siete tubos en la dirección del flujo y ocho tubos en la
dirección transversal al flujo. La longitud de cada tubo es de 4pies. El diámetro
exterior de los tubos es de ¾ de pulgada, la separación longitudinal es
SL=1.5pulg. Y la separación transversal ST=1.125pulg. La temperatura del aire
que entra en el intercambiador del calor es de 400°F y la temperatura de la
superficie de los tubos se puede considerar como 200°F. Las propiedades del
aire alas siguientes temperaturas son:
1. Diagrama de flujo
280
280
280
2. Cálculo del coeficiente h
S D 1.875; ()S2T() D S D D
0.751.125
1.125
5x 15 pies / seg
1.125 0.75
max
RE max max xD /
15x0.75 / 12
RE max 3060.946
1.1026x1 / 3600
k rs
C = 0.27; y n = 0.63
hD 0.701
0.27(3060.946)(0.063.701)()370.3.62 6 1/ 4
k 0.705
2.4 despejando (h)
37.26 x0.01995
h 11.89W / h. pieF
0.75 / 12
281
281
281
h 11.89 x0.955 11.3549
282
282
282
3. Para evaluar la transferencia de calor y la temperatura de salida del
fluido se procede:
Tf 1 Tf 2
Q hA(T
)() mCp T T
s f1 f2
2
0.75
A NT L 8x7 x x x4 43.9824 pie 2
12
m N S 4 0.04614 5 8
1.125
4 0.6921
f t T
12
4. Reemplazando
BTU 1h 400 T f 2
Q 11.3545 43.984(200 )
h. pie F 3600seg
2
2
lb BTU
Q 0.6921 0.2429 (T
400
f2)
s lbF
Tf 2 283.16F
()
t m a x
2 f
f C
C 2 R m
1
S C
(1 ) T D n
1.13D
n 0.43 m 0.15
ST
C2 0.34 ()0.68
1.15
0.75
1.13 0.75
n 0.45 ()1018 33
1.125
0.68
f 0.176 (3060.946) 0.15
0.4403
f 0.516
283
283
283
1 0 .0 5 1 9 9
P 0 .5 1 6 0 0 .0 4 6 1 4 1 5 2 ( 0 .1 4
)
2 0 .0 6 2 8 0
P 1 8 .2 5 9 6 lb f / p u lg 2 0 .1 2 6 lb f / p u lg 2
P.7.10.1
P.7.10.2
En el interior de un haz de tubos se utiliza vapor que condensa a 150 a.C. para
calentar una corriente de CO2, que entra a 3 atm,, 35 °C y 5 m/s. El haz de
tubos consta de 100 tubos de 1,25 cm de diámetro exterior en disposición
cuadrada en línea con espaciamientos Sn = Sp = 1,875 cm. Los tubos tienen
una longitud de 60 cm. Suponiendo que la temperatura de la pared de os tubos
se mantiene constante a 150 °C, Calcúlese el calor total transferido al CO2 y su
temperatura de salida del CO2 del banco de tubos.
P7.10.4
284
284
284
longitud del tubo necesaria para elevar la temperatura del volumen másico de
mercurio a 230°C, sin que se genere vapor de mercurio,
281
281
281
y determine el flujo de calor correspondiente. El punto de ebullición de mercurio
es de 355°C.
P7.10.5
P7.10.6
Fluye agua a una velocidad másica de 2.6Kg/seg. A través de una tubería liza
de 2 pulgadas de diámetro cedula 40 y es calentada de 20°C a 40°C. La pared
interna del tubo se mantiene a 90°C. Calcular:
a) La longitud de la tubería.
b) Cual es la caída de presión que experimenta el agua en la tubería.
Datos:
A la temperatura; Tm 30C
D0 6.033cm
para D 2 pu lg D
i 5.25cm
8.1 INTRODUCCION
280
280
280
Un Intercambiador de calor es un sistema mecánico, construido para transferir
calor entre dos fluidos a diferente temperatura que están separados por una
pared que puede ser metálica.
8.2 CLASIFICACION
281
281
281
Compactos: tubos – aletas, placas – aletas.
282
282
282
largo del intercambiador. En el flujo en contracorriente la temperatura
final del fluido más frió (que es lo que se calienta) puede superar la
temperatura de salida del fluido mas caliente (que se enfría), puesto que
existe un gradiente de temperatura favorable a todo lo largo del
intercambiador de calor.
283
283
283
Figura N° 8.2 Intercambiador de flujo cruzado
Fuente : Alan Chapman, Fundamentos de Transferencia de calor, 3era
Edición
284
284
284
Figura N° 8.3 Intercambiador de carcasa y tubos (1-2)
Fuente : Alan Chapman, Fundamentos de Transferencia de calor, 3era
Edición
285
285
285
Figura N° 8.4 Distribución de trmperaturas en :
a. Condensadores de un solo paso en los tubos
b. Vaporizadores de un solo paso de tubos
c. Intercambiador de flujos de calor en equicorrientes
Fuente: Alan Chapman, Fundamentos de Transferencia de calor, 3era
Edición
286
286
286
8.4 COEFICIENTE GLOBAL DE TRANSFERENCIA DE CALOR
i 1 hc A kA hF A
1
UA (8.4)
r
1n e
1 ri 1
hci Ai 2 kL hFe Ae
288
288
288
Mientras que si viene referido a la superficie interior Ai será:
1
U i (8.6)
re
A ln
i
1 ri Ai
hci 2 kL Ae h Fe
289
289
289
1 1 1
RSucio RFunc RLimpio U Fun (8.7)
U Func U Limpio 1
RSucio
U
Limpio
1 1 T T
... ( 1n ) T 2
A a a bT 1
290
290
290
U1T2 U 2 T1
DMLC
U T
1n 1 2 (8.10)
U 2 T1
DMLC (Diferencia media logatrimica cruzada)
U1T2 U 2 T1
QA (8.11)
U T
1n 1 2
U 2 T1
290
290
290
La diferencia de temperaturas media aritmética no difiere de la (LMTD) es
más de 1% y se puede utilizar en lugar de ella para simplificar los cálculos.
291
291
291
8.9 FACTOR DE CORRECIÓN DE LA DIFERENCIA MEDIA LOGARITMICA
DE TEMPERATURAS (LMTD)
mF C pF CF TC1 TC 2 1
R (8.17)
mC C pC CC TF 2 TF 1
292
292
292
8.10 EFECTIVIDAD Y NUMERO DE UNIDADES DE TRANSMISION
Para resolver este problema Nusselt desarrollo el método que lleva por
nombre Número de Unidades de Transmisión (NTU). Este método fue
perfeccionándose después por Kays y London. Consiste en determinar el
intercambio de calor por cada grado de diferencia de temperatura, que
según la expresión (18) es el producto UA.
Coeficiente de capacidad,
()mC p min
CR (8.20)
()mC p max
293
293
293
Efectividad,
T T
fs fe , m C mC (8.21)
Tce T fe
f pf c pc
Tcc Tcs
,m C m C (8.22)
Tcc T fe
f pf c pc
Problema Nº 1
Aceite caliente (Cp = 2 200 j/Kg. °C) se va a enfriar por medio de agua (Cp = 4
180 j/kg. °C) en un intercambiador de calor de dos pasos por el casco y 12
pasos por los tubos. Estos son de pared delgada y están hechos de cobre con
un diámetro de 1,8 cm. La longitud de cada paso de los tubos en el
intercambiador es de 3 m y el coeficiente de transferencia de calor total es de
340 W/m2°C. Por los tubos fluye agua a una razón total de 0,1 Kg./s y por el
casco fluye el aceite a razón de 0, 2 Kg./s. El agua y el aceite entran a las
temperaturas de 18 .C. y 160 .C. respectivamente. Determine la velocidad de
transferencia de calor en el intercambiador y las temperaturas de salida de las
corrientes del agua y del aceite.
Diagrama de flujo
294
294
294
FIGURA Nº 8.7 Intercambiador de calor de carcasa y tubos de fluidoa agua y aceite
Solución.-
Cálculo de NTU y Rc
W W
340 2 2, 0357m 2 418
NTU m .C 1, 655 Rc C 0, 95
W W
418 440
C C
295
295
295
T2 F -18
0, 61 T 104, 62C
160 18
F2
El flujo de calor:
W
Q=C T -T 418 104, 62 18 F 36207,16W 36, 20716KW
F
F F2 F1
Problema Nº 2
Solución.-
296
296
296
297
297
297
FIGURA Nº 8.7 Intercambiador de calor de tubos concéntricos, para fluidos agua
y alcohol etílico en corrientes paralelos
2. Datos: U = 568 w/m. °C; Do = 25,4 mm.; Cp.= 3810 j/Kg.K (alcohol
etílico);
Ln
T2 T2 TC2 TF2
Ln
39, 4 36, 23
T1 65, 6 10 1
TML 18, 3 C
Q 691766(/j )s
Q U .A.TML A 66, 55 m 2
A 66, 55m 2
A .D.L L 834 m
D 0, 0254m
298
298
298
FIGURA Nº 8.8 Intercambiador de calor de tubos concéntricos, para fluidos agua
y alcohol etílico en contracorriente
T
Ln 2
T2 TC2 TF1
Ln
39, 4 10
T1 65, 6 36, 23
0
TML , las diferencias de temperaturas en los terminales son
0
3. Puede ocurrir que las capacidades caloríficas de los fluidos son iguales,
las diferencias de temperaturas en contracorriente resultan iguales, ΔT=
ΔT1 = ΔT2, por que se ha de aplicar la regla de L´Hôtipal.
Q UATML
T T 0 T (1x)
TML 2 T xT 1 L´Hôpital xT T
T2 0
2 1
Ln x
1 2
Ln
T1
TML (T
)6C5
1
, 6TF 236, 23 29, 37 C
A 41, 47m2
A .D.L L 519, 69 m
D 0, 0254m
299
299
299
c. Intercambiador de calor de carcasa y tubos, dos pasos por la carcasa y
cuatro pasos por los tubos (2-4)
300
300
300
FIGURA Nº 8.9 Intercambiador de calor de carcasa y tubos, para fluidos agua y
alcohol etílico
5. De la ecuación de diseño
Q 691766w
Q U .A.F . T A
w m2 C C
T ML
U .F . T
T ML
A 42, 75m 2
A L 1.86m
L
4 72 D
301
301
301
4
2
,
7
5
m
2
4
7
2
0
,
0
2
5
4
m
302
302
302
d. Intercambiador de flujo cruzado con un paso de tubos y un paso por la
carcasa, siendo con mezcla de fluido por la carcasa
.
FIGURA Nº 8.9 Intercambiador de calor de de flujo cruzado, para fluidos agua y
alcohol etílico
3. De la ecuación de diseño
Q 691766w
Q U .A.F . T A
w m2 C C
T ML
U .F . T
T ML
A 47, 39m 2
Problema Nº 3
303
303
303
En una planta textil se va a usar el agua de desecho del teñido (Cp = 4 290
j/kg. °C) que está a 75 .C. para precalentar agua fresca (Cp = 4 180 j/kg. °C) a
301
301
301
15 °C, con el mismo gasto de masa., en un intercambiador de calor de tubo
doble y a contraflujo. El área superficial de transferencia de calor del
intercambiador es de 1,65 m2 y coeficiente de transferencia de calor total es de
625 W/m2. °C. Si la velocidad de la transferencia de calor en el intercambiador
es de 35 KW, determine la temperatura de salida y el gasto de masa de cada
corriente de fluido
Solución.-
m2
1. Datos: Área = 1,65 ; U = 625 W/m2. °C ; Q = 35 000 W
T C
TF TC TF
Q U A T U A 1 2 2 1
ML
T C TF
1 2
Ln
T C TF
2 1
Cp
T T F TF2 TF1
C2 C1
CpC
301
301
301
Luego mediante la ecuación de diseño se halla el calor transferido, si
no es igual al flujo de calor dado se procede a realizar otras
suposiciones, hasta encontrar un valor cercano.
TC TF 75 15
Asumiendo que TF2 = 45 °C; T 1 1
45C
F
2
2 2
j
4180
kg.C
TC2 75C
j
45 15 C 45, 80C
4295
kg.C
75 45 45, 80 15
TML 30, 4002C
75 45
Ln
45, 80 15
W
Q U A T 625 1, 65m 2 30, 4002C 31350, 235W
ML
m 2 .C
Como no se alcanza el valor del calor transferido, se produce a realizar
otras suposiciones, los cálculos realizados se muestran en la siguiente
tabla:
F F F2 F1 F
Cp (T
)41
80T(41, 4 15) s
302
302
302
F F2 F1
Problema Nº 4
303
303
303
En un intercambiador de calor con flujo en contracorriente, por el que circulan
5 Kg. de agua por minuto, y 8 Kg. de aceite por minuto, el agua entra a 20 °C y
sale a 40 °C, mientras que el aceite entra a 90 °C. El calor específico del agua
es Cpagua = 1 Kcal./kg. °C, el calor específico del aceite obedece a la siguiente
relación, Cpaceite = 0,8+0,002Taceite, con Taceite = °C.
Determinar:
Kcal 10Taceite
Dado por: U , con T en C
min .m .C
2
Taceite Tagua
Solución.-
5. Integrando la expresión:
T
0, 002T 2 2
C
mac 0, 8T
mag .Cpag TF TF mag .Cpag TF TF
2 TC1
2 1 1 2
0, 002T 2
0, 002 902
8 0, 8T 0, 8 90 5 1 40 20
C2
C2
2 2
TC2 77, 07 C
Qreal ,agua mag Cpag Tag 5 1 (40 20)100 Kca/ lmin
306
306
306
La velocidad máxima posible de transferencia de calor para el aceite
( Qmax ,aceite ) es cuando: Tc2 = Tf1
0, 8 0, 002Tac dTac
C1 C1
Qmax,aceite macCp acdTac m ac
F F
1 1
kg 902 202
8 0, 8 90 0, 002 0, 8 20
ac
Qmax,
min
aceite
2 2
Kcal
Q
max,aceite 509, 6
min
100
0, 2857 28, 7%
350
Cálculo de la ΔTML
T2 T1 T2 TC1 TF2 90 40 50C 50 57, 07
T
De la ecuación de diseño
Kcal Kcal
Q U .A.TML 100 U .A 53, 45C U .A 1, 87
min min .C
U .A 1, 87
NTU 0, 374
C min 5
308
308
308
La ecuación para determinar la eficiencia , para un intercambiador de
calor a contracorriente, es:
C 5
NTU min
1 0,374 1
1 e Cmax
1 e 7,736
C
5
0, 2857 28, 57%
NTU min
1 5 0,374 1
7, 736
)(F) TF
CF (T TF TF 40 20
2 1
2 1
0, 2857 28, 57 %
C T )90T 20 T T
min ()(C1 F1 C1 F1
Tac Tag
ac ac ac ac ag ac ag ac
TC
(0, 8 0, 002Tac )dTac mac
TC1
Am 0, 8Ln 1
0, 002 T T
ac
TC2 10Tac
10 TC2
C1 C2
90
A 0, 8 0, 8Ln 0, 002 90 77, 07
77, 07
A 0,11988m 2
309
309
309
10. También se puede determinar, de la siguiente forma:
U A 1, 87
Kcal
A
1, 87 Kcal min .C
0,11987m 2
min .C 15, 6
Kcal
m2 . min .C
310
310
310
Problema Nº5
Datos:
Para el NH3:
ρ = 580 kg/m3; Cp = 5 KJ/Kg.°C; k = 0,50 w/mK; ν = 0,34x10-6 m2/s;
Pr = 2
Para el agua:
ρ = 985 kg/m3; Cp = 4,186 Kj/kg.°C; k = 0,66 w/mK; ν = 0,484x10-6 m2/s;
Pr=3
Solución.-
1. Diagrama de flujo
311
311
311
Fuente: Elaboración propia, Ing. Alberto Emilio Panana Girio
kamon
D32 D22
D 4 4 D D 100 60 40mm
D D
H 3 2
3 2
313
313
313
3.3 Cálculo del número de Nusselt, para luego hallar el coeficiente por
convección por el lado del agua (ho), para esto se hará uso de la
ecuación de Dittus-Boelter, este fluido sufre enfriamiento
kagua
1 1
Uo
r2 r2 r2 1 30 0, 03 30 1
Ln Ln
w
Uo 2400
m 2 .K
Para el amoniaco:
m.Cp
kg Kj Kj Kj
C 3 12300 5 61500 17, 08
NH NH 3
h kg.C h.C s.C
Para el agua:
314
314
314
Cálculo del flujo másico del agua
D32 D22 0,12 0, 062 m2 m kg
mH 2O Qv v.A. v. 1, 5 985
4 4 s m3
kg kg
mH 2O 7, 4267 26736
s h
NTU 2, 6486
Cmin Kj
17, 08
s.C
Cmin
17,03
0, 8361
NTU 1 2,6486 1
17,03
Cmin
1 e Cmax
1 e 31,088
Cmax 31,088
315
315
315
5.5 Las temperaturas de salida de los fluidos serán :
C
TC2 TC1 TC1 TF1 min 80 (80 20)0, 8361 0, 5494
52, 5 C
C
max
316
316
316
TF2 TF1 TC1 TF1 Cmin
C
20 (80 20)0, 8361 1
70,17 C
F
T2 T1
Q UA
C min T T
T2 C1 F1
Ln
T1
T2 TC TF 80 71,17 9, 83
1 2
T1 TC TF 52, 5 20 32, 5
2 1
Kj
Q 0, 836117, 08 80 20 856, 8Kw
s.C
Problema N° 6
Solución:
1. Diagrama de flujo
310
310
310
311
311
311
FIGURA Nº 8.13 Intercambiador de calor de flujos cruzados
A = Q / U T FT
(80-30) (225-30)
P= = 0.26 ; R= = 2.5
(225-30) (80-30)
FT = 0.92
A = Q / U T FT = 627.600/(200*0.92*103)= 33.1 m2
PROLEMA N°1
312
312
312
a. El calor transferido.
b. La eficiencia del intercambiador.
PROBLEMA N°2
PROBLEMA N°3
PROBLEMA N°4
b. Con todas las demás condiciones iguales, pero teniendo en cuenta los
cambios en el coeficiente global. Grafique la temperatura de salida del
313
313
313
agua de enfriamiento (Tc) y el flujo de condensación de vapor ( h) como
Problema N°5
340 W . Por los tubos fluye el agua a una razón total de 0.1Kg/s y por la
m 2 C
coraza fluye aceite a razón de 0.2Kg/s. El agua y el aceite entran a las
temperaturas de 18°C y 160°C respectivamente.
Determine:
314
314
314
IX. REFERENCIALES
edición, 2001
edición., 2002.
1986
315
315
315
10. OCÓN GARCÍA JOAQUÍN, TOJO BARREIRO GABRIEL.
Edición, 2005
INTERNET
15. http://es.wikipedia.org/wiki/intercambiador_de_calor.
16/09/2007.
17. http://www.scielo.org.ve/cielo.php?pid=S0378-
18442001000900003&script=sci_artt... Efecto en la
316
316
316
hidrodinámica y transferencia de calor del desfasamiento entre
30/03/2005.
18. http://www.scielo.cl/scielo.php?pid=S=718-
07642004000400007&script=sci_arttext. Evaluación de
16/09/2007.
19. http://www.itsimo.edu.mx/metalmecanica/tema1intercambiadores
decalorudemavarra.pdf
20. htpp://torresdeenfriamiwento.blogspot.com/2007/02/torres-de-
enfriamiento.html.
317
317
317
APENDICE
317
317
317
NOMENCLATURA
Sistema
318
318
318
l longitud de la aleta ( distancia media entre las placas) m
N número de placas
P Perímetro de la aleta m
Q cantidad de calor W
exterior. m
t espesor de aletas mm
T Temperatura °K, °C
Letras griegas
espesor m
Eficiencia de la aleta
diferencia de valores
µ Viscosidad dinámica N.s/m2
319
319
319
ρ densidad kg/m3
φ Factor de escala
Números adimensionales
320
320
320
ANEXOS
321
321
321
FIGURA N° 8.14 Diagrama del factor de corrección F para intercambiadores
de calor comunes de casco y tubos y de flujo cruzado
(tomado de Bowman, Mueller y Tagle, Ref.2)
322
322
322
Fuente: Cengel Yunus A, Transferencia de calor, Tercera Edición, 2004
323
323
323
FIGURA N° 8.15 EFECTIVIDAD PARA LOS INTERCAMBIADORES DE
CALOR (TOMADO DE KAYS Y LONDON, Ref. 5)
324
324
324
325
Fuente: Cengel Yunus A, Transferencia de calor, tercera Edición, 2004
326
326
326
Figura N° 8.16 Para un NTU y una relación de capacidades C dados, el
intercambiador de calor a contraflujo tiene la efectividad más
alta y el de flujo paralelo, la más baja
327
327
327
Tabla Nº 6 Factores de resistencia por ensuciamiento normales
Asfalto 0,0050
Gasolina 0,0010
Queroseno 0,0010
328
329
UNIVERSIDAD NACIONAL DEL CALLAO
Señor:
Presente.-
De mi mayor consideración:
Atentamente,
330
330
330
Ing. Alberto Emilio Panana Girio
Investigador responsable
Señor:
331
331
331
Presente.-
De mi mayor consideración:
Atentamente,
Investigador responsable
332
332
332
TQ a===04=
qLamax
Q 2,950°07Ch
.3528.= .8*
Kcal/h
Kcal/hm
A *(T -T ) 2
QmáxO F
333
333
333