Universidad de Chile

Facultad de Ciencias Químicas y Farmacéuticas

Departamento de Ciencia de los Alimentos y Tecnología Química
Operaciones Unitarias I (QF)

PRÁCTICO Nº3
Transferencia de calor
Intercambiadores de calor multitubulares

Integrantes: Juan Osorio, Valentina Osses y Pamela Ponce

Fecha del trabajo práctico: lunes 3 de junio, 2019

Fecha de entrega: miércoles 12 de junio, 2019

Grupo: G6

Profesor Coordinador: Carlos Basualto F.

Ayudantes: Claudia Carrasco F., Felipe Soto Silva,

José Gaete y Matías Palma M.
Introducción
Calor es la transferencia de energía a través de la frontera de un sistema debido a una diferencia
de temperatura entre el sistema y su entorno. 1

Cuando dos objetos que están a diferentes temperaturas se ponen en contacto térmico, el calor
fluye desde el objeto de mayor temperatura hacia el de menor temperatura. 2 Esta transferencia de
calor puede ser mediante 3 diferentes mecanismos dependiendo del medio a través del cual fluya,
estos son:

- Conducción: Cuando se transmite por sólidos

- Convección: Cuando se transmite por fluidos

- Radiación: Cuando se transmite a través del vacío.

Las aplicaciones a nivel industrial son variadas, pero se pueden resumir básicamente en dos:
calentamiento y enfriamiento. Por ejemplo, en el caso de mezclado o secado de productos se debe
aplicar calor, es decir, calentar; mientras que, si se necesita preservar un producto o aplicar
diferentes métodos de extracción, como por ejemplo condensación, es necesario enfriar.

q=U x Á x ∆ T ln Ecuación N°1 Para medir la transferencia de calor de los intercambiadores de

calor tubulares

Donde:

q= flujo de calor cedido por el fluido caliente

U= Coeficiente integral de transferencia térmica, que puede ser interno o externo, dependiendo
del área que se considere.

A= área intercambiadora de calor, que corresponde a la del manto del cilindro (tubo).

∆ Tln = Diferencia de temperatura media logarítmica

∆ T 2−∆ T 1
∆ Tln= Ecuación N°2 Diferencia de temperatura media logarítmica
ln ( ∆ T 2 /∆ T 1 )
Donde:

∆ T 1= Diferencia de temperaturas terminales fluido caliente

∆ T 2= Diferencia de temperaturas terminales fluido frío

Para este trabajo experimental se utilizó un intercambiador de calor multitubular para docencia,
con el fin de comparar la eficiencia de los diferentes tipos de flujos en el intercambio de calor.
 T 1−T 2
Eficiencia= x 100 % Ecuación N°3 Eficiencia para enfriar en contracorriente
T 1−T 3

Donde:

T 1: Temperatura del fluido caliente entrante

T 2: Temperatura del fluido caliente saliente

T 3: Temperatura del fluido frío entrante

T 4: Temperatura del fluido frio saliente

T 2−T 1
Eficiencia= x 100 % Ecuación N°4: Eficiencia para enfriar en paralelo
T 2−T 3
Donde:

T1: temperatura del fluido caliente saliente

T2: Temperatura del fluido caliente entrante

T3: temperatura del fluido frio entrante

T4: Temperatura del fluido frio saliente

m
q= xCpx ∆T Ecuación N°5: Calor sensible
θ
Donde:

m
= Velocidad volumétrica
θ
C p= capacidad calorífica del agua

T= diferencia de temperatura del fluido caliente (entrada y salida)

D xV x ρ
R= Ecuación N°6: Número de Reynolds
µ

Donde:

D = Diámetro interno

V = Velocidad lineal

ρ = Densidad del fluido

µ = Viscosidad dinámica

Q=VxA Ecuación N°7: Caudal

Donde: A = Área

2
π xD
A= Ecuación N°8: Área de pasada del fluido caliente
4

π x De 2
A= Ecuación N°9: Área de pasada fluido frio
4

Donde:

De = Diámetro equivalente

Ai−7 x Ae
De= Ecuación N°10: Diámetro equivalente
π x Di+7 x De x π

Donde:

Ai = Área interna de la carcasa

Ae = Área externa de los tubos internos

D i= Diámetro interno de la carcasa

De = Diámetro externo de los tubos internos

Objetivos

I. Estudiar el funcionamiento de un intercambiador de calor multitubular completamente

automatizado y controlado remotamente mediante un software relacionando las medidas
experimentales con las ecuaciones de transporte.
 II. Comparar el intercambio de calor para un proceso en función de la configuración de
trabajo (paralelo y en contracorriente).
III. Comparar el intercambio de calor en función de los caudales de los fluidos participantes.
IV. Comparar resultados de los datos calculados por Ud. empleando las ecuaciones de
transporte, con los resultados que entrega el software del equipo.
V. Determinar el coeficiente integral de transferencia térmica (U) para el proceso empleando
ambas configuraciones: paralelo y en contra corriente.
VI. Determinar los números de Reynolds de los fluidos participantes y analizar su influencia en
las ecuaciones de transporte estudiadas en clase.

Materiales, equipos y métodos

I. Materiales
i. 4 termocuplas

II. Equipo
i. Intercambiador de calor para docencia marca Armfield modelo HT-30XC, del tipo
multitubular modelo HT-33.

Especificaciones técnicas del equipo:

Número de tubos= 7 en disposición triangular

Diámetro interno de tubo = 5,15 mm

Espesor de pared de tubo = 1,2 mm

Largo de tubo= 144 mm

Diámetro interno carcasa= 40,61 mm

Espesor de la carcasa= 5 mm

ii. Método

i. Se conectó el equipo al software mediante el cable USB que tiene el equipo (asegurar de
que las 4 termocuplas estén correctamente instaladas en los conectores específicos que
poseen).
ii. Se verificó que la bomba externa esté acoplada y correctamente instalada al
intercambiador, y que las válvulas estén en su posición correcta.
iii. En el computador de trabajo, se abrió el software HT33 Shell and Tube Heat Exchanger
versión 2.1 instalado, el cual corresponde a un intercambiador de calor multitubular.
 iv. Para el primer experimento se eligió el modo en contracorriente, luego de ello se
seleccionó el ícono que está en el panel de control y se conectó el software al equipo
(verificar si el equipo está efectivamente en contacto con el software, mediante el botón
que tiene un teléfono rojo de referencia)
v. Se hizo llegar lentamente al equipo el agua fría y se abrió lentamente la válvula hasta que
se llegó al caudal deseado en la opción cold flow
vi. Se reguló lentamente el flujo del agua caliente en la opción hot Flow y su temperatura de
trabajo en el calentador a través de la opción heater en el panel de control. El agua
caliente se recircula continuamente en el equipo, mientras que el agua fría entra y sale sin
recirculación.
vii. Se inició la captura de datos en el botón GO, regulado en este caso para obtener datos
cada 10 segundos para un tiempo total de 20 min.
viii. Se analizaron y registraron las 4 temperaturas que dan las termocuplas en línea. Verifique
en función del tiempo todas las columnas de la planilla interna del software.
ix. Para cambiar al modo en paralelo se debe salir del software, abrirlo nuevamente y elegir
la opción en paralelo. Observe la nueva variación de los datos que entrega el equipo.

Se realizarán experimentos para estudiar el efecto de dos variables:

1) Efecto de la configuración en paralelo y en contracorriente.

En este caso empleé para ambos fluidos un caudal aproximado de 1 L/min, sin modificarlos
(relación de caudales 1/1). Permita un tiempo de equilibrio de 20 minutos antes de cambiar la
configuración. Al final del tiempo y previo al cambio de configuración, registre las 4 temperaturas
finales y con ellas haga los cálculos solicitados. Repita lo mismo para la nueva configuración,
también permitiendo un equilibrio de 20 minutos.

2) Efecto de la relación de caudales de transporte de los fluidos.

En este caso compararemos el efecto de la relación de caudales de los fluidos frío y caliente.
Manteniendo la configuración del experimento anterior, disminuya el caudal del fluido frío hasta
un valor aproximado de 0,5 L/min, manteniendo la del caliente en 1 L/min, lo que generará una
relación de caudales caliente/frío de 2/1. Registre las temperaturas para sus cálculos. Después de
20 minutos, cambie la configuración (paralelo/contracorriente) y vuelva a registra las
temperaturas. Compare todos sus resultados.
Esquema
Resultados
i. Construya los diagramas temperatura v/s largo del intercambiador para ambos
sistemas de operación del intercambiador.

Sistema en paralelo, relación de flujo 1:1

Temperatura Promedio de las últimas 10 temperaturas

[°C]

T1 36.968

T2 40.007

T3 19.090

T4 21.022

T2

T1

T4

T3

Largo del tubo

Sistema en paralelo, relación de flujo 2:1

Temperatura Promedio de las últimas 10 temperaturas

[°C]
 T1 15.820

T2 16.067

T3 12.602

T4 13.060

T2

T1

T4

T3

Largo del tubo

Sistema en contracorriente, relación de flujo 1:1

Temperatura Promedio de las últimas 10 temperaturas

[°C]

T1 19.366

T2 17.979

T3 12.105

T4 12.605

T1

T2

T4
 T3

Largo del tubo

Sistema en contracorriente, relación de flujo 2:1

Temperatura Promedio de las últimas 10 temperaturas

[°C]

T1 17.304

T2 16.628

T3 12.940

T4 13.550

T1

T2

T4

T3

ii. Calcule la eficiencia del intercambiador para el caso de flujo en contracorriente y

en paralelo, basándose en la variación de temperatura del fluido caliente.

Para el cálculo de la eficiencia del sistema en paralelo se utilizará la Ecuación N°4,

mientras que para el cálculo de la eficiencia del sistema en contracorriente se utilizará la
Ecuación N°3.
 T 2−T 1
Eficiencia= x 100 % Ecuación N°4 para sistemas en paralelos
T 2−T 3

Eficiencia sistema en paralelo, relación de flujo 1:1

40.007−36.968
Eficiencia= x 100 % = 14.528%
40.007−19.090

Eficiencia sistema en paralelo, relación de flujo 2:1

16.067−15.820
Eficiencia= x 100 %=7.128 %
16.067−12.602

T 1−T 2
Eficiencia= x 100 % Ecuación N°3 para sistemas en contracorriente
T 1−T 3

Eficiencia sistema en contracorriente, relación 1:1

19.366−17.979
Eficiencia= x 100 %=19.102 %
19.366−12.105

Eficiencia sistema en contracorriente, relación 2:1

17.304−16.628
Eficiencia= x 100 %=15.273 %
17.366−12.940

iii. Determine los coeficientes integrales de transferencia térmica (U) (interior y

exterior) para flujo en contracorriente y en paralelo. Realice estos cálculos
basándose en el calor cedido por el fluido caliente. Analice estos resultados.

m
q= x C p x ∆ T =U x Á x ∆ T ln
θ
Á INTERNA=7 x π x ( Dinternotubo ) x L

Á INTERNA = 7 x π X 5,15 mm x 144 mm= 16300,37 mm2

1 m2
16300,37 mm2 x ( 6
2
2 ) = 0,016 m 3
10 mm

Á EXTERNA =7 x π x ( D internotubo +2 xespersor ) x L

Á EXTERNA =7 x π x ( 5,15 mm+2 x 1,2 mm ) x 144 mm = 23896,65mm2

Sistema en paralelo, relación de flujo 1:1

2
1m
23896,65mm2 x ( 6 2 ) = 0,023 m
2

10 mm

Datos
3 3
L m m
 Caudal (Q) = 1 =0,001 =0,06
min min h
kg
 ρ = 992,8
m3
m m3 kg kg
 = Q x ρ =0,06 x 992,8 3 =59,568
θ h m h
Kj
 C p=4,178
kgxK

ΔT =313,157 K – 310,118 K =3,039 K

kg kJ kJ
q=59,568 x 4,178 x 3,039 K=756,331
h kg x K h

∆ T 1−∆T 2 (313,157 K −292,24 K )−( 310,118 K −294,172 K )

∆ T ln= = =18.319 K
ln
∆T 1
∆T 2( ) ln
( ( 313,157 K−292,24 K )
(310,118 K −294,172 K ) )
q
U EXT =
A EXT x ΔT ln

kJ
756.331
h kJ
U EXT = 2
=1795.074 2
0.023 m x 18.319 K m x hx K

U q
∫ ¿= A ¿
∫ ¿ x Δ Tln ¿

U kJ
756.331
h kJ
∫ ¿= 2
0.016 m x 18.319 K
=2580.419 2
m xhxK
¿

Sistema en contracorriente, relación 1:1

 Datos
L m3 m3
 Caudal (Q) =1 =0,001 =0,06
min min h
Kg
 ρ = 998,5
m3
m m3 Kg Kg
 = Q x ρ = 0,06 x 998,5 3 =59,91
θ h m h
Kj
 C p=4,182
Kg x K

ΔT = 292,516 K – 291,129 K = 1,387 K

kg kJ kJ
q=59.91 x 4.182 x 1.387 K =347,504
h kg x K h

∆ T 1−∆T 2 (292,516 K −285,755 K )−( 291,129 K−285,255 K )

∆ T ln= = =6.307 K
ln
( )
∆T 1
∆T 2
ln(( 292,516 K −285,755 K )
( 291,129 K−285,255 K ) )
kJ kJ
347.504 347.504
h kJ h kJ
Uext = 2
=2395.572 2 U Interna = 2
=3443.634 2
0.023 m❑× 6.307 K ❑ m ×h×K 0.016 m❑ ×6.307 K ❑ m ×h × K

Sistema en paralelo, relación de flujo 2:1

Datos
L m3 m3
 Caudal (Q) = 1 =0,001 =0,06
min min h
Kg
 ρ = 999,0 3
m
3
m m Kg Kg
 =Q x ρ=0,06 x 999,0 3 =59,94
θ h m h
Kj
 C p=4,185
Kg x K
 ΔT = 289,217 K – 288,97 K = 0,247 K

kg kJ kJ
q=59.94 ∗4.185 ∗0.247 [ K ] =61.960
h kg × K h
 ∆ T 1−∆T 2 ( 289.217−285.752 )−( 288.97−286.21 )
∆ T ln= → =3.106 ° K
ln
( )
∆T 1
∆T 2
ln
(
( 289.217−285.752 )
( 288.97−286.21 ) )
kJ
61.960
h kJ
Uext = 2
=867.326 2
0.023 m❑× 3.106 K ❑ m ×h × K
kJ
61.960
h kJ
U Interna= 2
=1246.780 2
0.016 m❑ × 3.106 K❑ m ×h×K

Sistema en contracorriente, relación 2:1

Datos
L m3 m3
 Caudal (Q) = 1 =0,001 =0,06
min min h
Kg
 ρ = 998.8 3
m
m m3 Kg Kg
 =Q x ρ=0,06 x 998,8 3 = 59,928
θ h m h
Kj
 C p=4,184
Kg x K
 ΔT = 290,454 K – 289,778 K = 0,676 K

kg kJ kJ
q=59.928 x 4.184 x 0.676 K =169,499
h kg x K h

∆ T 1−∆T 2 (290,454 K −286,7 K )− ( 289,778 K −286,09 K )

∆ T ln= = =3,721 K
ln
( )
∆T 1
∆T 2
ln
(
( 290,454 K −286,7 K )
( 289,778 K −286,09 K ) )
kJ
169,499
h kJ
U EXT = 2
=1980,522 2
0,023 m x 3,721 K m xhxK
U kJ
169,499
h kJ
∫ ¿= 0,016 m2 ×3.721 K =2847,000 m2 x h x K ¿
Análisis de resultados

iv. Compare los resultados anteriores para fluidos en contracorriente y en paralelo

con los que entrega el software.

Fluido Uext Uint Uexperimental

Paralelo 1:1 1795.074 2580.419 586.290

Contracorriente 1:1 2395.572 3443.634 816.854

Paralelo 2:1 867.326 1246.780 289.509

Contracorriente 2:1 1980.522 2847.000 646.278

v. Determine los números de Reynolds para ambos fluidos (considere para fines
prácticos densidad abs= 1g/cm³; µ= 1 mPa·s).

 Para este cálculo primero debe conocerse la velocidad de cada fluido, que se calcula a
partir del caudal y del área circular del tubo que lo contiene.

•Número de Reynolds para el fluido caliente con caudal de 1 L/min

3
−5 m
1.6667 x 10 x4
Qx4 kg s
ρx x Di 1000 3 x x 5,15⋅10−3 m
2 −3 2
D xV x ρ πx Di m πx ( 5,15⋅10 )
R= = =
µ µ kg
1 x 10−3
m∙ s

R=4120,598
 Los resultados indican que el fluido estudiado se mueve con un flujo de carácter
turbulento.

•Número de Reynolds para el fluido frío con caudal 1 L/min

−3 2 −3 2
( 40.61 x 10 m) −7 x (7.55 x 10 m)
De=4 x =0.0053 m
( 40.61 x 10−3 m ) +7 x (7.55 x 10−3 m)

−5 m3
1.6667 x 10 x4
Qx 4 kg s
ρx 2x
x De 1000 3 x 2
x 0.0053 m
D xV x ρ πx De m πx ( 0.0053 m )
R= = =
µ µ kg
1 x 10−3
m∙s
R=4003,978

Los resultados indican que el fluido estudiado se mueve con un flujo de carácter
turbulento.

•Número de Reynolds para el fluido frio con caudal 0,5 L/min

3
m −6
8.333 x 10 x4
Qx4 kg s
ρx 2x
x D e 1000 x x 0.0053 m
D xV x ρ πx De m3 πx ( 0.0053 m )2
R= = =
µ µ −3 kg
1 x 10
m∙s

R=2001,8688

Los resultados indican que el fluido estudiado se mueve con un flujo de carácter laminar.

Discusión
El caudal es directamente proporcional al coeficiente integral de transferencia térmica debido a
que dicho coeficiente es mayor a medida que el caudal se incrementa.

Al realizar el cálculo del coeficiente integral se obtuvo que el coeficiente interno era mayor que el
externo, esto se debe a que, en un flujo en contracorriente, la transferencia de calor para el fluido
caliente va a ser mayor en el área interna que en el área externa de los tubos, por lo que se
cumple que el coeficiente integral interno es mayor que el externo. Esto también se debe al
diámetro, ya que este último determina el área de una manera inversamente proporcional al
coeficiente integral de transferencia térmica, de esta forma, al ser el diámetro interno menor que
el diámetro externo, el coeficiente interior será mayor que el exterior. Esta diferencia también se
podría deber a que en el caso del coeficiente externo se produce una leve pérdida de calor debido
a que la conductividad térmica varía.

Hay concordancia entre lo calculado, ya que hay una eficiencia mayor en los intercambiadores de
contracorriente comparado con el de paralelo, debido a que en contracorriente la temperatura de
salida del agua fría puede acercarse a la mayor temperatura del agua caliente y en cambio en los
intercambiadores de flujo paralelo hay una mayor diferencia de temperatura a la que entra en el
dispositivo.

Con respecto al cálculo del número de Reynolds se obtuvo que el fluido frio presenta un flujo del
tipo laminar, y el fluido caliente presenta un flujo del tipo turbulento. A mayor caudal, obtenemos
un mayor número de Reynolds. según la eficiencia calculada resulto más optimo que el sistema
funcione como calefactor. Es importante señalar que los resultados obtenidos podrían sufrir algún
grado de error, lo que puede conducir a una evidente interpretación errónea de los datos y así a
un consecuente índice de error considerable.

Conclusión
Se analizó el funcionamiento del equipo como sus partes. Se llevó a cabo el proceso de
transferencia de calor con el equipo configurado para flujo en paralelo y en contracorriente, y se
logró apreciar que la eficiencia del proceso era distinta, siendo mejor para la configuración en
contracorriente.

Se realizaron los cálculos de U y de Re a partir de los datos obtenidos en el proceso. Para U se

observó que los valores calculados que se asemejaban más a los entregados por el software eran
los del proceso en contracorriente con caudales 1/1, considerando solo el fluido caliente. A partir
de los Recalculados, se dedujo que los fluidos estudiados se movían con flujo laminar.

Bibliografía
MCCABE, Warren L., SMITH, Julian C., HARRIOTT, Peter, Operaciones unitarias en
ingeniería química. 7ª ed. Ciudad de méxico, McGraw-Hill, 2007. 311p.

2. SERWAY, R., JEWETT, J., Física para ciencias e ingeniería. Vol. 1. 7ªed.,MÉXICO, D.F.,
Cengage Learning Editores,S.A. de C.V., 2008. 554p.