Memorias de Calculo Grupo 2

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 19

UNIVERSIDAD DE ANTIOQUIA

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Omar Xavier Fonseca Serrano, Karen Jaramillo Meneses, Paula Andrea Cespedes Grattz,
Laura Puerta Gonzales, Juan José Vidal Vargas.
Memorias de calculo
Fluido 1 (caliente): dietanolamina (45%9), propilbenceno (15%), estireno (40%).
Fluido 2 (frío): fenol (25%), alcohol bencílico (65%), diisopropil (10%).
29 de Julio del 2022
Cálculo de la temperatura de salida del fluido 2.
Inicialmente se usan las ecuaciones que se muestran a continuación para calcular el calor
transferido por cado fluido, luego considerando que el sistema está idealmente aislado se
usa la ecuación 3 para calcular la t2 con ayuda de la herramienta “buscar objetivo” de Excel.
𝑄 = 𝑊 ∗ 𝐶𝑝1 ∗ (𝑇1 − 𝑇2 ) 𝑬𝒒. 𝟏
𝑞 = 𝑤 ∗ 𝐶𝑝2 ∗ (𝑡1 − 𝑡2 ) 𝑬𝒒. 𝟐
𝑄 = 𝑞 𝑬𝒒. 𝟑
Donde,
𝑙𝑏
𝑊 = 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑐𝑎𝑙𝑖𝑒𝑛𝑡𝑒 [ ] = 1000000

𝑙𝑏
𝑤 = 𝐹𝑙𝑢𝑗𝑜 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟í𝑜 [ ] = 730000

𝑇1 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑓𝑙𝑢𝑖𝑑𝑜 𝑐𝑎𝑙𝑖𝑒𝑛𝑡𝑒 [°𝐹] = 298
𝑇2 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑓𝑙𝑢𝑖𝑑𝑜 𝑐𝑎𝑙𝑖𝑒𝑛𝑡𝑒 [°𝐹] = 197
𝑡1 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑑𝑒 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟í𝑜 [°𝐹] = 90
𝐵𝑇𝑈
𝐶𝑝1 = 𝐶𝑎𝑝𝑎𝑐𝑖𝑑𝑎𝑑 𝑐𝑎𝑙𝑜𝑟í𝑓𝑖𝑐𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑐𝑎𝑙𝑖𝑒𝑛𝑡𝑒, 𝑇𝑝𝑟𝑜𝑚 [ ] = 0.5189
𝑙𝑏 ∗ °𝐹
𝐵𝑇𝑈
𝐶𝑝2 = 𝐶𝑎𝑝𝑎𝑐𝑖𝑑𝑎𝑑 𝑐𝑎𝑙𝑜𝑟í𝑓𝑖𝑐𝑎 𝑑𝑒𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟í𝑜, 𝑇𝑝𝑟𝑜𝑚 [ ] = 0.5212
𝑙𝑏 ∗ °𝐹
Se obtiene que,
𝑡2 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎 𝑑𝑒 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟í𝑜 [°𝐹] = 227.7414
Es importante resaltar que los flujos fueron calculados de acuerdo con el NG
correspondiente para nuestro grupo, el cual es 2, es decir,
𝑙𝑏 𝑙𝑏
𝑊 = 𝑁𝐺 ∗ 500000 = 1000000 𝑬𝒒. 𝟒
ℎ ℎ
𝑙𝑏 𝑙𝑏
𝑤 = 2 ∗ 365000 = 730000 𝑬𝒒. 𝟓
ℎ ℎ
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

19 de agosto del 2022


Cálculo de los espesores de los tubos y camisa
Se escogió BWG 16 para tube y cedula 80.
Especificaciones de los tubos
Diámetro externo de los tubos: De acuerdo con las recomendaciones del profesor se
escogió un diámetro externo de 3/4 pulgadas. Luego, con la tabla 10 del libro de Kern y con
los datos proporcionados por el profesor el diámetro interno de los tubos (IDt) es de 0.62
pulgadas.

Figura 1. Elección del ID de acuerdo con la tabla 10 del libro de Kern.


Longitud de los tubos: De acuerdo con el módulo 4, las longitudes que más se usan son
de 8, 2, 15 y 20 ft. Se decidió utilizar una longitud de 20 ft.
Arreglo de tubos: El tipo de arreglo que se usará es triangular.
Pitch de los tubos: De acuerdo con el módulo 4, recomiendan para un diámetro exterior
de 3/4 in y un arreglo triangular un pitch de 15/16 in.
Espesor de los tubos. El espesor del tubo se determina por medio de la tabla 10 del libro
Kern y los datos iniciales proporcionados, 𝑡 = 0.065
Diferencia media logarítmica de temperatura LMTD.
Teniendo en cuenta que el arreglo del intercambiador es de flujo en contracorriente se
calcula el LMTD con la siguiente ecuación
(𝑇1 − 𝑡2 ) − (𝑇2 − 𝑡1 )
𝐿𝑀𝑇𝐷 = 𝑬𝒒. 𝟔
𝑇 −𝑡
ln (𝑇1 − 𝑡2 )
2 1
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

𝐿𝑀𝑇𝐷 = 87.345
Especificaciones para la camisa
Diámetro interno de la camisa. Este diámetro se calcula con la tabla 9 de Kern, tomando
el valor maximo. Teniendo en cuenta los valores obtenidos anteriormente, ODt, Pt y m se
lee el IDs. Se determina que el IDs es igual a 39 in.

Figura 2. Elección del ID de acuerdo con la tabla 9 del libro de Kern


Diámetro equivalente de la camisa. Se calcula con la ecuación 7.4 del libro de Kern.

𝑂𝐷 2
4 ∗ (𝑃𝑡2 − 𝜋 ∗ 4 𝑡 )
𝑑𝑒 = 𝑬𝒒. 𝟕
𝜋 ∗ 𝑂𝐷𝑡
Donde,
𝑑𝑒
𝐷𝑒 = 𝐷𝑖𝑎𝑚𝑒𝑡𝑟𝑜 𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑚𝑖𝑠𝑎 [𝑓𝑡] = = 0.04443
12
15
𝑃𝑡 = 𝑃𝑖𝑡𝑐ℎ [𝑖𝑛] =
16
3
𝑂𝐷𝑡 = 𝐷𝑖𝑎𝑚𝑒𝑡𝑟𝑜 𝑒𝑥𝑡𝑒𝑟𝑛𝑜 𝑑𝑒 𝑙𝑜𝑠 𝑡𝑢𝑏𝑜𝑠 [𝑖𝑛] =
4

22 de agosto del 2022

Iteraciones

Primera iteración.

Condiciones iniciales de los fluidos. (Todos los datos se encuentran en la página 1).

𝑊 = 1000000 𝑙𝑏/ℎ 𝑤 = 730000 𝑙𝑏/ℎ


𝑇1 = 298 °𝐹 𝑇2 = 197 °𝐹
𝑡1 = 90 °𝐹 𝑡2 = 227.7414 °𝐹
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

𝐵𝑇𝑈 𝐵𝑇𝑈
𝐶𝑝1 = 0,5189 𝐶𝑝2 = 0,5212
𝑙𝑏 ∗ °𝐹 𝑙𝑏 ∗ °𝐹

Temperatura promedio.
Para el fluido caliente.
298 + 197
𝑇𝑝𝑟𝑜𝑚 = = 247.5 °𝐹 𝑬𝒒. 𝟖
2
Para el fluido frío.
90 + 227.7414
𝑇𝑝𝑟𝑜𝑚 = = 158.87 °𝐹 𝑬𝒒. 𝟗
2
UD,inicial.
Para iniciar el proceso iterativo, se debe de suponer un UD,inicial (coeficiente global sucio de
transferencia de calor). Se toma el valor más alto de los reportados por la literatura, tomados
de la tabla 8 del libro de Kern.

Figura 3. Elección UD,inicial de acuerdo a la tabla 8 del libro de Kern.


Teniendo en cuenta las viscosidades de los fluidos se determinó que el fluido caliente es
un orgánico medio (0.5 < µ < 1.0) y el fluido frío es un orgánico pesado (µ > 1.0), así se
escoge el mayor de todo el rango que sería entre 10 y 60.
𝐵𝑇𝑈
𝑈𝐷,𝑖𝑛𝑖𝑐𝑖𝑎𝑙 = 60
𝑓𝑡 2 °𝐹
Distribución de los fluidos
Se determinó que el fluido caliente va a pasar por los tubos y el fluido frio pasará por la
camisa de acuerdo con la matriz de evaluación de criterios para la ubicación de los fluidos.
Parámetros.
𝐵𝑇𝑈
𝑄 = 52405243.1

UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

𝐵𝑇𝑈
𝑘𝑚 = 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑣𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑚𝑎𝑡𝑒𝑟𝑖𝑎𝑙 𝑑𝑒 𝑙𝑜𝑠 𝑡𝑢𝑏𝑜𝑠 [ ] = 9,4242
ℎ ∗ 𝑓𝑡 ∗ °𝐹
𝑃1 = 𝑝𝑟𝑒𝑠𝑖ó𝑛 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑓𝑙𝑢𝑖𝑑𝑜 𝑐𝑎𝑙𝑖𝑒𝑛𝑡𝑒 [𝑃𝑆𝐼] = 280
𝑝1 = 𝑝𝑟𝑒𝑠𝑖ó𝑛 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 𝑓𝑙𝑢𝑖𝑑𝑜 𝑓𝑟í𝑜 [𝑃𝑆𝐼] = 250
Dimensiones del intercambiador.
La resistencia al incrustamiento recomendada interna (Rdi,requerido) y externa (Rdo,requerido) de
los tubos se toman de las recomendaciones de la tabla 12 del Kern. Los fluidos de trabajos
se toman como líquidos industriales orgánicos.

Figura X. Elección Rdo,requerido y Rdi,requerido de acuerdo a la tabla 12 del libro de Kern.

Para calcular el Rd,requerido se parte de la siguiente ecuación.


𝑅𝑑,𝑟𝑒𝑞𝑢𝑒𝑟𝑖𝑑𝑜 = 𝑅𝑑𝑜,𝑟𝑒𝑞𝑢𝑒𝑟𝑖𝑑𝑜 + 𝑅𝑑𝑖,𝑟𝑒𝑞𝑢𝑒𝑟𝑖𝑑𝑜 𝑬𝒒. 𝟏𝟎

𝑅𝑑,𝑟𝑒𝑞𝑢𝑒𝑟𝑖𝑑𝑜 = 0.001 + 0.001 = 0.002

Se debe tener en cuenta que esta tabla está hecha para un tiempo de operación entre 1.5
años, esto quiere decir que se debe hacer limpieza cada vez que se alcanza este tiempo.
Área del flujo liquido por un tubo.
𝜋(𝐼𝐷𝑡 )2
𝑎𝑡′ = = 0.3019 𝑖𝑛2 𝑬𝒒 𝟏𝟏.
4
Área de transferencia de calor por unl tubo.
𝜋(𝑂𝐷𝑡 ) 𝑓𝑡 2
𝑎𝑡′′ = = 0.19635 𝑬𝒒. 𝟏𝟐
12 𝑓𝑡
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Claridad entre los tubos.


3
𝐶 ′ = 𝑃𝑡 − (𝑂𝐷)𝑡 = 𝑖𝑛 𝑬𝒒. 𝟏𝟑
16
Tipo de intercambiador.
El intercambiador es de tubos y coraza donde el número de pasos por los tubos es 4 y el
número de pasos por la camisa es de 2, es decir, n = 4 y m = 2 y el flujo es en
contracorriente.
Separación entre bafles.
La separación entre bafles está en un rango entre el dímetro interno de la camisa y la quinta
parte de este, como se muestra a continuación.
𝐵𝑚á𝑥 = 𝐼𝐷𝑠 𝑬𝒒. 𝟏𝟒
𝐼𝐷𝑠
𝐵𝑚𝑖𝑛 = 𝑬𝒒. 𝟏𝟓
5
Se decidió que la separación entre bafles será:
𝐵 = 𝐼𝐷𝑠 = 39 𝑖𝑛
Factor de corrección por temperatura FT.
Para el cálculo del factor de corrección se debe tener en cuenta que la configuración del
intercambiador es I2-4, con esto se calcula S y R como se muestra en las siguientes
ecuaciones.
𝑇1 − 𝑇2
𝑅= = 0.7333 𝑬𝒒. 𝟏𝟔
𝑡2 − 𝑡1
𝑡2 − 𝑡1
𝑆= = 0.6622 𝑬𝒒. 𝟏𝟕
𝑇1 − 𝑡1
De acuerdo con los resultados obtenidos para S y R se lee en la figura 19 del libro de Kern
para determinar FT.
𝐹𝑇 ≈ 0,94
Área de transferencia de calor basada en el UD seleccionado.
Este valor depende del UD seleccionado anteriormente y se calcula como se muestra a
continuación:
𝑄
𝐴𝑖𝑛𝑐𝑖𝑎𝑙 = 𝑬𝒒. 𝟏𝟖
𝑈𝐷 ∗ 𝐹𝑇 ∗ 𝐿𝑀𝑇𝐷
𝐴𝑖𝑛𝑐𝑖𝑎𝑙 = 10637.92 𝑓𝑡 2
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Número de tubos inicial.


𝐴𝑖𝑛𝑖𝑐𝑖𝑎𝑙
𝑁𝑡,𝑖𝑛𝑖𝑐𝑖𝑎𝑙 = 𝑬𝒒. 𝟏𝟗
𝑎´´𝑡 ∗ 𝐿
𝑁𝑡,𝑖𝑛𝑖𝑐𝑖𝑎𝑙 = 2708.92

Número de tubos corregido.


De acuerdo con el diámetro interno de la camisa, (ID)s, encontrado inicialmente, el arreglo
que se escogió para el intercambiador que es triangular y con la tabla 9 del Kern el número
de tubos corregido es el valor máximo que se puede tomar.

Figura X. Elección de Nt, corregido de acuerdo con la tabla 9 del Kern


𝑁𝑡,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜 = 1258

Coeficiente sucio de transferencia de calor corregido.


Teniendo en cuenta el nuevo número de tubos corregidos.
𝑄
𝑈𝐷,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜 = 𝑬𝒒. 𝟐𝟎
𝐹𝑇 ∗ 𝐿𝑀𝑇𝐷 ∗ 𝑁𝑡,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜 ∗ 𝑎´´𝑡 ∗ 𝐿
𝐵𝑇𝑈
𝑈𝐷,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜 = 129.2
𝑓𝑡 2 °𝐹
Área de transferencia de calor corregida.
Teniendo en cuenta el nuevo número de tubos corregidos.
𝐴𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑎 = 𝑁𝑡,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜 ∗ 𝐿 ∗ 𝑎´´𝑡 𝑬𝒒. 𝟐𝟏

𝐴𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑎 = 4940.15 𝑓𝑡 2
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Área total de flujo por los tubos.


𝑁𝑡,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜 ∗ 𝑎´𝑡
𝑎𝑡 = 𝑬𝒒. 𝟐𝟐
144 ∗ 𝑛
𝑎𝑡 = 0.6593 𝑓𝑡 2
Velocidad másica por los tubos.
𝑊
𝐺𝑡 = 𝑬𝒒. 𝟐𝟑
𝑎𝑡
𝑙𝑏
𝐺𝑡 = 1516591
ℎ 𝑓𝑡 2
Área total de flujo por la camisa.
𝐼𝐷𝑠 ∗ 𝐶 ′ ∗ 𝐵
𝑎𝑠 = 𝑬𝒒. 𝟐𝟒
144 ∗ 𝑃𝑡 ∗ 𝑚
𝑎𝑠 = 1.05625 𝑓𝑡 2
Velocidad másica por la camisa.
𝑤
𝐺𝑠 = 𝑬𝒒. 𝟐𝟓
𝑎𝑠
𝑙𝑏
𝐺𝑠 = 691124.26
ℎ 𝑓𝑡 2
Propiedades de los fluidos.
Las condiciones de entrada y salida del fluido que va por los tubos muestran a continuación.
𝑙𝑏 𝑙𝑏
𝜇𝐸 = 1.028 𝜇𝑆 = 1.692
ℎ ∗ 𝑓𝑡 ℎ ∗ 𝑓𝑡
𝐵𝑇𝑈 𝐵𝑇𝑈
𝐶𝑝𝐸 = 0.541 𝐶𝑝𝑆 = 0.511
𝑙𝑏 ∗ °𝐹 𝑙𝑏 ∗ °𝐹
𝐵𝑇𝑈 𝐵𝑇𝑈
𝑘𝐸 = 0.073 𝑘𝑆 = 0.079
ℎ ∗ 𝑓𝑡 ∗ °𝐹 ℎ ∗ 𝑓𝑡 ∗ °𝐹
𝑙𝑏 𝑙𝑏
𝜌𝐸 = 53.517 𝜌𝑆 = 56.578
𝑓𝑡 3 𝑓𝑡 3
Las condiciones de entrada y salida del fluido que va la camisa muestran a continuación.
𝑙𝑏 𝑙𝑏
𝜇𝐸 = 5.273 𝜇𝑆 = 1.385
ℎ ∗ 𝑓𝑡 ℎ ∗ 𝑓𝑡
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

𝐵𝑇𝑈 𝐵𝑇𝑈
𝐶𝑝𝐸 = 0,496 𝐶𝑝𝑆 = 0.533
𝑙𝑏 ∗ °𝐹 𝑙𝑏 ∗ °𝐹
𝐵𝑇𝑈 𝐵𝑇𝑈
𝑘𝐸 = 0.089 𝑘𝑆 = 0.086
ℎ ∗ 𝑓𝑡 ∗ °𝐹 ℎ ∗ 𝑓𝑡 ∗ °𝐹
𝑙𝑏 𝑙𝑏
𝜌𝐸 = 62.56 𝜌𝑆 = 58.485
𝑓𝑡 3 𝑓𝑡 3
Reynolds de los tubos.
𝐼𝐷𝑡 ∗ 𝐺𝑡
𝑅𝑒𝑡 = 𝑬𝒒. 𝟐𝟔
12 ∗ 𝜇
𝑅𝑒𝑡,𝐸 = 76222.978, 𝑅𝑒𝑡,𝑆 = 46310.414

Reynolds de la camisa.
Teniendo en cuenta que de es el diámetro equivalente de la camisa en pulgadas.
𝑑𝑒 ∗ 𝐺𝑠
𝑅𝑒𝑠 = 𝑬𝒒. 𝟐𝟕
12 ∗ 𝜇
𝑅𝑒𝑠,𝐸 = 5823.67, 𝑅𝑒𝑠,𝑆 = 22171.98

j H.
Este valor se toma de la figura 24 del libro de Kern.
(𝑗𝐻 )𝑡,𝐸 = 205

(𝑗𝐻 )𝑡,𝑆 = 132

(𝑗𝐻 )𝑠,𝐸 = 40

(𝑗𝐻 )𝑠,𝑆 = 88

Coeficientes de transferencia de calor.


Luego, teniendo en cuenta la información previa se calcula los coeficientes de convectivo
de transferencia de calor del fluido dentro de los tubos y del fluido fuera de los tubos, hi y
ho, respectivamente con la siguiente ecuación.
Coeficiente convectivo de transferencia de calor por dentro de los tubos.
1
𝐶 𝜇 3
(𝑗𝐻 )𝑡 ∗ 𝑘 ( 𝑝 )
𝑘
ℎ𝑖 = 𝑬𝒒. 𝟐𝟖
𝐼𝐷𝑡
12
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

𝐵𝑇𝑈 𝐵𝑇𝑈
ℎ𝑖,𝐸 = 569.931 ℎ𝑖,𝑠 = 448.114
ℎ 𝑓𝑡 2 °𝐹 ℎ 𝑓𝑡 2 °𝐹
Coeficiente promedio de transferencia de calor por los tubos referido a Ao.
𝐼𝐷𝑡
ℎ𝑖𝑜 = ℎ𝑖 ∗ 𝑬𝒒. 𝟐𝟗
𝑂𝐷𝑡
𝐵𝑇𝑈 𝐵𝑇𝑈
ℎ𝑖𝑜,𝐸 = 471.14 ℎ𝑖𝑜,𝑠 = 370.441
ℎ 𝑓𝑡 2 °𝐹 ℎ 𝑓𝑡 2 °𝐹
Coeficiente convectivo de transferencia de calor por fuera de los tubos. (µ/µw)0.14
Inicialmente es 1.
1
𝐶 𝜇 3
(𝑗𝐻 )𝑠 ∗ 𝑘 ( 𝑝 ) 𝜇 0.14
𝑘
ℎ𝑜 = ∗( ) 𝑬𝒒. 𝟑𝟎
𝐷𝑒 𝜇𝑤
𝐵𝑇𝑈 𝐵𝑇𝑈
ℎ𝑜,𝐸 = 247.294 ℎ𝑜,𝑠 = 348.746
ℎ 𝑓𝑡 2 °𝐹 ℎ 𝑓𝑡 2 °𝐹
Coeficiente limpio de transferencia de calor.
1
𝑈ℎ = 𝑬𝒒. 𝟑𝟏
𝑂𝐷 𝐷𝑂
1 ∗ 𝑙𝑛 ( )
𝐷𝐼 + 1
+ 12
ℎ𝑖𝑜,𝐸 2𝑘𝑚 ℎ𝑜,𝐸
𝐵𝑇𝑈
𝑈ℎ = 177.901
𝑓𝑡 2 °𝐹
1
𝑈𝑐 = 𝑬𝒒. 𝟑𝟐
𝑂𝐷 𝐷𝑂
1 12 ∗ 𝑙𝑛 ( 𝐷𝐼 ) 1
+ +
ℎ𝑖𝑜,𝑆 2𝑘𝑚 ℎ𝑜,𝑆
𝐵𝑇𝑈
𝑈𝑐 = 135.59
𝑓𝑡 2 °𝐹
Calculo Kc.
𝑈ℎ − 𝑈𝑐
𝐾𝑐 = 𝑬𝒒. 𝟑𝟑
𝑈𝑐
𝐾𝑐 = 0.312
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Calculo r.
∆𝑡𝑐 𝑻𝟐 − 𝒕𝟏
𝑟= = 𝑬𝒒. 𝟑𝟒
∆𝑡ℎ 𝑻𝟏 − 𝒕𝟐
𝑟 = 1.523
Fc.
Se calcula el Fc a partir de la figura 17 del libro de Kern y teniendo en cuenta los datos del
Kc y el r que se obtuvieron anteriormente.
𝐹𝑐 = 0.55
Temperaturas calóricas.
𝑇𝑐ℎ = 𝑇2 + 𝐹𝑐 ∗ (𝑇1 − 𝑇2 ) 𝑬𝒒. 𝟑𝟓
𝑇𝑐ℎ = 248.75 °𝐹
𝑡𝑐𝑐 = 𝑡1 + 𝐹𝑐 ∗ (𝑡2 − 𝑡1 ) 𝑬𝒒. 𝟑𝟔
𝑡𝑐𝑐 = 160.57 °𝐹
Cálculos por los tubos.
Las propiedades de los fluidos se calcularán a partir de las temperaturas calóricas halladas
anteriormente.
𝑙𝑏 𝐵𝑇𝑈
𝜇𝑡 = 1.19 𝐶𝑝,𝑡 = 0.519
ℎ 𝑓𝑡 𝑙𝑏 °𝐹
𝐵𝑇𝑈 𝑙𝑏
𝑘𝑡 = 0.077 𝜌𝑡 = 55.68
ℎ 𝑓𝑡 °𝐹 𝑓𝑡 3
1
𝐶𝑝 𝜇 3
𝑠𝑡 = 0.8909 𝑘( ) = 0.1538
𝑘 𝑡
La velocidad lineal por lo tubos se calcula con la siguiente ecuación.
𝐺𝑡
𝑣𝑡 = 𝑬𝒒. 𝟑𝟕
3600 ∗ 𝜌𝑡
𝑓𝑡
𝑣𝑡 = 7.5656
𝑠
El Reynolds por los tubos se calcula a partir de la ecuación 26.
𝑅𝑒𝑡 = 65845.4182
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

El jH por los tubos se definirá de acuerdo con la figura 24 del libro de Kern.
𝑗𝐻,𝑡 = 180

La viscosidad inicial del fluido por los tubos es igual a la viscosidad hallada con la
temperatura calórica.
𝑙𝑏
𝜇𝑤,𝑡,𝑖𝑛𝑖𝑐𝑖𝑎𝑙 = 𝜇𝑡 = 1.19
ℎ 𝑓𝑡
Inicialmente se asume un valor de 1 para 𝜙𝑡 .
0.14
𝜇𝑡
𝜙𝑡 = ( ) 𝑬𝒒. 𝟑𝟖
𝜇𝑤,𝑡

El coeficiente promedio de transferencia de calor por los tubos se determina por la ecuación
28 y 29 con los nuevos valores obtenidos para el fluido que pasa por los tubos.
𝐵𝑇𝑈 𝐵𝑇𝑈
ℎ𝑖 = 535.699 ℎ𝑖𝑜 = 442.845
ℎ 𝑓𝑡 2 °𝐹 ℎ 𝑓𝑡 2 °𝐹
Caída de presión por los tubos.
Primero se debe tener en cuenta el facto de fricción (Fanning) por los tubos usando la figura
26 del libro de Kern.
𝑓𝑡 2
𝑓𝑡 = 0.00017 2
𝑖𝑛
La caída de presión para los tubos se calcula a partir de las ecuaciones 7.45 y 7.46
corregida por el profesor del libro de Kern, así
𝑓𝑡 ∗ 𝐺𝑡2 ∗ 𝐿 ∗ 𝑛 𝑣𝑡2 𝜌𝑡
∆𝑃𝑐𝑎𝑙𝑐𝑢. = 10
+ 4𝑛 ∗ ( ′ )∗( ) 𝑬𝒒. 𝟑𝟗
52.2 ∗ 𝐼𝐷𝑡 ∗ 𝑠𝑡 ∗ 𝜙𝑡 2∗𝑔 144

Donde,
𝑓𝑡
𝑔′ = 𝑔𝑟𝑎𝑣𝑒𝑑𝑎𝑑 = 32.2
𝑠2
∆𝑃𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 18.52 𝑃𝑆𝐼
La caída de presión permitida es de 5 PSI tanto para los tubos como para la camisa.
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Cálculos por la camisa.


Las propiedades de los fluidos para la camisa se calculan a partir de las temperaturas
calóricas halladas anteriormente.
𝑙𝑏 𝐵𝑇𝑈
𝜇𝑠 = 2.2128 𝐶𝑝,𝑠 = 0.5216
ℎ 𝑓𝑡 𝑙𝑏 °𝐹
𝐵𝑇𝑈 𝑙𝑏
𝑘𝑠 = 0.0864 𝜌𝑠 = 60.1455
ℎ 𝑓𝑡 °𝐹 𝑓𝑡 3
1
𝐶𝑝 𝜇 3
𝑠𝑠 = 0.9623 𝑘( ) = 0.2051
𝑘 𝑠
La velocidad lineal por la camisa se calcula con la siguiente ecuación.
𝐺𝑠
𝑣𝑠 = 𝑬𝒒. 𝟒𝟎
3600 ∗ 𝜌𝑠
𝑓𝑡
𝑣𝑠 = 3.1919
𝑠
La viscosidad inicial del fluido por la camisa es igual a la viscosidad hallada con la
temperatura calórica.
𝑙𝑏
𝜇𝑤,𝑠,𝑖𝑛𝑖𝑐𝑖𝑎𝑙 = 𝜇𝑠 = 2.2128
ℎ 𝑓𝑡
Inicialmente se asume un valor de 1 para 𝜙𝑠 .
0.14
𝜇𝑠
𝜙𝑠 = ( ) 𝑬𝒒. 𝟒𝟏
𝜇𝑤,𝑠

El Reynolds por la camisa se calcula a partir de la ecuación 42.


𝐷𝑒 ∗ 𝐺𝑠
𝑅𝑒𝑠 = 𝑬𝒒. 𝟒𝟐
𝜇𝑠
𝑅𝑒𝑠 = 13877.4341
El jH por los tubos se definirá de acuerdo con la figura 28 del libro de Kern (válida para
hidrocarburos, compuesto orgánicos, agua, soluciones acuosas y gases, y corte bafles del
25%).
𝑗𝐻,𝑠 = 68

El coeficiente promedio de transferencia de calor por la camisa se determina por la ecuación


40.
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

1
𝐶 𝜇 3
(𝑗𝐻 )𝑠 ∗ 𝑘 ( 𝑝 )
𝑘
ℎ𝑜 = 𝑬𝒒. 𝟒𝟑
𝐷𝑒
𝐵𝑇𝑈
ℎ𝑜 = 313.8196
ℎ 𝑓𝑡 2 °𝐹
Verificación de que los valores de ф asumidos fueron correctos.
La temperatura de pared se calcula a partir de la siguiente ecuación
ℎ𝑜
𝑡𝑤 = 𝑇𝑐 − ∗ (𝑇𝑐 ∗ 𝑡𝑐) 𝑬𝒒. 𝟒𝟒
ℎ𝑖𝑜 + ℎ𝑜
𝑡𝑤 = 212.1796 °𝐹
Teniendo la temperatura de la pared se calcula la viscosidad del fluido en los tubos y en la
camisa.
𝑙𝑏 𝑙𝑏
(𝜇𝑤 )𝑡 = 1.4047 (𝜇𝑤 )𝑠 = 1.4455
ℎ 𝑓𝑡 ℎ 𝑓𝑡
Luego, se calcula ф para el tubo y la camisa con las ecuaciones 40 y 43 respectivamente
para revisar si es igual o no a 1.
𝜙𝑡 = 0.98 𝜙𝑠 = 1.061
Caída de presión por la camisa.
Inicialmente, se debe tener en cuenta el factor de fricción (Fanning) por la camisa usando
la figura 29 del libro de Kern.
𝑓𝑡 2
𝑓𝑡 = 0.0021
𝑖𝑛2
La caída de presión para la camisa se calcula a partir de las ecuaciones 7.44 del libro de
Kern, así
𝑓𝑠 ∗ 𝐺𝑠2 ∗ 𝐷𝑠 ∗ (𝑁 + 1)
∆𝑃𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 𝑬𝒒. 𝟒𝟓
5.2210 ∗ 𝐷𝑒 ∗ 𝑠𝑠 ∗ 𝜙𝑠
Donde,
12 ∗ 𝑚 ∗ 𝐿
(𝑁 + 1) = 𝑛ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑐𝑟𝑢𝑐𝑒 𝑝𝑜𝑟 𝑙𝑎 𝑐𝑎𝑚𝑖𝑠𝑎 =
𝐵
∆𝑃𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜 = 16.9360 𝑃𝑆𝐼
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Cálculo del Rd y del exceso de área de transferencia de calor.


El coeficiente limpio de transferencia de calor se calcula con la ecuación 48 y debe ser
mayor que el coeficiente sucio de transferencia de calor UD,corregido
1
𝑈𝐶 = 𝑬𝒒. 𝟒𝟔
𝑂𝐷𝑡
1 1 𝑂𝐷𝑡 𝑙𝑛 ( 𝐼𝐷𝑡 )
+ + 12 ∗
ℎ𝑜 ℎ𝑖𝑜 2 ∗ 𝑘𝑚
𝐵𝑇𝑈
𝑈𝐶 = 164.5856
𝑓𝑡 2 °𝐹
Luego se calcula el área limpia de transferencia de calor, siguiendo la ecuación 14 y
reemplazando UD por el UC que se calculó anteriormente.
𝐴𝑟𝑒𝑎 𝑙𝑖𝑚𝑝𝑖𝑎 = 3878.07 𝑓𝑡 2
Verificar si el Rd,calculado es menor al Rd,requerido.
𝑈𝐶 − 𝑈𝐷,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜
𝑅𝑑,𝑐𝑎𝑙𝑢𝑙𝑎𝑑𝑜 = 𝑬𝒒. 𝟒𝟕
𝑈𝐷,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜 ∗ 𝑈𝐶

𝑅𝑑,𝑐𝑎𝑙𝑢𝑙𝑎𝑑𝑜 = 0.0016

Finalmente, para aceptar el diseño, en cuanto al objetivo de la transferencia de calor, el


exceso de área en el intercambiador debe tener un máximo de 30%, valor que se podría
optimizar.
𝐴𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑎 − 𝐴𝑙𝑖𝑚𝑝𝑖𝑎
𝐸𝑥𝑐𝑒𝑠𝑜 𝑑𝑒 á𝑟𝑒𝑎 = ∗ 100 𝑬𝒒. 𝟒𝟖
𝐴𝑙𝑖𝑚𝑝𝑖𝑎

𝐸𝑥𝑐𝑒𝑠𝑜 𝑑𝑒 á𝑟𝑒𝑎 = 27.39%


¿Qué está mal?
- ΔP calculado > ΔP permitido, tanto para la camisa como para el tubo.
- Gt está por encima del valor permitido
- Vt no se encuentra en el rango requerido
- ∆Pcalculado>∆Ppermitido
- El Rdcalculado < Rdrequerido
- Nt calculado es muy grande.
Se sigue con la iteración 2.
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

25 de agosto del 2022


Hay que tener en cuenta que para que el modelo sea aceptado debe cumplir las siguientes
restricciones.
a) Separación entre los bafles.
7.4 < 𝐵 < 39
b) Velocidad lineal por los tubos.
3 < 𝑣𝑡 < 5
c) Velocidad lineal por la camisa.
0.9 < 𝑣𝑠 < 3.2
d) Caída de presión por los tubos y la camisa.
∆𝑃 ≤ 5
e) Coeficiente limpio de transferencia de calor.
𝑈𝐶 > 𝑈𝐷,𝑐𝑜𝑟𝑟𝑒𝑔𝑖𝑑𝑜

f) Resistencia de incrustamiento combinada.


𝑅𝑑,𝑐𝑎𝑙𝑢𝑙𝑎𝑑𝑜 > 𝑅𝑑,𝑟𝑒𝑞𝑢𝑒𝑟𝑖𝑑𝑜 = 0.002

g) Exceso de área.
𝐸𝑥𝑐𝑒𝑠𝑜 𝑑𝑒 𝑎𝑟𝑒𝑎 ≤ 30%
Segunda iteración.
Cambios realizados: disminuir el caudal del flujo, aumentar el valor inicial de UD y
aumentar B.
Restricciones.
a) Si se cumple
b) Si se cumple
c) Si se cumple
d) No cumple
e) Si se cumple
f) Si se cumple
g) No cumple
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

¿Qué está mal?


- El exceso de área es mayor al 30%
- ∆Pcalculado>∆Ppermitido para los tubos
26 de agosto del 2022
Tercera iteración.
Cambios realizados: aumentar el UD inicial y disminuir el IDs
Restricciones.
a) Si se cumple
b) Si se cumple
c) Si se cumple
d) No cumple
e) Si se cumple
f) Si se cumple
g) Si se cumple
¿Qué está mal?
- ∆Pcalculado>∆Ppermitido para los tubos
29 de agosto del 2022
Cuarta iteración.
Cambios realizados: aumentar L, aumentar los diámetros y aumentar B.
Restricciones.
a) Si se cumple
b) Si se cumple
c) Si se cumple
d) Si se cumple
e) Si se cumple
f) Si se cumple
g) Si se cumple
¿Qué está mal?
- Todo está bien.
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Diámetro interno de las boquillas


Para calcular los diámetros internos de las boquillas en la entrada y la salida de los tubos y
la camisa, se utiliza como guía la ecuación 9-76 y 9-77 del libro Plant Design and Economics
for Chemical Engineers que se puede aplicar para flujos turbulentos.

Cuando el diámetro es ≥ 0.0254 𝑚


0.363 ∗ 𝑚0.45 ∗ 𝜇𝑐0.025
𝐷𝑖,𝑜𝑝𝑡𝑖𝑚𝑜 =
𝜌0.32

Cuando el diámetro es < 0.0254 𝑚


𝐷𝑖,𝑜𝑝𝑡𝑖𝑚𝑜 = 0.49 ∗ 𝑞𝑓0.49 ∗ 𝜌0.14 ∗ 𝜇𝑐0.027
Donde,

𝑘𝑔
𝑚 = 𝑓𝑙𝑢𝑗𝑜 𝑚𝑎𝑠𝑖𝑐𝑜 [ ]
𝑠
𝑘𝑔
𝜌 = 𝑑𝑒𝑛𝑠𝑖𝑑𝑎𝑑 [ 3 ]
𝑚
𝜇𝑐 = 𝑣𝑖𝑠𝑐𝑜𝑠𝑖𝑑𝑎𝑑 [𝑃𝑎 ∗ 𝑠]
𝑚3
𝑞𝑓 = 𝑓𝑙𝑢𝑗𝑜 𝑣𝑜𝑙𝑢𝑚𝑒𝑡𝑟𝑖𝑐𝑜 [ ]
𝑠
En la siguiente tabla, se muestran las características de los flujos de entrada y salida

Tabla 1. Propiedades de los flujos de entrada y salida de los tubos y camisa.


Propiedades Unidades Entrada tubos Salida tubos Entrada camisa Salida camisa
m [Kg/s] 62,99 62,99 45,98 45,98

ρ [Kg/m3] 857,262 906,295 1.002,117 936,842

qf [m3/s] 0,073 0,070 0,046 0,049

μc [Pa. S] 0,0004250 0,0006994 0,0021797 0,0005725

Como el diámetro interno de los tubos y de la camisa es mayor a 0.0254 m, se utiliza la ecuación ..
En la tabla, se encuentran registrados los diámetros calculados y los comerciales que se leyeron de
la tabla 11 del libro de Kern.
UNIVERSIDAD DE ANTIOQUIA
FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA QUÍMICA
CURSO OPERACIONES DE TRANSFERENCIA DE CALOR
DOCENTES: MARIANA PEÑUELA, LUIS ALBERTO RÍOS

Tabla 2. Diámetros calculados y comerciales internos de las boquillas de entrada y salida


de las corrientes.

IDE.calculado [in] IDE.comercial [in] IDS.calculado [in] IDS.comercial [in]

Tubos 11,205 12,09 11,138 12,09

Camisa 7,524 7,625 7,435 7,625

Selección de tipo de brida


Brida ciega (Blind), están destinadas a cerrar extremos de tubería, válvulas o aberturas de
recipientes, sometidos a variadas presiones de trabajo. Desde el punto de vista técnico, este
tipo de bridas es el que soporta condiciones de trabajo más severas (particularmente las de
mayores dimensiones), ya que, al esfuerzo provocado por la tracción de los bulones, se la
adiciona el producido por la presión existente en la tubería. En los terminales, donde la
temperatura sea un factor de trabajo o actúen esfuerzos variantes o cíclicos, es aconsejable
efectuar los cierres mediante el acople de bridas con cuello y ciegas(ASME B16.5).

Referencias.

Kern, D. Q. (1983). Process heat transfer (International student ed.).

L. Ríos, M. Peñuela (2018, abril) “Ingeniería Química, Universidad de Antioquia”.

M. Peters. et al. “Plant Design and Economics For Chemical Engineers”

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy