Diseño de equipos de almacenamiento e

intercambiadores de calor
“Remoción de CO​2​”
​ 20122)
Proyecto de Ingeniería ​(4
Ingeniería Civil Química

DEPARTAMENTO DE INGENIERÍA EN MADERAS

Facultad de Ingeniería

Profesor: Dr. Luis E. Arteaga P.

Alumnos: ​ Matías Alarcón S.

Samuel Cárdenas G.
Cristian Espejo E.
Jorge Herrera S.
Marcelo Inzunza V.
Rodrigo Regla M.

Concepción, Viernes 12 de Octubre de 2018

Índice

Resumen 2

Objetivos 3

Descripción y parámetros de diseño de intercambiadores de calor: 3

Tanques de almacenamiento 7

Conclusión 7

Referencias 8

Anexo 9

1
Resumen

El proceso de absorción y desorción de CO​2 necesita condiciones específicas para una

correcta operación. Una de ellas es brindar una temperatura adecuada a los flujos de
productos que llegan a las torres de absorción y desorción por medio de intercambiadores
de calor, los cuales proporcionan el requisito energético necesario para que se produzcan
las operaciones unitarias. La Figura 1 muestra en rojo los equipos de intercambio de calor
presente en la remoción de CO​2​. Para poder dimensionar los intercambiadores de calor, se
basa en el método de Kern [1].

Figura 1:​​ Diagrama de flujo del proceso de remoción de CO​2

2
Objetivos

● Desarrollar un método de cálculo para el diseño de los intercambiadores de calor

propuestos de acuerdo a la norma TEMA.
● Determinar parámetros fundamentales en el diseño de los intercambiadores de
calor, tales como la superficie de intercambio, dimensión y diseño de tubos y
corazas.

Descripción y parámetros de diseño de intercambiadores de calor:

Diseño de intercambiador W-302

Recipiente: ​W-302

Descripción: ​El equipo W-302 de la Figura 1, conocido como intercambiador

amina-amina, busca elevar la temperatura de la amina enriquecida en CO​2 entrante a la
torre de desorción, aprovechando el calor de la amina pobre en CO​2 destinada a
recirculación. La amina rica en CO​2 debe ingresar al regenerador a una temperatura de
90°C y la amina que sale de la torre de desorción tiene una temperatura que oscila entre
110°C y 116°C.
Debido al cruce de temperatura, se tienen dos intercambiadores en serie (W1 y W2),
visibles en la Figura 4 del Anexo. Las características y dimensiones se aprecian en las
siguientes tablas. La determinación de estas variables se encuentran como ejemplo de
cálculo en el Anexo.

W1 W2

Según norma TEMA AES AES

Tipo de tubo 1 ¼” 1 ¾”

BWG 20 18

Diámetro interno Tubo 29,97 mm 42,02 mm

Diámetro externo Tubo 31,75 mm 44,5 mm

3
 Largo tubos 7,32 m 6,1 m

Número de tubos 3348 2158,94 ≈ 2159 tubos

Material de tubos Acero inoxidable 316 Acero Inoxidable 316

Pitch 23,8125 mm 55,625 mm

Arreglo de tubos Cuadrado 90° Cuadrado 90°

Diámetro interno carcasa 2500 mm 2893,09 mm

Número de pasos 1:2 1:2

Carcasa:Tubos

Tipo de deflector Segmentación único Segmentación único

Espaciado entre deflectores 1500 mm 2893,09 mm

Corte de deflector 25% 25%

Fluido por tubos Corriente fría de MDEA Corriente fría de MDEA

pobre en CO​2 pobre en CO​2

Flujo másico por tubo 707,242 kg/s 707,242 kg/s

Fluido por coraza Corriente caliente de MDEA Corriente fría de MDEA

enriquecida en CO​2 enriquecida en CO​2

Flujo másico por coraza 652,197 Kg/s 652,197 Kg/s

Velocidad fluido por tubos 1,22 m/s 1,83 m/s

Velocidad fluido por 0,9 m/s 1,05 m/s

carcasa

Tamaño de boquillas 10” 10”

Caída de presión en tubos 0,071 bar 0,052 bar

Caída de presión en coraza 0,47 bar 0,44 bar

Tabla 1:​​ Características intercambiador W-302, W1 y W2.

● Se seleccionó la configuración AES. Cabezal extremo frontal tipo A (canal y

cubierta desmontable), coraza tipo E (coraza de un paso) y cabezal extremo
posterior S (cabezal flotante con dispositivo de apoyo). Debido a que es el modelo

4
 más común y permite una configuración de un paso por carcasa y dos pasos por
tubos que es el que se dimensionó, incorporación de deflectores y permite
desmontaje para limpieza [2].
● Se consideró un arreglo cuadrado de tubos en 90°, dado que el fluido que circula
tiene acción corrosiva e incrustante. Este tipo de arreglo es recomendado cuando
la limpieza mecánica es crítica, cuando el flujo está en régimen turbulento y para
minimizar la caída de presión.
● Este intercambiador es con flujo a contracorriente. Se asignó la corriente de
MDEA enriquecida en CO​2​ por criterio de corrosión, para facilitar limpieza [2].
● La apertura de las boquillas es función del diámetro interno, de acuerdo a la
Figura 16, para diámetro interno de coraza mayor a 39 pulg se recomiendan
boquillas de 10 pulg.

Diseño de intercambiador W-301

Recipiente: ​W-301

Descripción​​: El equipo W-301 de la Figura 1 tiene como principal función extraer calor de
la corriente para satisfacer la necesidad térmica de la torre de absorción, la cual está
fijada a 35°C, sin embargo, antes de entrar al intercambiador la corriente 10 lleva una
temperatura de 75,3°C, motivo por el cual se busca disminuir su temperatura por medio
de un intercambiador con uso de 700 Kg/s de agua de servicio a 23°C a 1 bar, como
refrigerante. Tras análisis de perfil de temperatura, se experimentó cruce de
temperatura, debido a aquello especifica los siguientes intercambiadores en serie.

U1 U2

Según norma TEMA AES AES

Tipo de tubo 1 ¼” 2”

BWG 18 11

Diámetro interno Tubo 29,26 mm 44,70 mm

Diámetro externo Tubo 31,75 mm 50,80 mm

Largo tubos 6,1 m 6,1 m

Número de tubos 3471 2140

Material de tubos Acero inoxidable 316 Acero Inoxidable 316

5
 Pitch 39,69 mm 63,50 mm

Arreglo de tubos Triangular 30° Triangular 30°

Diámetro interno carcasa 2400 mm 3079 mm

Número de pasos 1:2 1:2

Carcasa:Tubos

Tipo de deflector Segmentación único Segmentación único

Espaciado entre deflectores 1454 mm 1864 mm

Corte de deflector 25% 25%

Fluido por tubos Corriente de MDEA-H20 Corriente de MDEA-H20

Flujo másico por tubo 66,0319 kg/s 66,0319 kg/s

Fluido por coraza Agua de servicio Agua de servicio

Flujo másico por coraza 700 kg/s 700 kg/s

Velocidad fluido por tubos 0,58 m/s 0,4 m/s

Velocidad fluido por carcasa 0,99 m/s 0,6 m/s

Tamaño de boquillas 10” 10”

Caída de presión en tubos 0,4 bar 0,16 bar

Caída de presión en coraza 0,0018 bar 0,0005 bar

Tabla 2:​​ Características intercambiador W-301, U1 y U2

● Se seleccionó la configuración AES. Cabezal extremo frontal tipo A (canal y

cubierta desmontable), coraza tipo E (coraza de un paso) y cabezal extremo
posterior S. Debido a que es el modelo más común y permite una configuración de
un paso por carcasa y un pasos por tubos que es el que se dimensionó,
incorporación de deflectores y permite desmontaje para limpieza [2].
● Se consideró un arreglo triangular 30°, dado que el fluido que por la carcasa es
agua limpia en un régimen turbulento. La elección de configuración triangular se
sugiere cuando hay gran cantidad de tubos y se desea disminuir el tamaño de la
carcasa [2].
● Este intercambiador es con flujo a contracorriente.
● La apertura de las boquillas es función del diámetro interno, de acuerdo a la Figura
16, para diámetro interno de coraza mayor a 39 pulg se recomiendan boquillas de
10 pulg.

6
Tanques de almacenamiento

Los tanques de almacenamiento o depósitos tienen como función principal el

almacenamiento inicial para materias primas, corrientes intermedias o finales
correspondientes al producto final de los fluidos del proceso. Para nuestra sección de
planta en particular, no se necesita la utilización de estos depósitos, debido a que los
productos principales como el gas dulce en la torre absorbedora y el CO​2 en la torre
regeneradora, son líneas de productos que se destinan a otras secciones de una planta de
amoniaco. Además, los flujos que circulan en la sección de remoción de CO​2 son
aprovechados energéticamente en los equipos de intercambio de calor, para luego ser
transportados con las condiciones necesarias en donde se requieran.

Conclusión

Se concluye que una metodología efectiva para el dimensionamiento de intercambiadores

de tubo y coraza es el método de kern, seleccionando configuraciones de cabezales y
coraza de acuerdo a la norma TEMA. En la determinación de los parámetros
fundamentales en los intercambiadores W-301 y W-302, se determinó que dada las
condiciones de operación, se experimenta cruce de temperatura, es por esto que se hace
necesario una nueva configuración de dos intercambiadores de calor de tubos y coraza
conectados en serie, en ambos casos. Sin embargo, las magnitudes determinadas escapan
de los estándares comerciales, esto se debe a la gran cantidad de área de intercambio
requerida, que exceden los 1800 [m​2​]. Los arreglos comerciales no lograron suplir esta
demanda, conllevando al diseño con proporciones poco convencionales.

7
Referencias

[1] R.K.Sinnott, Coulson & Richardson’s Chemical Engineering Design,

ELSEVIER Coulson Richardson’s Chem. Eng. Ser. 6 (2005) 634–709.
[2] INTERCAMBIADORES DE CALOR. (TEMA 7). Material complementario
Adecca.
[3] TEMA. (2007). Standards of the Tubular Exchanger Manufacturers
Association. New York.

8
 ANEXO

1.1 Marco Teórico

Para poder determinar el área de transferencia, se procede del balance de energía. Donde
Th es el fluido frio, y Tc es el fluido caliente.

𝑑𝐸 𝑑( 𝑈 + 𝐸𝑐 + 𝐸𝑝 )
= ̇ (𝐸̂𝑐,𝑖𝑛 + 𝐸̂𝑝,𝑖𝑛 + 𝐻
= 𝑚𝑐,𝑖𝑛 ̂𝑖𝑛 ) − 𝑚𝑐,𝑜𝑢𝑡 ̂𝑜𝑢𝑡 ) − 𝑄̇ + 𝑊̇ + ∆𝐻
̇ (𝐸̂𝑐,𝑜𝑢𝑡 + 𝐸̂𝑝,𝑜𝑢𝑡 + 𝐻 ̂𝑟𝑥
𝑑𝑡 𝑑𝑡

Considerando las siguientes suposiciones:

𝑑𝐸
- Estado estable 𝑑𝑡 = 0
- ∆𝐸̂𝑐 = ∆𝐸̂𝑝 = 0
- No hay trabajo 𝑊̇ = 0
- Conservación de masa 𝑚𝑐,𝑖𝑛 = 𝑚𝑐,𝑜𝑢𝑡 = 𝑚𝑐
- No hay reacción ∆𝐻 ̂𝑟𝑥 = 0
- El calor intercambiado es negativo, ya que es cedido al fluido frio −𝑄̇
̂𝑐 + 𝑄̇
0 = 𝑚̇ 𝑐 ∆𝐻

Si no existe cambio de fase, de acuerdo con las condiciones de operación y naturaleza de los
componentes, entonces la variación de entalpia se encuentra basada en el calor sensible:
𝑇2
̂ 𝑐 ∙ 𝑑𝑇
̂𝑐 = ∫ 𝐶𝑝
∆𝐻
𝑇1

De esta manera, el calor neto intercambiado para el fluido caliente es:

𝑇𝑜𝑢𝑡
𝑄̇ = − 𝑚 ̂𝑐 = − 𝑚
̇ 𝑐 ∆𝐻 ̇ 𝑐∫ ̂ 𝑐 ∙ 𝑑𝑇
𝐶𝑝
𝑇𝑖𝑛

De manera análoga para el fluido frio con recepción de calor + 𝑄̇, se tiene que:

9
 𝑇𝑜𝑢𝑡
𝑄̇ = 𝑚̇ ℎ ∆𝐻
̂ℎ = 𝑚̇ 𝑓 ∫ ̂ ℎ ∙ 𝑑𝑇
𝐶𝑝
𝑇𝑖𝑛

Considerando que el calor que entrega el fluido frio es el mismo que el caliente, se tiene
que:

𝑇𝑜𝑢𝑡 𝑇𝑜𝑢𝑡
𝑄̇ = 𝑚̇ ℎ ∫ ̂ ℎ ∙ 𝑑𝑇 = − 𝑚
𝐶𝑝 ̇ 𝑐∫ ̂ 𝑐 ∙ 𝑑𝑇
𝐶𝑝
𝑇𝑖𝑛 𝑇𝑖𝑛

Por otro lado, se tiene que la ecuación que gobierna a un intercambiador de calor es la
siguiente

𝑄̇ = 𝑈 𝐴𝑠 ∆𝑇𝑚𝑙 𝐹

𝑈: 𝐶𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑡𝑜𝑡𝑎𝑙 𝑑𝑒 𝑡𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟

𝐴𝑠: Á𝑟𝑒𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑎𝑙

∆𝑇𝑚𝑙: 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑚𝑒𝑑𝑖𝑎 𝑙𝑜𝑔𝑎𝑟í𝑡𝑚𝑖𝑐𝑎
𝐹: 𝐹𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑐𝑜𝑟𝑟𝑒𝑐𝑐𝑖ó𝑛 (𝐴𝑝𝑙𝑖𝑐𝑎𝑏𝑙𝑒 𝑎 𝑇𝑢𝑏𝑜 𝑦 𝑐𝑜𝑟𝑎𝑧𝑎𝑠 𝑑𝑒 𝑝𝑎𝑠𝑜𝑠 𝑚𝑢𝑙𝑡𝑖𝑝𝑙𝑒𝑠 𝑦 𝑓𝑙𝑢𝑗𝑜 𝑐𝑟𝑢𝑧𝑎𝑑𝑜)

De esta manera la ecuación que gobierna al intercambiador de calor es la siguiente.

𝑇𝑜𝑢𝑡 𝑇𝑜𝑢𝑡
𝑄̇ = 𝑈 𝐴𝑠 ∆𝑇𝑚𝑙 𝐹 = 𝑄̇ = 𝑚̇ ℎ ∫𝑇𝑖𝑛 𝐶𝑝
̂ ℎ ∙ 𝑑𝑇 = − 𝑚 ̂ 𝑐 ∙ 𝑑𝑇
̇ 𝑐 ∫𝑇𝑖𝑛 𝐶𝑝

De la ecuación anterior se puede determinar el área superficial por medio del balance del
fluido frio o caliente.
Donde el coeficiente total de transferencia de calor U se determina por medio de una
ecuación que cuantifica las resistencias:
𝑑
1 1 1 𝑑0 ln ( 𝑜 ) 𝑑 1 𝑑0 1
𝑑𝑖 0
= + + + ∙ + ∙
𝑈 ℎ0 ℎ𝑜𝑑 2 ∙ 𝑘𝑤 𝑑𝑖 ℎ𝑖𝑑 𝑑𝑖 ℎ𝑖
Donde
𝑈: Coeficiente global basado en el área exterior del tubo, W/m2 °C.
ℎ0 : Coeficiente de la película del fluido exterior W/m2 °C.
ℎ𝑖 : Coeficiente de la película del fluido interior W/m2 °C.

10
ℎ𝑜𝑑 : Coeficiente de suciedad exterior (factor de ensuciamiento), W/m2 °C.
ℎ𝑖𝑑 : Coeficiente de suciedad interior, W/m2 °C.
𝑘𝑤 : Conductividad térmica del material de la pared del tubo, W/m2 °C.
𝑑𝑖 : Diámetro interior del tubo, m.
𝑑𝑜 : Diámetro exterior del tubo, m.

Capacidad Calorífica a presión constante de una mezcla

El fluido caliente es una mezcla compuesta por MDEA, piperazina y otros componentes, que
pasa de T hasta T’, donde como mezcla permanece en estado líquido bajo la presión de 3
bar.
𝑇𝑜𝑢𝑡
𝑄̇ = − 𝑚
̇ 𝐶∫ ̂ 𝐶 ∙ 𝑑𝑇
𝐶𝑝
𝑇𝑖𝑛

Se tiene que, para una mezcla, la capacidad calorífica es la ponderada.

𝑇𝑜𝑢𝑡
𝑄̇ = −𝑚̇ 𝐶 ∫ ̂ 1 + 𝑋2 𝐶𝑝
( 𝑋1𝐶𝑝 ̂ 3 + ⋯ . +𝑋𝑛 𝐶𝑝
̂ 𝑛 ) ∙ 𝑑𝑇
𝑇𝑖𝑛

Donde a su vez, cada capacidad calorífica se puede extraer de tablas, nomogramas o por
medio de ecuaciones polinómicas con coeficientes factibles en un intervalo de tiempo de la
forma:
𝐶𝑝 = 𝐴 + 𝐵𝑇 + 𝐶𝑇 2 + 𝐷𝑇 3

2.- Ejemplo de cálculo intercambiador W-302

Este intercambiador de tubo y coraza cumple con la función de enfriar la amina MDEA
pobre en CO2 que sale por el fondo de la torre de desorción y transferir el calor a la amina
entrante a la torre de desorción y para realizar los cálculos de diseño es necesario tener
en cuenta la ecuación donde se obtienen los parámetros del intercambiador de calor:

Q=U* A* LMTD* FT
Para ello se tienen que obtener las propiedades de los distintos fluidos.

Corriente Flujo másico [kg/s] Temperatura de entrada [°C] Temperatura de salida [°C]
Amina tubo 707,243 55 99
Amina coraza 652,196 138 86
Tabla 3: Flujos y temperaturas de ambas corrientes.

Los datos de la tabla anterior se obtuvieron a través del simulador Aspen Plus, Sin
embargo, el perfil de temperatura experimenta un cruce, que es el que se aprecia en la
Figura 3:

11
 Figura 3: Perfil de temperatura en el intercambiador de tubo y coraza inicial

Por la razón anterior, se decidió considerar dos intercambiadores de coraza conectados

en serie, para evitar este cruce. Esta modificación es visible en la Figura 4, donde los
perfiles que se experimentan en cada intercambiador se aprecian en la Figura 5 y 6.

Figura 4: Conexión en serie de los intercambiadores de calor W1 y W2.

12
 Figura 5: Perfil de temperatura en Intercambiador W1.

Figura 6: Perfil de temperatura del intercambiador W2.

13
2.1- Dimensionamiento intercambiador W1.
El procedimiento para poder dimensionar el intercambiador de calor es el descrito por el
libro Diseño en Ingeniería Química. Ray Sinnott, Gavin Towler. Que se sustenta en el
método de Kern.

Figura 7: Procedimiento para dimensionar un intercambiador.

Para comenzar a dimensionar se debe saber que, de los balances energéticos, se puede
determinar el calor transferido del fluido caliente al frío, independiente de cual se analice,
según lo visto en el marco teórico. De esta manera, se procede a trabajar con la corriente
MDEARE2 y MDEARE1, para esto se deben conocer las temperaturas de las corrientes, los
flujos másicos y algunas propiedades termodinámicas, las que son extraídas por el simulador
Aspen Plus.

14
 MDEARE1 MDEARE2 RICHAW1 RICHAW2
T °C 112,585 86,116 54,855 80
Flujo másico [kg/s] 652,197 652,197 707,242 707,242
Presión [bar] 3 3 3 3
Cp [J/kg K] 3373,47 3308,59 - -
Tabla 4: Condiciones de operación del intercambiador W1.

𝑇𝑜𝑢𝑡
De la ecuación 𝑄 = − 𝑚 ̂ 𝑐 ∙ 𝑑𝑇, se puede considerar un promedio en la
̇ 𝑖 ∫𝑇𝑖𝑛 𝐶𝑝
capacidad calorífica debido a que la corriente es una mezcla y para la determinación por
polinomio se necesitan los valores de los coeficientes del polinomio, que se dificulta
encontrarlos por las sustancias que se manejan como lo son el MDEA y la piperazina,
mencionando que además, dependen de las condiciones de operación. Por lo anterior, se
tiene:
𝐶𝑝𝑀𝐷𝐸𝐴𝑅𝐸1 + 𝐶𝑝𝑀𝐷𝐸𝐴𝑅𝐸2 𝐽
𝐶𝑝𝑚𝑒𝑧𝑐𝑙𝑎 = = 3341,03 [ ]
2 𝑘𝑔 𝑘
De esta manera el calor se determina resultando en:
𝐾𝑔 𝐽
𝑄 = −6207,753 [ ] ∙ 3341,03 [ ] ∙ (86,116 − 112,585)°𝐶 = 57676238 [𝑤]
𝑠 𝑘𝑔 𝐾

Luego se toma como valor de coeficiente global U, el que arroja el simulador Aspen:

Figura 8: Resultados de la simulación en Aspen Plus.

Es así como se trabaja con un coeficiente global que incluye el ensuciamiento equivalente
a 850 [w/m2 k].

15
Sin embargo, para determinar el área, se necesita determinar la temperatura media
logarítmica y el factor de corrección F, el que puede ser extraído de gráficos según la
configuración del intercambiador de calor a través de gráficos.

∆𝑇𝑎 − ∆𝑇𝑏 (112,585 − 80) − (86,116 − 54,855)

∆𝑇𝑚𝑙 = = = 31,92°𝐶
∆𝑇𝑎 112,585 − 80
𝑙𝑛 ( ) 𝑙𝑛 ( )
∆𝑇𝑏 86,116 − 54,855

Figura 9: Factor de corrección de temperatura para un paso en coraza y 2 o n pasos par

en tubo.

Ingresando con los parámetros S y R a través de los siguientes cálculos:

𝑡2 − 𝑡1 80 − 54,855
𝑆= = = 0,435
𝑇1 − 𝑡1 112,585 − 54,855
𝑡2 − 𝑡1 112,585 − 86,116
𝑅= = = 1,052
𝑇1 − 𝑡1 80 − 54,855
𝐹 = 0,87

De esta manera se puede determinar el área:

57676230 [𝑤]
𝐴= 𝑤 = 2443,4 [𝑚2 ]
850 [𝑚2 𝐾 ] ∙ 31,92 𝑘 ∙ 0,87

16
Podemos ver que el área encontrada es consecuente con el área determinada por el
software, Aspen Plus.

Para cumplir con los parámetros determinados, se ha seleccionado las siguientes

dimensiones de tubos:

Longitud [m] Diámetro externo (do) BWG Diámetro Interno (di)

[pulg] [pulg]
7,32 1 ¼” 20 1,18
Tabla 5: Dimensiones de los tubos del intercambiador W1.

Las siguientes ecuaciones son extraídas del libro Diseño en Ingeniería Química. Ray
Sinnott, Gavin Towler del Sinnot.

• Número de tubos:

𝐴𝑡 𝐴𝑡 2443,4 [𝑚2 ]
𝑁𝑡 = = = = 3348 𝑡𝑢𝑏𝑜𝑠
𝐴𝑡𝑢𝑏𝑜 𝜋 ∙ 𝑑𝑜 ∙ 𝐿 0,7298 [𝑚2 ]

• Diámetro haz de tubos:

1 1
𝑁𝑡 𝑛 3348 2,291
𝐷𝑏 = 𝑑𝑜 ∙ ( ) = 31,75 [𝑚𝑚] ∙ ( ) = 2468,73 [𝑚𝑚]
𝐾𝑖 0,156

Los valores de n y Ki se extraen de la siguiente tabla, según sea el número de pasos y

arreglo de tubos que para este caso es un arreglo cuadrado con 2 pasos por tubo:

Figura 10: Valores de constantes para la determinación del diámetro del haz de tubos.

17
2.1.1 COEFICIENTE LADO DEL TUBO

Para poder determinar el coeficiente de convección asociado al fluido interior de los tubos,
se debe contar con algunas propiedades físicas y termodinámicas de las corrientes, a
temperatura media con ayuda del Software Aspen Plus.

54,855 + 80 °𝐶
𝑇𝑚 = = 67,4275°𝐶
2
𝐾𝑔
𝜌 = 492,558 [ ]
𝑚3
𝑤
𝑘 = 0,231679 [ ]
𝑚𝐾
𝑠
𝜇 = 0,000627754 [𝑁 ]
𝑚2
𝐽
𝐶𝑝 = 3078,12 [ ]
𝑘𝑔 𝐾

• Área total para el flujo en tubos:

𝑁𝑇 𝜋 3348 𝜋
𝐴𝑓 = ∙ ∙ 𝑑𝑖 2 = ∙ ∙ (29,97𝑥10−3 [𝑚])2 = 1,18 [𝑚2 ]
𝑁𝑃𝑎𝑠𝑜𝑠 4 2 4

• Velocidad del flujo:

𝑘𝑔
𝑚̇ 707,242 [ 𝑠 ] 𝑚
𝑣= = = 1,22 [ ]
𝜌 ∙ 𝐴𝑓 𝑘𝑔 𝑠
492,558 [ 3 ] 1,18 [𝑚2 ]
𝑚
• Número de Reynolds:
𝜌 𝑣 𝑑𝑖
𝑅𝑒 = = 28689
𝜇
• Número de Prandtl:
𝐶𝑝 𝜇
Pr = = 8,34
𝑘𝑓

• Factor de transferencia de calor 𝐽ℎ : Se puede omitir el termino de corrección de

viscosidades, ya que cuando la viscosidad del fluido es muy baja, como es el caso,
dicha relación tiende a uno.
𝜇 −0,14
𝐽ℎ = 𝑆𝑡 𝑃𝑟 0,67 ( ) = 3,46 𝑥 10−3
𝜇𝑤

𝐸 = 0,0225𝐸𝑥𝑝(−0,0225(𝑙𝑛𝑃𝑟)2 )

18
 𝑆𝑡 = 𝐸 𝑅𝑒 −0,205 𝑃𝑟 −0,505 (𝑠𝑜𝑙𝑜 𝑒𝑛 𝑅𝑒 > 10 000)
• Determinación de coeficiente ℎ𝑖 :
𝐾𝑓 𝜇 0,14 𝑤
ℎ𝑖 = 𝐽ℎ 𝑅𝑒 𝑃𝑟 0,33 ( ) = 1545 [ 2 ]
𝑑𝑖 𝜇𝑤 𝑚 𝑘

2.1.2 COEFICIENTE DEL LADO DE LA CARCASA

Al igual que en tubos, se debe contar con algunas propiedades físicas y termodinámicas
de las corrientes, a temperatura media con ayuda del Software Aspen Plus.

86,116 + 112,585 °𝐶
𝑇𝑚 = = 99,3505°𝐶
2
𝐾𝑔
𝜌 = 968,928 [ ]
𝑚3
𝑤
𝑘 = 0,233669 [ ]
𝑚𝐾
𝑠
𝜇 = 0,366425𝑥10−3 [𝑁 ]
𝑚2
𝐽
𝐶𝑝 = 3339,32 [ ]
𝑘𝑔 𝐾

• Pitch: Se tiene un arreglo cuadrado 90°, la siguiente es la distancia entre centros.

𝑃𝑡 = 1,25 ∙ 𝑑𝑜 = 1,25 ∙ 31,75 [𝑚𝑚] = 39,6875 [𝑚𝑚]

• Espaciado deflectores: Tras realizar dos iteraciones previas, para evitar tener una
caída de presión significativa y tener velocidad de flujo entre 1-2 [m/s], se estima
lo siguiente:
3
𝑙𝐵 = 𝐷𝑠 = 1500[mm]
5

• Diámetro interno de carcasa: debido a el diámetro del haz 𝐷𝑏 , se debió extrapolar

la recta de la Figura 15 para determinar el espaciado.

𝐷𝑠 = 𝐷𝑏 + 𝑒𝑠𝑝𝑎𝑐𝑖𝑎𝑑𝑜 = 2468,73 + 32,69 = 2501,42 [𝑚𝑚]

• Área total para el flujo cruzado:

(𝑃𝑡 − 𝑑𝑜) 𝐷𝑆 𝑙𝐵
𝐴𝑆 = = 0,75 [𝑚2 ]
𝑃𝑡

19
 • Diámetro equivalente:
1,27
𝑑𝑒𝑞 = (𝑃𝑡 2 − 0,785 𝑑𝑜2 ) = 0,03135 [𝑚]
𝑑𝑜

• Velocidad del flujo:

𝑘𝑔
𝑚̇ 652,197 [ 𝑠 ] 𝑚
𝑣= = = 0,9 [ ]
𝜌 ∙ 𝐴𝑠 𝑘𝑔 𝑠
968,928 [ 3 ] 0,5 [𝑚2 ]
𝑚
Esta cercano al límite inferior, pero es factible.
• Número de Reynolds:
𝜌 𝑣 𝑑𝑖
𝑅𝑒 = = 71324
𝜇
• Número de Prandtl:
𝐶𝑝 𝜇
𝑃𝑟 = = 5,2365
𝑘𝑓

• Factor de transferencia de calor 𝐽ℎ : Este término se extrae del siguiente gráfico,

considerando el número de Reynolds y el % de corte del bafle (25%).

Figura 11: Factores de transferencia de calor del lado de la carcasa, deflectores

segmentados.

20
 𝐽ℎ = 2𝑥 10−3

• Determinación de coeficiente ℎ𝑜 : Omitiendo el factor de corrección de la

viscosidad, se tiene lo siguiente.

𝐾𝑓 𝜇 0,14 𝑤
ℎ𝑜 = 𝐽ℎ 𝑅𝑒 𝑃𝑟 0,33 ( ) = 1820 [ 2 ]
𝑑𝑒𝑞 𝜇𝑤 𝑚 𝑘

2.1.3 DETERMINACIÓN DE CAÍDA DE PRESIÓN

Los valores determinados deben cumplir con las especificaciones del apartado A.2 en la
sección de caída de presión.

• Factor de fricción: El valor de 𝐽𝑓 se determina de las siguientes Figuras, 12 y 13,

ingresando con el valor de Reynolds.

Figura 12: Factor de fricción del lado de los tubos.

21
 Figura 13: Factor de fricción del lado de la carcasa, deflectores segmentados.

• Caída de presión en tubos: Se omite el término de corrección de viscosidad. Se

determina por la Figura 6 que 𝐽𝑓 = 3,7𝑥10−3

𝐿 𝜇 −𝑚 𝜌 ∙ 𝑣2
∆𝑃 = 𝑁𝑝 ∙ (8 ∙ 𝐽𝑓 ∙ ( ) ∙ ( ) + 2,5) = 7133 [𝑃𝑎] = 7,13 [𝐾𝑃𝑎]
𝑑𝑖 𝜇𝑤 2
La caída de presión esta dentro del rango permisible.

• Caída de presión en coraza: Se omite el termino de corrección de viscosidad. Se

determina por la Figura 12 que 𝐽𝑓 = 3,9𝑥10−2

𝐷𝑠 𝐿 𝜇 −0,14 𝜌 ∙ 𝑣 2
∆𝑃 = 8 ∙ 𝐽𝑓 ∙ ( )∙( )∙( ) = 47673[𝑃𝑎] = 47 [𝑘𝑃𝑎]
𝑑𝑒𝑞 𝑙𝐵 𝜇𝑤 2
La caída de presión está sobre el limite superior, pero aun así es factible.

22
2.2- Dimensionamiento intercambiador W2.

MDEAHOT MDEARE1 RICHAW2 RICHAC2

T °C 137,651 112,585 80 98,89
Flujo másico [kg/s] 652,197 652,197 707,242 707,242
Presión [bar] 3 3 3 3
Cp [J/kg K] 3451,21 3373,47 - -
Tabla 6: Condiciones de operación del intercambiador W1.

𝐶𝑝𝑀𝐷𝐸𝐴𝑅𝐸1 + 𝐶𝑝𝑀𝐷𝐸𝐴𝑅𝐸2 𝐽
𝐶𝑝𝑚𝑒𝑧𝑐𝑙𝑎 = = 3412,34 [ ]
2 𝑘𝑔 𝑘
De esta manera el calor se determina resultando en:
𝐾𝑔 𝐽
𝑄 = − 652,197 [ ] ∙ 3412,34 [ ] ∙ (112,585 − 137,651)°𝐶 = 55555555,6 [𝑤]
𝑠 𝑘𝑔 𝐾

Luego se toma como valor de coeficiente global U, el que arroja el simulador Aspen:

Figura 14: Resultados de la simulación en Aspen Plus W2.

Es así como se trabaja con un coeficiente global que incluye el ensuciamiento equivalente
a 850 [w/m2 k].

23
 ∆𝑇𝑎 − ∆𝑇𝑏 (137,65 − 98,89) − (112,585 − 80)
∆𝑇𝑚𝑙 = = = 35,5°𝐶
∆𝑇𝑎 137,65 − 98,89
𝑙𝑛 ( ) 𝑙𝑛 ( )
∆𝑇𝑏 112,585 − 80

Ingresando a La Figura 7 con los parámetros S y R a través de los siguientes cálculos:

𝑡2 − 𝑡1 98,89 − 80
𝑆= = = 0,327
𝑇1 − 𝑡1 137,65 − 80
𝑡2 − 𝑡1 137,65 − 112,585
𝑅= = = 1,3268
𝑇1 − 𝑡1 98,89 − 80
𝐹 = 0,95

De esta manera se puede determinar el área:

55555555,6 [𝑤]
𝐴= 𝑤 = 1841,112 [𝑚2 ]
850 [𝑚2 𝐾 ] ∙ 35,5 𝑘 ∙ 0,95
Podemos ver que el área encontrada es consecuente con el área determinada por el
software, Aspen Plus.

Para cumplir con los parámetros determinados, se ha seleccionado las siguientes

dimensiones de tubos:

Longitud [m] Diámetro externo (do) BWG Diámetro Interno (di)

[pulg] [pulg]
6,1 1 3⁄4” 18 1,6543
Tabla 7: Dimensiones de los tubos del intercambiador W1.

Las siguientes ecuaciones son extraídas del libro Diseño en Ingeniería Química. Ray
Sinnott, Gavin Towler.

• Número de tubos:

𝐴𝑡 𝐴𝑡 1841,112 [𝑚2 ]
𝑁𝑡 = = = = 2158,94 𝑡𝑢𝑏𝑜𝑠 ≈ 2159𝑡𝑢𝑏𝑜𝑠
𝐴𝑡𝑢𝑏𝑜 𝜋 ∙ 𝑑𝑜 ∙ 𝐿 0,852353 [𝑚2 ]

• Diámetro haz de tubos:

1 1
𝑁𝑡 𝑛 2159 2,291
𝐷𝑏 = 𝑑𝑜 ∙ ( ) = 44,5 [𝑚𝑚] ∙ ( ) = 2857,09[𝑚𝑚]
𝐾𝑖 0,156

24
Los valores de n y Ki se extraen de la Figura 9, un paso por carcasa y dos pasos por tubo.

2.2.1 COEFICIENTE LADO DEL TUBO

Para poder determinar el coeficiente de convección asociado al fluido interior de los tubos,
se debe contar con algunas propiedades físicas y termodinámicas de las corrientes, a
temperatura media con ayuda del Software Aspen Plus.

98,89 + 80 °𝐶
𝑇𝑚 = = 89,445°𝐶
2
𝐾𝑔
𝜌 = 257,142 [ ]
𝑚3
𝑤
𝑘 = 0,231932 [ ]
𝑚𝐾
𝑠
𝜇 = 0,000440645 [𝑁 ]
𝑚2
𝐽
𝐶𝑝 = 3060,09 [ ]
𝑘𝑔 𝐾

• Área total para el flujo en tubos:

𝑁𝑇 𝜋 2159 𝜋
𝐴𝑓 = ∙ ∙ 𝑑𝑖 2 = ∙ ∙ (42,02𝑥10−3 [𝑚])2 = 1,5 [𝑚2 ]
𝑁𝑃𝑎𝑠𝑜𝑠 4 2 4

• Velocidad del flujo:

𝑘𝑔
𝑚̇ 707,242 [ 𝑠 ] 𝑚
𝑣= = = 1,83 [ ]
𝜌 ∙ 𝐴𝑓 𝑘𝑔 𝑠
257,142 [ 3 ] 1,5 [𝑚2 ]
𝑚

• Número de Reynolds:
𝜌 𝑣 𝑑𝑖
𝑅𝑒 = = 44874
𝜇
• Número de Prandtl:
𝐶𝑝 𝜇
Pr = = 5,814
𝑘𝑓

25
 • Factor de transferencia de calor 𝐽ℎ : Se puede omitir el termino de corrección de
viscosidades, ya que cuando la viscosidad del fluido es muy baja, como es el caso,
dicha relación tiende a uno.
𝜇 −0,14
𝐽ℎ = 𝑆𝑡 𝑃𝑟 0,67 ( ) = 3𝑥 10−3
𝜇𝑤

𝐸 = 0,0225𝐸𝑥𝑝(−0,0225(𝑙𝑛𝑃𝑟)2 )

𝑆𝑡 = 𝐸 𝑅𝑒 −0,205 𝑃𝑟 −0,505 (𝑠𝑜𝑙𝑜 𝑒𝑛 𝑅𝑒 > 10 000)

• Determinación de coeficiente ℎ𝑖 :
𝐾𝑓 𝜇 0,14 𝑤
ℎ𝑖 = 𝐽ℎ 𝑅𝑒 𝑃𝑟 0,33 ( ) = 1328,3 [ 2 ]
𝑑𝑖 𝜇𝑤 𝑚 𝑘

2.2.2 COEFICIENTE DEL LADO DE LA CARCASA

Al igual que en tubos, se debe contar con algunas propiedades físicas y termodinámicas
de las corrientes, a temperatura media con ayuda del Software Aspen Plus.

137,651 + 112,585 °𝐶
𝑇𝑚 = = 125,118°𝐶
2
𝐾𝑔
𝜌 = 948,368 [ ]
𝑚3
𝑤
𝑘 = 0,230988 [ ]
𝑚𝐾
𝑠
𝜇 = 0,270224𝑥10−3 [𝑁 ]
𝑚2
𝐽
𝐶𝑝 = 3409,11 [ ]
𝑘𝑔 𝐾

• Pitch: Se tiene un arreglo cuadrado 90°, la siguiente es la distancia entre centros.

𝑃𝑡 = 1,25 ∙ 𝑑𝑜 = 1,25 ∙ 44,5 [𝑚𝑚] = 55,625 [𝑚𝑚]

• Espaciado deflectores: Tras realizar dos iteraciones previas, para evitar tener una
caída de presión significativa y tener velocidad de flujo entre 1-2 [m/s], se estima
lo siguiente:
3
𝑙𝐵 = 𝐷𝑠 = 1157,236 [mm]
5

• Diámetro interno de carcasa: debido a el diámetro del haz 𝐷𝑏 , se debió extrapolar

la recta de la Figura 15 para determinar el espaciado.

26
 𝐷𝑠 = 𝐷𝑏 + 𝑒𝑠𝑝𝑎𝑐𝑖𝑎𝑑𝑜 = 2857,09 + 36 = 2893,09 [𝑚𝑚]

• Área total para el flujo cruzado:

(𝑃𝑡 − 𝑑𝑜) 𝐷𝑆 𝑙𝐵
𝐴𝑆 = = 0,6696 [𝑚2 ]
𝑃𝑡

• Diámetro equivalente:
1,27
𝑑𝑒𝑞 = (𝑃𝑡 2 − 0,785 𝑑𝑜2 ) = 0,04394 [𝑚]
𝑑𝑜

• Velocidad del flujo:

𝑘𝑔
𝑚̇ 652,197 [ 𝑠 ] 𝑚
𝑣= = = 1,05 [ ]
𝜌 ∙ 𝐴𝑠 𝑘𝑔 𝑠
948,368 [ 3 ] 0,75 [𝑚2 ]
𝑚
(Es válido, pero se trabaja en el límite)
• Número de Reynolds:
𝜌 𝑣 𝑑𝑒𝑞
𝑅𝑒 = = 158065
𝜇
• Número de Prandtl:
𝐶𝑝 𝜇
𝑃𝑟 = = 3,988 ≈ 4
𝑘𝑓

• Factor de transferencia de calor 𝐽ℎ : Este término se extrae de la Figura 11,

considerando el número de Reynolds y el % de corte del bafle (25%).

𝐽ℎ = 1,7 𝑥 10−3

• Determinación de coeficiente ℎ𝑜 : Omitiendo el factor de corrección de la

viscosidad, se tiene lo siguiente.

𝐾𝑓 𝜇 0,14 𝑤
ℎ𝑜 = 𝐽ℎ 𝑅𝑒 𝑃𝑟 0,33 ( ) = 2232 [ 2 ]
𝑑𝑒𝑞 𝜇𝑤 𝑚 𝑘

27
2.2.3 DETERMINACIÓN DE CAÍDA DE PRESIÓN
Los valores determinados deben cumplir con las especificaciones del apartado A.2 en la
sección de caída de presión.

• Factor de fricción: El valor de 𝐽𝑓 se determina de la Figura 12 y 13, según sea el

caso, ingresando con el valor de Reynolds.

• Caída de presión en tubos: Se omite el término de corrección de viscosidad. Se

determina por la Figura 12 que 𝐽𝑓 = 3,1𝑥10−3

𝐿 𝜇 −𝑚 𝜌 ∙ 𝑣2
∆𝑃 = 𝑁𝑝 ∙ (8 ∙ 𝐽𝑓 ∙ ( ) ∙ ( ) + 2,5) = 5253,14 [𝑃𝑎] = 5,2[𝐾𝑃𝑎]
𝑑𝑖 𝜇𝑤 2
La caída de presión está dentro del rango permisible.

• Caída de presión en coraza: Se omite el termino de corrección de viscosidad. Se

determina por la Figura 13 que 𝐽𝑓 = 3𝑥10−2

𝐷𝑠 𝐿 𝜇 −0,14 𝜌 ∙ 𝑣 2
∆𝑃 = 8 ∙ 𝐽𝑓 ∙ ( )∙( )∙( ) = 44263,36[𝑃𝑎] = 44[𝑘𝑃𝑎]
𝑑𝑒𝑞 𝑙𝐵 𝜇𝑤 2

La caída de presión excede el limite estimado, pero es aceptable.

3.- VARIABLES CONSIDERADAS EN EL DIMENSIONAMIENTO

El decidir el diámetro y el largo del tubo, influirá en la posterior determinación del

diámetro interno de carcasa, las velocidades de flujo y caída de presión.

• Diámetros de tubos: La siguiente tabla muestra las especificaciones para tubos de

acero comercial según BWG.

28
 Tabla 8: Diámetros de tubos de acero comercial.

• Longitud de tubos: Por otro lado, en cuando a largos de tubos, se indica que los
largos preferentes para intercambiadores son 1.83 [m], 2.44 [m], 3.66 [m], 4.88
[m], 6,10 [m] y 7,32 [m].

• Diámetro coraza: El estándar británico BS3274 abarca intercambiadores desde 6

in (150 mm) hasta 42 in (1067 mm). Por otro lado, los estándares TEMA
comprenden hasta los 1520 mm de diámetro. Sin embargo, existen compañías1
proveedoras de intercambiadores con carcazas de hasta 100 pulgadas (2.54 m).
Por otro lado, existen tablas que comprenden hasta 3 m de diámetro:

1
http://www.wermac.org/equipment/heatexchanger_part2.html

29
 Tabla 9: Dimensiones de intercambiadores de calor hasta 3 metros de diámetro.

Se debe considerar que el diámetro interno de la carcasa es igual al diámetro del haz de
tubos más un espaciado, la que puede ser estimada por la siguiente gráfica. En este caso,
es de interés el cabezal fijo y tubos en U.

30
 Figura 15: Espaciado en función del diámetro del haz para disposiciones de cabezales y
tubos.

• Deflectores: Según el libro Procesos de transferencia de calor de (Donald Q. Kern)

estima que para los deflectores segmentados, 25% corte, la distancia entre ellos
obedece a la siguiente relación en términos del diámetro interior de la coraza.
1
𝐷 < 𝑙𝐵 < 𝐷𝑖
5 𝑖

• Velocidad: En cuanto a la velocidad, se debe decir que el libro Diseño en Ingeniería

Química. Ray Sinnott, Gavin, recomienda para fluidos 1 – 2 [m/s], máximo 4 [m/s]
para evitar el ensuciamiento. En el caso del agua 1,5 a 2,5 [m/s] y en el lado de
carcasa 0,3 – 1 [m/s]. En el caso de vapores depende de la presión de trabajo y
densidad, por lo general

31
 Vacío 50 – 70 [m/s]
Presión atmosférica 10-30 [m/s]
Presión elevada 5 – 10 [m/s].

• Caída de Presión: Depende según la fase de la corriente, el libro Diseño en

Ingeniería Química. Ray Sinnott, Gavin, tiene las siguientes tablas de acuerdo a
criterios como la viscosidad en los líquidos, o a la presión a la que se encuentran
los gases y vapores.

Figura 16: Caída de presión permitida en el diseño según estado de fase de las corrientes.
• Ensuciamiento: Para la determinación del coeficiente global U de diseño, se debe
considerar el factor de ensuciamiento. La siguiente tabla expresa algunos valores
típicos.

Tabla 10: Valores de factor de ensuciamiento en distintos fluidos.

32
 Tabla 11: Valores de coeficiente U en intercambiador tubo y coraza según la naturaleza
de las corrientes.

• Tamaño Boquillas: la siguiente tabla es extraída del libro Process Heat Transfer,
D.Q. Kern.

Tabla 12: Tamaño de boquillas de acuerdo con el tamaño de coraza.

33