(Español) Chemical Process Design and Simulation of
(Español) Chemical Process Design and Simulation of
Juma Haydary
Departamento de Ingeniería Química y Bioquímica
Universidad de Tecnología de Eslovaquia
Bratislava, Eslovaquia
Machine Translated by Google
Una publicación conjunta del Instituto Americano de Ingenieros Químicos y John Wiley & Sons, Inc.
Reservados todos los derechos. Ninguna parte de esta publicación puede reproducirse, almacenarse en un sistema de recuperación o transmitirse de ninguna forma ni por ningún medio, ya
sea electrónico, mecánico, fotocopiado, grabación u otros, excepto según lo permita la ley. Los consejos sobre cómo obtener permiso para reutilizar material de este título están disponibles en
http://www.wiley.com/go/permissions.
El derecho de Juma Haydary a ser identificado como el autor de esta obra se ha hecho valer de conformidad con la ley.
Oficina registrada
John Wiley & Sons, Inc., 111 River Street, Hoboken, NJ 07030, EE. UU.
Oficina editorial
111 River Street, Hoboken, Nueva Jersey 07030, EE. UU.
Para obtener detalles sobre nuestras oficinas editoriales globales, servicios al cliente y más información sobre los productos Wiley, visítenos en www.wiley.com.
Wiley también publica sus libros en una variedad de formatos electrónicos y bajo demanda. Parte del contenido que aparece en las versiones impresas estándar de este libro puede no estar
disponible en otros formatos.
10 9 8 7 6 5 4 3 2 1
Machine Translated by Google
viii
Contenido
viii Contenido
3 Intercambiadores de calor
3.1 63 Modelos de calentadores y
3.2 enfriadores 63 Modelos de intercambiadores
3.3 de calor simples 66 Diseño simple y clasificación de
3.4 intercambiadores de calor 69 Diseño detallado y simulación de
intercambiadores de calor 72 3.4.1 Clasificación
dinámica de HYSYS 74 3.4.2 Diseño riguroso de intercambiadores de calor de carcasa y
tubos con EDR 76 3.5 Selección y costeo de intercambiadores de
calor 77 Referencias 82
Contenido ix
Contenido x
Contenido xi
15.2 Simulación de remoción de gas amargo por solución acuosa de aminas 355
15.3 Modelado basado en tasa de absorbentes con electrolitos 361
Referencias 365
dieciséis
Simulación de procesos de producción de polímeros 367
16.1 Descripción general del modelado del proceso de polimerización en Aspen Plus 367
16.2 Caracterización de componentes 368 Método
16.3 de propiedad 369
16.4 Cinética de reacción 370
16.5 Diagrama de flujo del proceso 375
16.6 Resultados 379
Referencias 383
Índice 387
Machine Translated by Google
XIII
Lista de tablas
1.1 Lista de los simuladores de procesos más conocidos 10 6.7 Cálculo de fracciones molares a partir de la
concentración másica 145
2.1 Algunas propiedades de los componentes del proceso de
acetato de etilo 18 2.2 6.8 Resultados del proceso de absorcióndesorción
Bancos de datos de propiedades físicas de ASPEN 26 Submodelos de gases ligeros 147 Resultados
disponibles en Aspen Plus 27 2.3 2.4 Algunas de la especificación de diseño 150 6.9 6.10
ecuaciones cúbicas de estado en Aspen Physical Property Costo del equipo 153 Datos de la curva de
System y Aspen HYSYS 28 Ecuación de modelos 7.1 secado 156 Distribución del
7.2 tamaño de partículas 156 Solubilidad de
de estado 29 Azeótropos
2.5 de acetato de etilo–etanol – 7.3 CuSO4 a diferentes temperaturas 161
2.6 Resultados del análisis
mezcla de agua 45 7.4 de sensibilidad: Efecto de temperatura 162
2.7 Puntos singulares de la mezcla Balance de materia de
acetato de etiloetanolagua 45 7.5 filtración 163 7.6 Distribución del tamaño de
3.1 Coeficiente de transferencia de calor del lado del tubo partículas de ceniza 164 8.1 Algunas propiedades de
correlaciones 73 los componentes puros del proceso de acetato de etilo 182 8.2
3.2 Correlaciones del coeficiente de transferencia Algunas propiedades de los componentes
de calor del lado de puros de
la carcasa 73 3.3 Geometría del intercambiador de calor el proceso de estireno 184
utilizado en el Ejemplo 8.3 Parámetros cinéticos del proceso de acetato de
3.4 73 4.1 Datos de la curva de rendimiento de la bomba etilo de diferentes fuentes 184
4.2 86 Composición del gas natural utilizado en el 8.4 Energía de activación de la deshidrogenación
Ejemplo 4.2 89 Correlaciones catalítica de etilbenceno 185
4.3 utilizadas para el cálculo de la caída de presión de la 8.5 Parámetros de las columnas de destilación 190
tubería en Aspen HYSYS 92 Rango de 8.6 Flujo de resultados para la columna de destilación
4.4 funcionamiento de algunos tipos de bombas 97 Información reactiva 191
4.5 de diseño mapeada de la simulación 98 8.7 Flujo de resultados para el separador de fase líquido
Resultados del líquido 192
4.6 cálculo de costos del compresor 99 8.8 Transmita los resultados de la columna de
5.1 Composición de reactivos y productos del proceso purificación de acetato de etilo 192
de acetato de etilo 114 Resultados del 8.9 Flujo de resultados para la columna de
6.1 modelo de separador HYSYS 127 Resultados del modelo recuperación de ácido acético 192
6.2 Aspen Plus FLASH3 128 Balance de materia 8.10 Resultados de flujo para la columna de destilación de
y energía del fase acuosa 193
6.3 columna 135 Condiciones y composiciones 8.11 Resultados de la simulación de la parte de reacción de
de las corrientes el proceso de estireno 196
6.4 de material en la destilación extractiva 139 Resultados 8.12 Resultados de la simulación de la parte de separación del
de la destilación azeotrópica 141 Resultados de la
proceso de estireno 199 Flujos másicos
6.5 simulación de destilación reactiva 144 9.1 y composiciones de las corrientes de
6.6 alimentación 205
Machine Translated by Google
Parámetros de las columnas 206 9.2 9.3 12.18 Costos operativos anuales totales 273 12.19
Geometría de los intercambiadores de calor 206 Balance Resumen del análisis de rentabilidad del proceso de
de materia de las columnas 210 9.4 9.5 metanol 273
Flujos de energía del proceso 211 9.6 12.20 Datos del diagrama de flujo de efectivo para el proceso
Comparación de la composición del producto medida y de metanol 274
simulada 212 9.7 Comparación de 13.1 Datos de destilación de TBP 290 13.2
escenarios con y sin columna D215 213 9.8 Composición de las luces de los extremos 291 13.3
Comparación de flujos de energía Marcadores laterales y bombeos
214 10.1 Material resultados de balance para MIXER1 y especificaciones 294 13.4
Resultados de la corriente para la columna preflash 297 13.5
MEZCLADOR2 220 Resultados de la corriente de la columna atmosférica 298 13.6
10.2 Resultados del balance de materia para el reactivo Destilación medida experimentalmente
columna de destilación (RD) 220 10.3 curvas 301
Resultados del balance de materia para el separador de fase 13.7 Modelos de reactores de refinería de Aspen HYSYS 308
LL (DEC) 221 13.8
10.4 Resultados del balance de materia para el etilo Parámetros cinéticos de hidrocraqueo VR 310 13.9 Curvas de
columna de purificación de acetato (C1) 221 10.5 destilación TBP y grados API de productos 311 13.10
Tolerancias internas de Aspen HYSYS 223 10.6 Resultados del Propiedades de
proceso de estireno después del reciclaje de etilbenceno alimentación FCC 315 14.1 Atributos de los
224 10.7 Resultados del cálculo del componentes en Aspen Plus 322 14.2 Tipos de clases de
requerimiento de vapor 232 11.1 Lista de componentes para el proceso corrientes en Aspen Plus 322 14.3 Características de los biomasa
de deshidrogenación de nheptano 240 11.2 322 14.4 Parámetros de la curva de secado de la biomasa
Requerimiento de calor de proceso utilizada 323 14.5 Resultados de la corriente para la biomasa
calculado 242 11.3 Corrientes de material flujo y condiciones 243
11.4 Información sobre las corrientes de calor 244 11.5 Condiciones proceso de combustión 330 14.6
de la corriente del proceso de estireno después de Constantes de equilibrio de las principales reacciones de
la recuperación de energía 247 11.6 Resultados de la corriente gasificación del carbón 334 14.7
de la combustión de residuos en Una composición elemental típica de RDF y sus componentes
334 14.8 Contenido de alquitrán
el proceso de estireno 250 11.7 de gas productor 334 14.9 Lista de componentes
Pureza de los productos, servicio del rehervidor y potencia para la gasificación de RDF 335 14.10 Resultados de flujo del
del compresor 253 12.1 Balance proceso de gasificación de RDF 338 14.11 Resumen de
de materia de la compresión de gas de síntesis 264 12.2 Balance de RDF resultados
materia de la sección de reacción del proceso de gasificación en condiciones óptimas 341
del proceso de metanol 266
12.3 Balance de materia de la destilación 14.12 Atributos de biomasa y char 341 14.13 Lista de
sección del proceso de metanol 267 12.4 componentes para la pirólisis de biomasa 342 14.14
Parámetros y costos de los intercambiadores de calor 268 12.5 Pseudocomponentes definidos en el proceso de
Parámetros y costos de las bombas 268 12.6 Parámetros pirólisis de biomasa 342 14.15 produce
y costos de los compresores 268 12.7 Parámetros y costos de las 342 14.16 Resultados del proceso
columnas de destilación 269 de pirólisis de biomasa 344 15.1 Parámetros de pares de
electrolíteos en HCL Removal de la transmisión de la
12.8 Costos de los separadores vaporlíquido (VL) 269 12.9 corriente de víndico de la Vinyl. 353 15.3 Flujo de resultados
Parámetros y costos del reactor 269 12.10 Costos adquiridos y para la eliminación de HCl en el
costos totales del equipo instalado 269 12.11 Costos
de capital indirecto 269 12.12 proceso de cloruro de vinilo 354 15.4
Inversión de capital total del proceso de Lista de componentes para la limpieza de gas de síntesis
metanol 270 12.13 Costos operativos 271 12.14 Costos de con amina 357
servicios públicos 15.5 Esquema de reacción del proceso de amina para la limpieza
del proceso de metanol 272 12,15 de gas de síntesis 358 15.6
Mano de obra operativa 272 12,16 Otros costos de fabricación Datos de especificación de la corriente de entrada 358
272 12,17 Gastos generales 273 15.7 Resultados del depurador de HCl (ABS1) 359 15.8
Resultados del absorbedor de amina (ABS2) ) 359 15,9
Resultados de la columna de regeneración (DC1) 360
Machine Translated by Google
Lista de Tablas xv
15.10 HETP calculado 363 15.11 16.1 Resumen de los métodos de propiedad de
Flujo de resultados del cálculo basado en tasa polímeros de Aspen 372
de un depurador de HCl 364 16.2 Resultados de la corriente de polimerización de estireno 379
15.12 Efecto del tipo y tamaño del empaque en la 16.3 Resultados de conversión y propiedades del
eliminación de HCl 365 polímero 380
Machine Translated by Google
xvii
Lista de Figuras
1.1 Niveles de jerarquía para el diseño de ingeniería 2.20 Comparación del diagrama x,y isobárico de la mezcla
química 4 binaria nheptano/tolueno 30 2.21 Selección de
1.2 Pasos principales del diseño de una planta química los métodos Asistencia 31 2.22 Selección del tipo de sistema
4 PFD simplificado para el diseño del componente 32 2.23 Métodos Recomendación de asistencia 32 2.24
proceso de producción de acetato de vinilo Especificación del método seleccionado en
1.3 7 1.4 PFD más complejo para la producción de
acetato de vinilo la hoja de Propiedades 33
1.5 8 Estructura de un simulador de proceso 2.25 Descripción del método NRTL 33 2.26 Descripción
9 1.6 Abrir una nueva simulación de Aspen Plus 11 de la ecuación de estado de HydenO'Connell 34
Selección del tipo de simulación 12 1.7 1.8
Entorno de propiedades de Aspen Plus 12 1.9 Abrir una 2.27 Parámetros de interacción binaria de la ecuación
nueva simulación de Aspen HYSYS 13 1.10 Entorno de NRTL 34 2.28 Parámetros
propiedades de Aspen HYSYS 13 2.1 Lista de componentes de interacción binaria de la ecuación de estado HOC
del proceso de acetato de etilo 16 2.2 Motor de 35 2.29 Recomendación de los
búsqueda de métodos
componentes de Aspen Plus 16 2.3 Renombrar el Asistencia para sistemas de hidrocarburos 35
componente 17 Parámetros escalares de 2.30 Selección de un paquete de propiedad en Aspen
de componentes puros 17 2.4 Creación de una lista HYSYS 36
2.5 componentes en Aspen HYSYS 18 2.31 Parámetros de interacción binaria del
Ecuación de estado de PengRobinson 36
2.6 Página de propiedades del componente en Aspen 2.32 Descripción de la ecuación de estado de Peng
HYSYS 19 Robinson en Aspen HYSYS Help 37 2.33 Hoja de
2,7 Creación de un hipocomponente 19 Creación análisis de propiedades de componentes puros 37 2.34 Dependencia
2,8 de un conjunto de hipocomponentes 20 Estructura de la temperatura del calor molar de vaporización de los
2,9 química del dibenzo (ah) antraceno 20 2.10 componentes del proceso de acetato de etilo 38
Página de estructura 2.35 Visualización de
molecular 20 2.11 Herramienta de dibujo de los resultados del análisis de propiedades de componentes
estructura molecular 21 2.12 Estructura molecular y puros 38 2.36 Dependencia
cálculo de enlace de la temperatura de viscosidad dinámica para los
página 21 componentes del proceso de acetato de
2.13 Estimación de parámetros de componentes puros etilo 38 2.37 Entrada
página 22 de datos de análisis binario página 39 2.38 Diagrama
2.14 Resultados de parámetros de componentes Txy de acetato de etiloetanol
estimación 22 mezcla 40
2.15 Selección del tipo de componente 23 2.16 2.39 Visualización de datos completos de equilibrio isobárico
Introducción de parámetros conocidos de un 40
pseudocomponente 23 2.40 Mostrar otros gráficos usando la barra de
2.17 Revisión de parámetros desconocidos 24 2.18 herramientas Gráfico 41
Parámetros desconocidos calculados 24 2.19 2.41 Uso del motor NIST ThermoData 42 2.42 Selección de
Comparación del diagrama isobárico t–x,y de la mezcla datos isobáricos VLE en el
binaria nheptano/tolueno 30 Motor NIST ThermoData 42
Machine Translated by Google
2.43 Comparación de datos VLE para el sistema 2.72 Diagrama de flujo de hidrodesalquilación de tolueno 58
binario acetato de etiloetanol 43 2.44 Diagrama de flujo del modelo de
Comparación de datos VLE de etanolácido acético 43 2.45 calentador 64 3.1 Especificación de la corriente
de
Comparación del entrada 65 3.2 Especificación del bloque
de datos VLE para el calentador por 3.3 temperatura y presión 65
sistema binario etanolácido acético 43 2.46 3.4 Visualización de los resultados del bloque del
Comparación de datos VLE para el sistema binario agua 3.5 calentador 66 Especificación del bloque del
ácido acético 43 2.47 calentador ajustando la presión y VF 67
Selección de mapas ternarios de síntesis de destilación 44 3.6 Selección del modelo de intercambiador de calor en
2.48 Página Aspen HYSYS 68
de entrada de parámetros de síntesis de destilación 45 3.7 Conexión de corrientes en el modelo de
2.49 intercambiador de calor 68
Diagramas ternarios del sistema acetato de etilo 3.8 Introducción de parámetros de entrada
etanolagua 46 2.50 Curva de utilizando la página de la
residuos mapa del sistema acetato de etiloetanol 3.9 hoja de trabajo 69 Introducción de la composición
agua 46 de la corriente en Aspen HYSYS 70
2.51 Vía de separación del acetato de etilo 46 2.52 Esquema 3.10 Especificación del intercambiador de calor en Aspen
de proceso propuesto para el Plus 71
separación del acetato de etilo de la mezcla 3.11 Especificación de la transferencia de calor total
ternaria 47 2.53 Página método del coeficiente 71
de análisis de envolvente PT 47 2.54 Resultados 3.12 Resultados del cálculo nominal simple 72 3.13 Detalles
de análisis de envolvente PT para todo el rango de del intercambiador de calor existente 74 3.14 Selección del
temperaturas y presiones 48 2.55 modelo de intercambiador de calor 75 3.15 División de la
Resultados de carcasa en zonas 75 3.16 Dimensionamiento
análisis de envolvente PT para un rango específico de detallado del intercambiador de calor 76 3.17 Selección del
temperaturas y presiones 48 2.56 Adición de un modelo riguroso de carcasa y tubos en Aspen HYSYS
conjunto de reacción en Aspen HYSYS 49 2.57 Adición de 77 3.18 Selección de asignación de
reacciones a un conjunto de reacciones y el conjunto de fluido caliente 77 3.19 Geometría de transferencia
reacciones al Fluid Package 50 definida en HYSYS 78 3.20 Resultados del modelo riguroso de
2.58 Definición de la estequiometría y conversión carcasa y tubos 78 3.21 Perfil de temperatura del
de reacciones 50 2.59 Definición de intercambiador
la estequiometría y método de cálculo de Ke 51 2.60 de calor calculado por EDR 79 3.22 Activación del
Creación de una nueva analizador económico en
química en Aspen Plus 51 2.61 Selección de la química Aspen Plus 80 3.23 Mapeo del modelo de operación unitaria
especificando el método y el tipo de reacción 52 2.62 81 3.24
Definición de la estequiometría de Dimensionamiento del equipo 81 3.25 Evaluación económica y
reacciones en Aspen Plus 52 costo del equipo 82 4.1 Ejemplo de
curvas de rendimiento de una bomba centrífuga 86 Diagrama de flujo
2.63 Estequiometría de disociación de electrolitos 53 2.64 Química para un ejemplo de bomba 87 Especificación
de neutralización definida en Aspen Plus 53 2.65 PFD para de la bomba en Aspen Plus si
proceso de acetato 4.2 se conoce la presión de descarga 87
de etilo CSTR en fase líquida 54 2.66 Selección de modelos 4.3 Resultados del bloqueo de la bomba si se conoce
de bloque la presión de descarga 87 4.5 Curva de
predefinidos en Aspen Plus 55 2.67 Selección múltiple de 4.4 rendimiento de la bomba especificación 88
un modelo de bloque 4.6 4.7 4.8
en Aspen Plus 55 2.68 Selección de material corrientes
56 2.69 PFD del Resultados de la simulación de bombas en Aspen Plus 88
proceso de acetato de etilo en Aspen Plus 56 2.70 Diagrama de flujo del compresor 90
Un diagrama de flujo simple del proceso de hidrodesalquilación Especificación del flujo de entrada del
de tolueno compresor 90
57 2.71 Selección y definición de conexiones 4.9 Especificación de los parámetros del compresor 91
4.10 Resultados del cálculo del compresor 91 4.11 Sistema
de tuberías para el cálculo de la caída de presión 94 4.12 Definición
de modelos predefinidos 57 de los segmentos de tubería 95
Machine Translated by Google
4.13 Especificación de las condiciones de transferencia de calor 5.27 Parámetros cinéticos de heterogéneos
del sistema de tuberías 95 reacción catalítica 117
4.14 Condiciones de la corriente de salida del 5.28 Diagrama de flujo del modelo PFR 117
sistema de tuberías 5.29 Introducción de datos del catalizador
96 4.15 Perfil de presión a lo largo de la tubería 118 5.30 Tamaño de PFR en Aspen HYSYS
96 4.16 Clasificación de compresores y 118 5.31 Selección de un bloque de hoja de
sopladores en APEA 97 cálculo 118 5.32 Definición de variables de
4.17 Opciones de mapeo del analizador económico 98 importación en el bloque de
4.18 Mapeo de equipos por APEA 98 hoja de cálculo 119 5.33 Cálculo de la selectividad al
5.1 Diagrama de flujo del modelo de reactor de estireno en el bloque de
conversión en hoja de cálculo 120 5.34 Temperatura y perfiles de
5.2 Aspen Plus 103 Especificación de presión en PFR para producción
Rstoic ingresando temperatura y de estireno 120 5.35 Perfil de composición en PFR
5.3 presión 103 Especificación de para producción de
estequiometría y conversión estireno 121 5.36 Conversión de etilbenceno
5.4 en Rstoic 104 Creación de un bloque de y selectividad a estireno 121
calculadora 5.37 Lista de agitadores y tanques agitados
5.5 en Aspen Plus 104 Definición del bloque de disponibles en APEA 122
calculadora 5.38 Selección del tipo de material 123 5.39
5.6 en Aspen Plus 105 Fórmula y secuencia de Relación entre el volumen del reactor, conversión
cálculo en el bloque de y costo del equipo 123 5.40 Costo de los
5.7 calculadora 105 Configuración de las opciones de servicios públicos en función del reactor
5.8 informe en Aspen Plus 106 Composición volumen y conversión 124
de los productos de 6.1 Procesos de separación en una planta
del hidratación del óxido de etileno 106 Configuración química 126
cálculo del calor de reacción 107 5.9 5.10 Resultados del calor de reacción
6.2 Esquema de una separación continua de
cálculo 107 5.11 líquidovapor de una etapa 126
Adición de un conjunto de reacciones al modelo de Diagrama de flujo de destilación flash de una etapa
reactor en Aspen HYSYS 127 6.3 6.4 Página de parámetros de
109 5.12 Inicio de un estudio de caso en Aspen HYSYS interacción binaria UNIQUAC en Aspen
109 5.13 Selección de variables para el caso de estudio 6.5 Plus 128 Bloque de operación de unidad flash
110 5.14 Especificación del rango de variable trifásica en Aspen
independiente para el caso de 6.6 Plus 128 Diagrama de flujo del modelo de
estudio 110 5.15 Composición de los productos de reacción frente a acceso directo
temperatura de reacción 111 6.7 en Aspen HYSYS 129 Especificación de
5.16 Diagrama de flujo CSTR en Aspen Plus 112 parámetros para el cálculo de
5.17 Especificación de RCSTR en Aspen Plus 112 5.18 6.8 destilación de acceso directo 130 Resultados del
Definición de reacciones para modelos de reactores 6.9 modelo de destilación de acceso
cinéticos 113 directo 130 Esquema general de
5.19 Selección de la clase de reacción y separación multietapa y multicomponente 131 6.10
estequiometría 113 Página de conexión de
5.20 Introducción de parámetros cinéticos 113 5.21 6.11 la columna de destilación de Aspen HYSYS
132
Introducción de la constante de equilibrio para inversión
reacción 114 Pestaña final del experto de entrada de la columna de destilación 133 6.12 Ejecu
5.22 Adición de reacción al modelo CSTR en Aspen columna 133
Plus 114 5.23 Inicio 6.13 Adición de nueva especificación de columna 134
de un análisis de sensibilidad en Aspen Plus 115 6.14 Definición de especificación de columna por
5.24 la fracción molar componente 134
Definición de parámetros para el análisis de 6.15 Cálculo de parámetros de columna con
sensibilidad 115 nuevas especificaciones
5.25 Definición de un parámetro local en el análisis de 135 6.16 Número de etapas teóricas versus
sensibilidad 116 relación de reflujo 136
5.26 Resultados del análisis de sensibilidad, conversión 6.17 Diagrama de flujo de la destilación extractiva sin
versus volumen del reactor 116
reciclado de disolventes 137
Machine Translated by Google
xx Lista de figuras
6.18 Especificación del bloque de operación de la 7.11 Definición de la estequiometría de cristalización 161 7.12
unidad Radfrac 138 Resultados de la simulación del proceso de
6.19 Especificación de corrientes de alimentación, cristalización 162
etapas y presión en la columna 138 7.13 Diagrama de flujo de un proceso de cristalización con
6.20 Pureza de nheptano frente a la etapa de alimentación filtración 163 7.14
de NMP 139 Especificación del bloque de operación de la unidad de filtro 163
6.21 Pureza del producto frente al requisito específico del 7.15 Diagrama de flujo del ciclón 164 7.16
solvente 139 6.22 Definición de una subcorriente sólida NC 165 7.17
Diagramas txy binarios isobáricos 140 6.23 Mapa Especificación de la operación de una unidad de ciclón
ternario de la acetona, benceno , y sistema de bloque 165
ciclohexano 140 6.24 Esquema de 7.18 Resultados de la simulación de ciclones 166
una columna de destilación azeotrópica en Aspen 7.19 Curva de eficiencia de separación de un ciclón 166 7.20
HYSYS 141 6.25 Diagrama de flujo de Esquema de tipos de secadores disponibles en APEA 167 Proceso
destilación reactiva del proceso de acetato de etilo 142 8.1 continuo de producción de acetato de etilo
6.26 Configuración de columna 177 Deshidrogenación
de destilación reactiva 143 8.2 adiabática de etilbenceno a estireno
178 Deshidrogenación isotérmica
6.27 Etapas reactivas y especificación de retención 143 6.28 8.3 de etilbenceno a estireno 179 Enfoque
Definición de una reacción química de tipo equilibrio 144 estándar en estireno 180 8.5 Enfoque
8.4 de Monsanto en la separación de estireno 181 8.6
6.29 Perfil de temperatura de la columna de destilación Presión de vapor de los componentes del proceso de estireno
reactiva 145 versus temperatura 184 Datos de equilibrio de fase
6.30 Perfil de composición de la columna de destilación líquidovapor de etilbencenoestireno a 5 kPa
reactiva 145 8.7 186 Datos de equilibrio líquidovapor
6.31 Definición de la eficiencia de la columna 146 isobárico de etilbencenoestireno 186
6.32 Diagrama de flujo del absorbentedesorbedor 146 8.8 Diagrama xy isobárico de etilbenceno
6.33 Especificación del desorbedor por la tasa de – sistema binario de estireno 186 8.10 Efecto
producto de cabeza 147 6.34 8.9 de la presión sobre la volatilidad relativa del
Diagrama ternario del nheptano, benceno y DMSO etilbenceno a estireno 187 8.11
sistema 148 6.35 Diagrama de flujo del Datos de equilibrio de vapor líquido isobárico para el
proceso de extracción multietapa 149 6.36 Selección sistema binario toluenoetilbenceno
de componentes 187 8.12 Diagrama xy isobárico del sistema binario
clave para el líquido toluenoetilbenceno 187 8.13 Diagrama de
conexión fases y corriente extractora 149 flujo preliminar
6.37 Definición de la especificación de diseño 150 6.38 del proceso de acetato de etilo 189
Definición de la variable a especificar y manipular 150 6.39
Componentes del
modelo de operación de la unidad RadFrac 151
6.40 Mapeo de una columna 8.14 Opciones para continuar con la especificación
de destilación 152 6.41 Selección del tipo de plato en una de flujo 189 8.15
columna de plato Definición de una especificación de diseño en Aspen Plus
y tipo de empaque en una columna empaquetada 153 190
7.1 Ejemplo de curva de secado 156 7.2 7.3 7.4 8.16 Definición de un parámetro variable en Aspen Plus
7.5 Ficha de caracterización de sólidos 191
157 Diagrama de flujo del proceso de 8.17 Diagrama de flujo de la primera etapa del proceso
secado 157 Definición de corrientes de estireno 194
sólidas 158 Especificación de un bloque de 8.18 Sustitución de flujos definidos previamente por nuevos
operación de flujos 194 8.19 Diagrama
7.6 la unidad de secado 158 Resultados de la simulación de flujo de la parte de reacción del proceso de estireno 195
7.7 de secado convectivo 159 Gráficos disponibles 8.20 Diagrama de flujo de
para el modelo de secador la parte de separación de
7.8 convectivo en Aspen Plus 159 Perfil de temperatura de un el proceso de estireno 198
7.9 secador convectivo 160 Perfil de humedad de un 9.1 Diagrama de flujo del proceso derivado de
convectivo secador 160 7.10 Diagrama de flujo del proceso de cristalización 161 un esquema tecnológico 204
Machine Translated by Google
9.2 Esquema de simulación obtenido de un PFD más 10.22 Definición de una hoja de cálculo en Aspen HYSYS 236 10.23
complejo 205 Constante Definición del operador de ajuste 237 11.1 Una
9.3 de equilibrio calculada y experimental de propano a simulación simple y directa de la deshidrogenación de n
pentano versus presión a 65.35 ◦C 207 Datos de heptano 240 11.2 Definición
equilibrio isobárico para nbutano/ de un bloque de cálculo en Aspen Plus 241 11.3 Cálculo del
9.4 npentano a 2,068 kPa 207 Datos requerimiento de combustible mediante la interconexión
de equilibrio isotérmico para el sistema C3/C4 a de bloques 241
9.5 53 ◦C 207 Cambio entre los modos solver 11.4 Resultados del cálculo del requerimiento de
activo y en espera en Aspen HYSYS metano sin recuperación de energía 242
9.6 208 9.7 Diagrama de flujo HYSYS para los gases 11.5 Uno de los arreglos posibles para la recuperación
ligeros de calor 243 11.6
Combustión de gas natural y uso del calor de los gases de
proceso de separación 209 combustión en el proceso de estireno 245
9.8 Parámetros de corrientes de pentano reciclado 11.7 Simulación de recuperación de energía del
y eliminado 210 proceso de producción de estireno 246
9.9 Resultados obtenidos para los intercambiadores de 11.8 Combustión de corrientes residuales en el proceso de
calor utilizados 211 estireno 249 11.9
9.10 Diagrama de flujo del proceso sin columna Ciclo de bomba de calor 251 11.10
D215 212 9.11 Posibles arreglos de bombas de calor en un
Hojas de trabajo del intercambiador de calor para el columna de destilación 251
proceso sin D215 213 10.1 11.11 Simulación de bomba de calor en el proceso de
Estrategia para el reciclaje de un flujo de material 216 10.2 Método separación de propileno/propano 252 11.12
de convergencia y parámetros en Aspen Plus 217 10.3 Curvas compuestas 254 11.13 Grandes
Proceso de acetato curvas compuestas 254 11.14 Costo objetivo
de etilo con bucles de reciclaje 218 10.4 Especificación de diseño antes del diseño HEN 255 11.15 Ejemplo de un diagrama
para el reciclaje de ácido acético 219 10.5 Conexión de cuadrícula 255 11.16 Sección del proceso de
de flujos de reciclaje deshidrogenación de nheptano para
de material en Aspen HYSYS 222 10.6 Página de PPA 256 11.1 7 Inicio del análisis energético 256
parámetros de 11.18 Objetivos y potencial de ahorro de
convergencia de Aspen HYSYS 222 10.7 Diagrama de flujo energía 257 11.19 Resultados del análisis energético
del proceso de proporcionados directamente en el entorno de
producción de estireno con reciclaje de etilbenceno 223 simulación 257 11.20 Resumen del rendimiento y
10.8 Diagrama de flujo simplificado de la diagrama de cuadrícula para el caso de simulación
producción de acetato de vinilo 227 10.9 Definición de un actual 258
bloque de cálculo para 11.21 Detalles de las conexiones del intercambiador de calor
composición de ácido acético 227 10.10 Bloque de cálculo en AEA 258
para flujo molar 11.22 Curva compuesta y objetivos estimados por el PPA
de acetileno 228 10.11 Sensibilidad bloque para la optimización de 259 11.23 Encontrar
la fracción dividida 228 10.12 Resultados de la un arreglo óptimo de HEN 259 11.24 Diseño casi óptimo de HEN
optimización de la recomendado por AEA 260 11.25 Diagrama
fracción dividida 229 10.13 Definición de la herramienta de de flujo del proceso después de la
optimización 229 10.14 Resultados de la optimización implementación optimizada de HEN 260
230 10.15 Cálculo del requisito de vapor de
las columnas de destilación 231 10.16 Especificación 12.1 Diagrama de flujo de la compresión de gas de síntesis 265
de un bloque calentador 231 12.2 Diagrama de flujo de la sección de reacción del
10.17 Simulación de un ciclo de refrigeración 232 10.18 proceso de metanol 265
Especificación de diseño para el requisito de amoníaco 233 12.3 Diagrama de flujo de la sección de destilación de
10.19 Definición de la reacción de conversión en el proceso de metanol 266 12.4
Aspen HYSYS 234 Pasos del análisis económico en la APEA integrada
10.20 Diagrama de flujo para el cálculo de los requisitos 267
de combustible 12.5 Ingreso de costos de materias primas en una
235 simulación de Aspen 272 12.6
Diagrama de flujo de efectivo calculado para el
10.21 Definición del operador Conjunto 235 proceso de metanol 274
Machine Translated by Google
12.7 Procesar la evaluación económica con un clic 275 12.8 Uso de 13.32 Cambio del método termodinámico para la columna
Economics Active para la evaluación económica 276 12.9 de vacío 305 13.33 Perfiles de
Hoja de Excel de temperatura y flujo en la columna de destilación
análisis de inversión en el método Economics Active 276 atmosférica 306
13.34 Visualización de curvas de destilación 306 13.35
12.10 Envío de una simulación a APEA 277 12.11 Curvas ASTM D86 de productos y alimentación de
Mapeo, dimensionamiento y evaluación de bloques de columna atmosférica 306 13.36
operaciones unitarias en APEA 277 12.12 Comparación de curvas de destilación experimental y modelo
Edición de dimensionamiento de equipos y visualización de ASTM D86 307
resultados de evaluación económica 278 13.37 Curvas de destilación ASTM D86 de productos de
13.1 Generación de pseudocomponentes basados en rangos columna de vacío 307 13.38
de puntos de ebullición 286 13.2 Modelos de refinería de Aspen HYSYS V.9 308 13.39 Esquema
Importación de una lista de componentes de reacción de hidrocraqueo VR 309 13.40 Esquema de un proceso
adecuados para petróleo típico de hidrocraqueo en lecho fijo con dos
286 13.3 Selección de un paquete de fluidos y apertura del reactores 310 13.41 Hidrocraqueo en lecho ebullido 310
administrador de ensayos de petróleo 13.42 Generación de pseudocomponentes a partir de
287 13.4 Adición de un nuevo ensayo de petróleo 287 un conjunto de ensayos 312
13.5 Selección de un ensayo basado en el país de origen
288 13.6 Ensayo 13.43 PFD de hidrocraqueo VR 313 13.44
caracterizado 288 13.7 Formateo de Selección de un modelo de usuario en Aspen Plus 314 13.45
gráficos de características del ensayo 289 13.8 Curva de Especificación del modelo de usuario 314 13.46 Hojas de
destilación TBP del petróleo crudo seleccionado 289 13.9 archivo de Excel para la comunicación con
Contenido Aspen 314
de PNA del crudo seleccionado 290 13.10 Apertura de una refinería 13.47 Rendimientos de productos de hidrocraqueo
de petróleo VR calculados por el modelo cinético 315
caso de simulación de fraccionamiento 13.48 Efecto de la temperatura del reactor en la
290 13.11 Una lista típica de componentes de refinería curva de destilación del producto da 315 13,49
de petróleo en Aspen Plus 291 Selección de un modelo FCC en Aspen
13.12 Especificación de la curva de destilación 292 13.13 HYSYS 316
Pestaña Composición final ligera 292 13.14 13.50 Conectación de una unidad FCC con alimentación
Generación de pseudocomponentes 293 13.15 Resultados interna solo 316
de la generación de pseudocomponentes 293 13.16 Conexión de un 13.51 Definición del alimento FCC 317 13.52
ensayo de petróleo a la simulación HYSYS 295 Especificación de alimentación de FCC 317
13.53 FCC Especificación del elevador 318
13.17 Especificación de columna por condensador 13.54 Especificación de regenerador FCC 318 13.55
temperatura y punto de corte del 95 % 296 Resultados de simulación de FCC 319 14.1
de corriente lateral 304 13.30 Especificación Variables en el bloque Calculator 328 14.12 Cálculo de
de pumparound 304 13.31 Especificación de columna de vacío 305 rendimientos en base húmeda 328
Machine Translated by Google
14.13 Especificación del separador de sólidos 329 15.13 Resultados de la geometría interna de la columna 363
14.14 Contracorriente (corriente ascendente) y 15.14 Perfil de concentración de HCl en la fase gaseosa 364
gasificadores de lecho móvil en cocorriente (b 15.15 Perfiles
downdraft) 332 de temperatura de las fases líquida y gaseosa 364 16.1
14.15 Gasificador de lecho fluidizado circulante 332 14.16 Inicio de la
Gasificador de flujo arrastrado 333 14.17 simulación de un proceso con polímeros 369 16.2 Definición
Esquema simplificado de gasificación de de un nuevo
combustible sólido para la producción de metanol 333 conjunto de unidades 370 16.3 Lista de
14.18 Gasificación RDF PFD 336 14.19 componentes para PS libre polimerización por radicales en
Variables definidas en el bloque de cálculo del masa 370 16.4 Especificación
gasificador 336 14.20 del tipo de segmento y grupo de
Cálculo de rendimientos másicos y secuencia de atributos del polímero 371 16.5
cálculo 337 14.21 Contenido de Atributos de los componentes del polímero 371 16.6
componentes principales y Definición de parámetros escalares de componentes
temperatura del gasificador versus oxígeno para puros 373 16.7
Relación de masa RDF 339 Selección del tipo de reacción de polimerización 374 16.8
14.22 Conversión, caudal volumétrico de gas específico, poder Especificación de las especies de reacción de
calorífico inferior (LHV) del gas y temperatura del polimerización
gasificador 340 14.23 375 16.9 Creación de reacciones de polimerización 375 16.10
Contenido de alquitrán de gas y fracción molar de CO2 , Lista de reacciones de polimerización por radicales
temperatura del gasificador y conversión 340 libres de estireno a granel 376 16.11
14.24 Efecto de la relación de masa de vapor a RDF en el gas Parámetros cinéticos de estireno a granel polimerización
composición y temperatura del reactor 340 14,25 Pirólisis por radicales libres 376 16.12 Parámetros
de biomasa PFD 343 14,26 Configuración de especiales de iniciación 376 16.13 Polimerización por
la mejora de bioaceite radicales libres en masa de estireno PFD 377 16.14 Definición
columna de destilación 343 del bloque
14.27 Especificación del extractor lateral 343 15.1 DS para ajustar la composición de alimentación constante 377
Selección del asistente de electrolitos 348 15.2 16.15 Lista de variables
Componentes básicos y opciones de generación de reacción definidas para análisis de sensibilidad 378 16.16 Cálculo
y especies generadas y de conversión y
reacciones 348 tabulado
15.3 Selección del enfoque de simulación de electrolitos 349 15.4 variables en el análisis de sensibilidad 378
Revisión de las constantes de equilibrio 349 15.5 Lista final 16.17 Visualización de los resultados del tamaño de la cadena
de componentes, incluidas las especies de electrolitos del polímero 380 16.18 Curvas de distribución del
350 15.6 PFD para la eliminación tamaño de la cadena 381 16.19 Efecto de la temperatura del
de HCl del vinilo segundo reactor sobre la conversión en los reactores 381
corriente de cloruro 353 16.20 Efecto de la temperatura del segundo reactor en el PDI
15.7 Efecto del flujo másico de solvente en la 381 16.21 Efecto de
eliminación de HCl 355 la temperatura del segundo reactor
15.8 Efecto de la temperatura en la eliminación de HCl 355 15.9 sobre polímero MWW 382
Diagrama de flujo del proceso de eliminación de gases ácidos 16.22 Efecto de la temperatura del segundo reactor
por solución acuosa de aminas 358 15.10 sobre polímero MWN 382
Selección del enfoque basado en la velocidad 361 15.11 16.23 Efecto de la temperatura del segundo reactor sobre el
Especificación de los componentes internos de la flujo másico del PS producido y reciclado 382
columna 362 15.12 Especificación de los parámetros basados en la velocidad 362
Machine Translated by Google
xiv
Sobre el Autor
Juma Haydary es profesor asociado de ingeniería química ing y separación de mezclas multicomponentes. Sus
en el Departamento de Ingeniería Química y Ambiental de la principales actividades de investigación, además de la
Universidad Tecnológica de Eslovaquia en Bratislava, simulación y optimización de procesos, incluyen procesos de
Eslovaquia. Tiene un doctorado en ingeniería química y pirólisis, craqueo y gasificación. Ha realizado contribuciones
control de procesos de la misma universidad. Está involucrado de gran alcance en estos campos durante los últimos 15
en la enseñanza del diseño de procesos químicos asistidos años. Hasta 2017, ha sido el investigador principal de nueve
por computadora, operaciones unitarias en ingeniería química. proyectos de investigación y educación.
Machine Translated by Google
xvii
Prefacio
El libro proporciona un estudio práctico, completo y orientado a se resaltan los equipos y se muestran los métodos para su
ejemplos sobre el diseño y la simulación de procesos químicos selección y costeo.
utilizando el software de simulación Aspen Plus y Aspen La tercera parte del libro describe el diseño de nuevas plantas
HYSYS. Este es el tema de interés actual de los ingenieros de y la simulación de plantas existentes donde se utilizan
procesos químicos, diseñadores y desarrolladores de procesos, productos químicos convencionales y mezclas de materiales
ingenieros energéticos, evaluadores económicos de procesos con composiciones medibles. En esta sección, la integración
e investigadores involucrados en ingeniería química. de materiales, el análisis energético y la evaluación económica
El libro tiene como objetivo proporcionar principios de diseño se describen y muestran en ejemplos en Aspen Plus y Aspen
de ingeniería química y, al mismo tiempo, ejemplos de su HYSYS, respectivamente. También se muestran herramientas
aplicación utilizando el software de simulación más utilizado. específicas integradas en Aspen Plus y Aspen HYSYS para el
Como no es posible brindar una discusión exhaustiva de todos análisis energético y económico de los procesos. El uso de
los principios básicos de ingeniería química requeridos en un programas, como Aspen Exchanger Design and Rating (EDR)
libro de extensión razonable, este libro requiere el conocimiento para el diseño detallado de intercambiadores de calor de
de los principios básicos de ingeniería química. El libro puede carcasa y tubos, Aspen Energy Analyzer (AEA) para el análisis
servir como libro de texto para un curso superior o de posgrado de puntos críticos y el diseño de redes de intercambiadores de
en ingeniería química, así como también como una referencia calor, y Aspen Economic Evaluation para el mapeo ,
para profesionales involucrados en ingeniería de procesos dimensionamiento y evaluación económica de los procesos, se introduce en el l
químicos o investigación y desarrollo de ingeniería química. La última parte del libro trata sobre el diseño de plantas y la
simulación de procesos utilizando componentes no
El libro está dividido en cuatro secciones: Parte I Introducción convencionales. En muchos procesos industriales, no se
al diseño y simulación, Parte II Diseño y simulación de conoce la composición exacta de los flujos de materiales. El
operaciones de una sola unidad, Parte III Diseño y simulación diseño y la simulación de tales procesos requieren enfoques específicos.
de plantas: componentes convencionales, y Parte IV Diseño y Pertenecen a este grupo la refinación de petróleo y la mejora
simulación de plantas: componentes no convencionales. La de los productos de destilación del petróleo crudo, la combustión,
sección Introducción al diseño y la simulación describe la gasificación, la pirólisis y el secado de diferentes sólidos, y
metodologías básicas para el diseño asistido por computadora muchos otros procesos. Los dos últimos capítulos de esta parte
seguidas de una descripción de los pasos básicos de la están dedicados a la modelización de procesos con electrolitos
simulación de procesos en Aspen Plus y Aspen HYSYS. Los y procesos de polimerización. Los ejercicios no resueltos se
temas principales de esta sección son la síntesis de conceptos agrupan en tres partes ubicadas al final de las Partes II, III y IV.
de procesos, la recolección de datos, el equilibrio termodinámico Las soluciones de todos los ejercicios en Aspen Plus y Aspen
de fases, las propiedades físicas, los componentes HYSYS están disponibles como materiales complementarios.
convencionales y no convencionales y los datos de reacciones Hay varios libros sobre diseño de procesos. Sin embargo, la
químicas. La Parte II está dedicada al diseño y simulación de mayoría se centran más en principios teóricos con ejemplos
operaciones de unidades individuales individuales. Luego de limitados de su aplicación o se limitan a tipos específicos de
representar un modelo matemático de la operación de cada procesos. El presente libro tiene como objetivo combinar los
unidad, tales como reactores, separadores, intercambiadores principios teóricos con ejemplos reales resueltos por el software
de calor, etc., se resuelven ejemplos de problemas reales utilizando Aspen Plus y Aspen
de simulación HYSYS,que
más utilizado respectivamente.
emplea la simulación de
Requisitos específicos para el modelo de un tipo dado de procesos en estado estacionario utilizando módulos secuenciales.
Machine Translated by Google
xxviii Prefacio
acercarse. Además de los procesos con productos químicos orgánicos con una estructura gráfica completamente diferente en comparación con
convencionales, también se incluyen en este libro procesos con materiales las versiones anteriores utilizadas en los libros publicados hasta el momento.
más complejos, como sólidos, mezclas de aceites, polímeros y electrolitos. Otra particularidad de este libro es que se presentan tanto Aspen Plus
Los ejemplos se resuelven utilizando una versión más nueva del software como Aspen HYSYS, lo que permite su comparación.
Aspen (ASPEN One 9)
Machine Translated by Google
xxix
Expresiones de gratitud
El autor desea agradecer a las personas e instituciones que hicieron en el libro. Expreso mi agradecimiento al Departamento de Ingeniería
posible este libro. Deseo agradecer a Barbora Duda´sova, Jakub Hus ar y Química y Bioquímica de la Universidad Tecnológica de Eslovaquia por
ˇ ´
Patrik
ˇ proporcionar instalaciones y apoyo. Mi profunda gratitud pertenece a
Suhaj por su ayuda en la revisión y corrección del manuscrito y Juraj Aspen Technology, Inc. (200 Wheeler Road, Burlington, MA 018035501,
´
Labovsky por su apoyo técnico. Deseo expresar mi profunda gratitud y EE. UU., www.aspentech.com) por permitirme usar algunas tablas de los
agradecimiento a muchos de mis colegas, nuestros estudiantes de manuales de usuario de Aspen plus y Aspen HYSYS.
doctorado y estudiantes de pregrado por proporcionar comentarios sobre
los materiales incluidos.
Machine Translated by Google
xxi
abreviaturas
abdominales
Acrilonitrilo butadieno estireno LHHW Langmuir–Hinshelwood–Hougen–
AAE Análisis de energía de Aspen modelo watson
CÓMO Analizador económico de procesos Aspen LLDPE Polietileno lineal de baja densidad
API Instituto Americano de Petróleo LLE Equilibrio líquidolíquido
COMO YO Sociedad Americana de Mecánica LP Baja presión
Ingenieros GLP Gas de petróleo licuado
PA Punto de ebullición MDEA metildietanolamina
RRC Anillos de cerámica Raschig COSA Amina de monoetanol
CS Acero carbono MALLA Materialequilibriosumacalor
CSTR Reactor de tanque agitado continuo MWN Número de peso molecular promedio
DDM Ndodecilmercaptano No convencional
CAROLINA DEL NORTE
xxxii Abreviaturas
xxxii
simbolos
Capitulo 2
A −1
factor de frecuencia [m3 mol−1 s ]
a Parámetro de la ecuación cúbica de PengRobinson
B Parámetro de ecuaciones viriales
b Parámetro de la ecuación cúbica de PengRobinson
Cd Parámetro de ecuaciones viriales
Exponente de temperatura
Y Energía de activación [J mol−1 ]
F Número de grados de libertad
f Fugacidad [Bien]
d factor de compresibilidad
Volatilidad relativa
Factor de separación
Coeficiente de actividad
Capítulo 3
xxxiv Símbolos
hacer
Diámetro exterior del tubo [metro]
F Factor de corrección
Capítulo 4
−1
cp Capacidad calorífica isobárica [J mol−1 K ]
−1
cv Capacidad calorífica isocórica [J mol−1 K ]
ffg Frecuencia de rotación [s−1 ]
r Fricción [Pa m−1 ]
−2
Aceleración gravitacional [m s ]
H Cabeza [m]
h.com cambio de entalpía por mol de gas [J mol−1 ]
Hs Cabeza estática hidráulica [m]
Δh Cambio de entalpía real por mol de gas [J mol−1 ]
L Longitud de la tubería [m]
norte
exponente politrópico
n55 Velocidad específica de succión [rpm]
NPSHA Altura neta de succión positiva disponible [m]
NPSHR Cabezal de succión positivo neto requerido [m]
PAG
Presión [Bien]
PAG ◦ Presión de vapor del líquido en condiciones de entrada. [Bien]
Símbolos xxxv
Pi el numero de ludolf
Densidad del líquido [kg m−3 ]
metro
Velocidad del eje [kg m−3 ]
angular de la densidad del fluido [rad s−1 ]
Capítulo 5
comer
Actividad del componente i
B Constante de ecuación para calcular la constante de equilibrio o el término de [K]
absorción de una reacción química
b Parámetro de la ecuación para calcular la constante de equilibrio de la reacción
química.
C Constante de ecuación para calcular la constante de equilibrio o el término de
absorción de una reacción química
ΔfHi
◦ Calor estándar de formación [J mol−1 ]
ΔrHn Entalpía de reacción de reacción n [J mol−1 ]
xxxvi Símbolos
Cuando
Constante de equilibrio de la reacción química.
a Constante del término de adsorción
ni, o
Flujo molar inicial (o moles) del componente i [mol s −1 (o mol)] (o
ni, ni
Flujo molar (o moles) del componente que reaccioné [mol s −1 mol)] (o
reaccionado , t
Flujo molar (o moles) del componente i en la mezcla de reacción en el tiempo t [mol s −1 mol)]
es lo
Flujo molar del componente i en la corriente de salida k [mol s −1 ]
Flujo molar de la corriente de entrada j [mol s −1 ]
que
haciendo
Presión del componente i [Bien]
Flujo de calor −1
qq [J s ]
Salida
Constante de gas −1
R [J mol−1 K ]
Velocidad de reacción química −1
r [mol m−3 s ]
◦ Entropía absoluta estándar [J mol−1
−1
K ]
Y
T Temperatura [K]
Capítulo 6
B Abajo
D Destilar
F Alimentar
Símbolos xxxvii
q Cantidad de líquido
R Relación de reflujo externo
Relación de reflujo mínima
rmin
−1
Velocidad de reacción n en la etapa j [J m−3 s ]
rjn
SL Líquido lateral
SV Vapor lateral
T Temperatura [K]
EN Vapor
Retención volumétrica de líquido en la etapa j [m3 ]
(VLH) j
X Fracción molar de la fase líquida Fracción
estequiométrico del
componente i en la
Capítulo 7
punto de acceso
Área de superficie de una sola partícula [m2 ]
un Altura de entrada del ciclón [m]
b Ancho de entrada del ciclón [m]
Concentración de sólidos en la corriente de entrada [kg m−3 ]
c0
Concentración de sólidos en la corriente de salida [kg m−3 ]
cout
Corriente continua
Diámetro del cuerpo del ciclón [m]
Diámetro de partícula [m]
Dp
GRAMO Gas
H entalpía de masa [J kg−1 ]
k Parámetro de configuración geométrica
L Longitud del secador [m]
−1
METRO
Tasa de evaporación de una sola partícula [kg s ]
P.ej
Número total de partículas
norte
Parámetro de eficiencia del ciclón
PAG Caída de presión [Pa]
Caudal volumétrico de gas −1
q [m3 s ]
S Sólido
T Temperatura [K]
−1
Afuera Velocidad del gas de entrada [m s ]
Machine Translated by Google
xxxviii Símbolos
En Agua
X Contenido de humedad en el sólido.
Xcr
Contenido crítico de humedad
GRAMO
Densidad de gases [kg m−3 ]
pag
Densidad de particula [kg m−3 ]
Capítulo 8
Capítulo 10
q Presión Parámetro de
qcom
aceleración Calor de la [kcal kmol–1]
S combustión del gas natural Parámetro calculado a
T partir de la [K]
X ecuación 10.4 Temperatura Parámetro
Machine Translated by Google
Símbolos xxxix
Capítulo 11
Capítulo 12
Capítulo 13
k Equilibrio constante
k Constante de la cinética de reacción [varía] [s]
t tiempo
En Fracción de peso
Coeficiente de actividad
Fugacidad [Bien]
Coeficiente de fugacidad
Machine Translated by Google
Símbolos xl
capitulo 14
−1
Capacidad calorífica [kJ K ]
CP
CON h Conversión de combustible [%]
entalpía específica [kJ kg−1 ]
ΔH◦ Entalpía estándar de reacción [kJ mol−1 ]
Δch Entalpía específica de combustión/calor de combustión [kJ kg−1 ]
n1 Relación molar de H2 y CO
R Relación de peso de oxígeno a RDF
R Número de átomos de azufre basado en un solo átomo de carbono en combustible
sólido
R1 Ración de peso de vapor a RDF
r
Número de átomos de azufre basado en un solo átomo de carbono en
acepta
CON
Número de átomos de nitrógeno basado en un solo átomo de carbono en combustible
sólido
'z Número de átomos de nitrógeno basado en un solo átomo de carbono en el alquitrán
Capítulo 15
Símbolos xli
−1
S Entropía molar [kJ mol−1 K ]
T Temperatura termodinámica [K]
X Fracción molar
factor de no aleatoriedad
Parámetro binario moléculamolécula
Energía libre de Mole Gibbs [kJ mol−1 ]
capitulo 16
A Exponente
A Agente
B Exponente
ESO Concentración de agente [mol m−3 ]
CC Concentración de iniciador [mol m−3 ]
CINI Coiniciador
Dn polímero muerto
Dn+m Cadena muerta única con un par de segmentos de cabeza a cabeza
Y Energía de activación [J mol−1 ]
novia
factor de efecto gel
I Iniciador
IDENTIFICACIÓN
Descomposición térmica del iniciador
k Tarifa constante Según el orden de la reacción.
METRO monómero
norte
Longitud de la cadena de polímero
No. Indicando la formación de 1 o 2 radicales
PAG
Cadena polimérica
PAG Presión [Bien]
−1
R Constante de gas [J mol−1 K ]
R Radical
RAD radical primario
ORZUELO monómero [m3 .mol−1 ]
T Temperatura [K]
FRENTE A
Transferencia de cadena a agente
CT Terminación por combinación
DE Iniciación térmica
TM Transferencia a monómero
xliii
www.wiley.com/go/Haydary/ChemDesignSimulation_Aspen
Los materiales complementarios contienen la solución de los ejercicios presentados al final de las partes II, III y IV del libro. Los
ejercicios se resuelven utilizando Aspen Plus y Aspen HYSYS, respectivamente. Las soluciones se proporcionan en forma de archivos
compuestos de Aspen Plus (.apwz) y archivos de casos de simulación HYSYS (.hsc). Para abrir estos archivos, es necesario instalar
Aspen Plus V9 y Aspen HYSYS V9 o versiones superiores en su computadora. Para utilizar archivos complementarios, siga los
siguientes pasos:
– Asegúrese de que Aspen Plus V9 o versión superior para abrir archivos compuestos de Aspen Plus (.apwz) y Aspen HYSYS V9 o
versiones superiores para abrir archivos de casos de simulación HYSYS (.hsc) estén correctamente instalados en su computadora.
– Asegúrese de tener una licencia válida para usar Aspen Software.
– Descarga el archivo que deseas revisar y ábrelo de manera regular.
– En el entorno de Propiedades comprobar:
◦ lista de componentes
◦ método termodinámico utilizado
◦ reacciones en caso de archivos
◦ HYSYS ensayo, mezclas y petrocaracterización en caso de refinería casos resueltos por Aspen Plus
◦ ensayos de petróleo y manejo de aceite en caso de casos de refinería resueltos por Aspen
HYSYS ◦ propiedades de componentes no convencionales en caso de ejercicios con componentes no
convencionales ◦ polímeros en caso de ejercicios
con polímeros ◦ química en caso de electrolitos resueltos por Aspen Plus
◦ Resultados, resultados de la transmisión y páginas de resultados de la transmisión (personalizados) debajo de cada bloque de Aspen Plus,
incluida la sensibilidad, las especificaciones de diseño.
y bloques de calculadora ◦ Hoja de trabajo y páginas de rendimiento de cada bloque de Aspen
HYSYS ◦ Resumen de resultados en el caso de Aspen Plus
◦ Libro de trabajo, estudios de casos, secuencias, etc. en el caso de Aspen HYSYS
– En caso de ejercicios de análisis de energía, verifique los resultados del análisis de energía en Análisis de energía
– Para verificar los resultados de la evaluación económica, active Economía y pase por mapeo, dimensionamiento y evaluación
pasos.
Machine Translated by Google
Parte I
En general, los ingenieros de procesos químicos se ocupan de dos incluida la selección del proceso químico, la tecnología, las
tipos de tareas: el diseño de un nuevo proceso y la simulación de condiciones del proceso, la recopilación de los datos requeridos, la
un proceso existente. Estas tareas pueden ser simples o muy emisión de diagramas de flujo del proceso, la selección, la
complejas. Para resolver algunos problemas simples, se puede especificación y los cálculos de ingeniería química del equipo y la
usar el cálculo manual. La ventaja del cálculo manual es un estimación preliminar de costos. La segunda fase es el diseño
conocimiento más profundo del problema. Sin embargo, en un básico de la planta, incluido el diseño mecánico detallado del
problema complejo, a menudo se requiere la solución de muchos equipo, el diseño detallado de las estructuras eléctricas y civiles,
miles de ecuaciones. Por lo tanto, el cálculo manual de este tipo de y las tuberías y los servicios auxiliares. Los pasos incluidos en la
problemas en tiempo real es prácticamente imposible y los primera fase suelen ser realizados por ingenieros químicos,
simuladores de procesos son una herramienta insustituible. Tanto mientras que los incluidos en la segunda fase son realizados por
las tareas de diseño como las de simulación requieren enfoques especialistas en diseño. En la Figura 1.2, los pasos de la primera
específicos. El diseño de un proceso químico comienza con un fase se muestran en detalle en rectángulos con bordes sólidos y
requisito para un producto y pasa por diferentes pasos de diseño. crean los temas de estudio en este libro. Los pasos del siguiente
Una tarea de simulación comienza con un requisito para la diseño básico solo se mencionan aquí y se muestran en rectángulos
modificación u optimización del proceso y continúa con un análisis con bordes discontinuos.
del estado del arte existente. El diseño de un proceso comienza con una idea para producir
Este capítulo introduce el concepto de este libro para el diseño un nuevo producto o para mejorar un producto existente. El
de procesos nuevos y la simulación de procesos existentes. requerimiento de un nuevo producto puede resultar de la demanda
Se explican los niveles de jerarquía, la profundidad y los pasos del mercado. Si el análisis inicial indica que la idea puede
básicos del diseño de procesos químicos. También se analizan el convertirse en un proyecto, se inicia una revisión de sus aspectos
concepto de la química del proceso, las variantes tecnológicas, la químicos, tecnológicos y económicos. La evaluación de la química
recopilación de datos y los pasos de desarrollo del diagrama de del proceso y las posibles materias primas para la producción del
flujo del proceso del diseño conceptual. Este capítulo también se producto requerido es una de las primeras investigaciones. Las
ocupa de los programas de simulación de procesos, el enfoque reacciones químicas y los catalizadores se seleccionan en el
modular secuencial y el enfoque orientado a ecuaciones. La última contexto de las condiciones locales, la disponibilidad de materias
parte está dedicada a iniciar una simulación con Aspen Plus y Aspen HYSYS.
primas y los efectos ambientales, económicos, de seguridad y de
salud. Un concepto químico a menudo se analiza simultáneamente
con el concepto tecnológico y la evaluación económica aproximada.
1.1 Diseño de procesos Se investigan los tipos de reactores, la fase de las reacciones, los
tipos de unidades de separación y otras condiciones tecnológicas
Para el desarrollo de un diseño de ingeniería química, primero se en relación con las especificaciones económicas, ambientales y
debe definir el nivel de jerarquía del problema de diseño. Bajo el locales. El siguiente paso es la recopilación de datos necesarios
término de diseño en ingeniería química, podemos imaginar desde para el diseño del proceso; los datos que se requieren con mayor
el diseño de una sola operación hasta el diseño completo de una frecuencia son los datos de las propiedades del material, el
fábrica química. Los niveles de jerarquía para las tareas de diseño equilibrio de la reacción química y los datos cinéticos, así como los
en ingeniería química se muestran en la Figura 1.1. El propósito de datos del equilibrio de fase.
este libro es presentar los principales aspectos del diseño de Un diagrama de flujo del proceso se desarrolla en una serie de
operaciones unitarias, así como plantas químicas. iteraciones iniciadas por una configuración simple de los bloques
de equipos principales. La selección de un modelo de equilibrio
Otro aspecto del diseño de ingeniería química es la profundidad de fase termodinámico adecuado es el momento crucial en el
del diseño. Por lo general, el diseño se divide en dos grandes diseño del equipo de separación, como columna de destilación,
fases. La primera fase es un diseño conceptual. separadores, absorbedores, extractores y otros.
Diseño y simulación de procesos químicos: aplicaciones Aspen Plus y Aspen HYSYS, primera edición. Juma Haydary.
© 2019 John Wiley & Sons, Inc. Publicado en 2019 por John Wiley & Sons, Inc.
Sitio web complementario: www.wiley.com/go/Haydary/ChemDesignSimulation Aspen
Machine Translated by Google
Recopilación de datos necesarios para el diseño (datos de propiedades del material, equilibrio
químico y datos cinéticos, datos de equilibrio de fases, etc.)
Construcción de la planta
los procesos físicos son el tema del diseño, por ejemplo, la separación alternativas de separación,
de un componente o un grupo de componentes de una mezcla natural. alternativas para el reciclaje de flujo de materiales,
La separación primaria de petróleo crudo es un buen ejemplo de tales concepto de integración energética y
procesos. efectos ambientales, de salud y seguridad.
Si el nuevo producto es el resultado de reacciones químicas, el
proceso de diseño comienza con la búsqueda de reacciones químicas En el paso de selección del reactor, una decisión muy importante
que conduzcan a este nuevo producto. Las condiciones locales iniciales que se debe tomar es la elección entre procesos continuos y por lotes,
juegan un papel importante en la selección de la química del proceso. que se ve influenciada principalmente por el carácter del proceso y la
El uso de una planta existente, la modificación de una planta existente capacidad de producción. Para muchos procesos con capacidades
y el diseño de una planta completamente nueva son tres situaciones considerables, se prefieren los procesos continuos. La fase de reacción
iniciales diferentes que afectan la selección de una química apropiada es otro tema importante.
para el nuevo producto. La búsqueda de una química potencial comienza La conversión y la eficiencia del proceso pueden verse fuertemente
con libros, libros de texto y enciclopedias. Se puede encontrar afectadas por la fase de reacción. Las condiciones de reacción, como
información más detallada en revistas, patentes y otras publicaciones. la temperatura, la presión y la transferencia de calor entre el reactor y
el entorno (condiciones isotérmicas o adiabáticas) son otros parámetros
importantes de los reactores. Las temperaturas y presiones más
El primer resultado de una búsqueda inicial es la respuesta a la preferibles para los reactores químicos son aquellas cercanas a la
pregunta, si se conoce la química exacta de la producción del producto temperatura ambiente y la presión atmosférica; muchos procesos
requerido. Las siguientes respuestas a esta pregunta son posibles: requieren diferentes condiciones. Sin embargo, para seleccionar
temperaturas y presiones altas o bajas tiene que haber un
En muchos casos, hay disponible más de una elección de materias Los reactores continuos más utilizados son CSTR (reactor de
primas para un nuevo producto. La materia prima apropiada y también tanque agitado continuo) y reactor tubular. Para modelar los reactores
el tipo de catalizador se seleccionan considerando los mejores CSTR, a menudo se aplica la teoría de la mezcla ideal; mientras que en
indicadores económicos que satisfacen las condiciones ambientales y el modelado de reactores tubulares se emplea la teoría del flujo pistón.
de seguridad. Una visión general económica y ambiental inicial ayuda Para cada variante de reactor para una tecnología dada, también se
a excluir la mayoría de los métodos inconvenientes de una mayor deben evaluar los aspectos ambientales, de seguridad y de salud.
investigación. La química con los mejores indicadores económicos
puede no ser siempre la misma; se ve fuertemente afectado por las Los productos de reacción suelen estar en forma de mezclas
condiciones locales, como la disponibilidad de materias primas, la homogéneas o heterogéneas, y el producto requerido debe separarse,
posibilidad de utilizar las tecnologías e infraestructuras existentes, las lo que normalmente requiere una serie de operaciones. Para mezclas
reglamentaciones ambientales, la disponibilidad de fuentes de energía, heterogéneas, procesos de separación como filtración, separación
etc. ciclónica, precipitación, sedimentación, etc. estan empleados.
El objetivo es seleccionar alternativas adecuadas sin una simulación 1.4.2 Datos de equilibrio de fase
detallada de todas las alternativas posibles. La tecnología debe diseñarse
Para un diseño satisfactorio del equipo de separación y reacción, la
con la tasa de reciclaje óptima de la materia prima, teniendo en cuenta
calidad del modelo para los cálculos de equilibrio de fase tiene un efecto
también las cuestiones ambientales, de seguridad y de salud. Tenga en
crucial. La calidad del modelo viene dada por su capacidad para describir
cuenta que el reciclaje de flujos de materiales en toda la gama puede no
el proceso real.
ser siempre el método más eficaz.
Los datos experimentales de equilibrio de fase permiten la verificación
del modelo termodinámico utilizado (ver Ejemplo 2.7). Se han publicado
Las variaciones de la integración de la energía del proceso también
datos experimentales de vaporlíquido y líquidolíquido para miles de
deben evaluarse como parte del concepto de tecnología. En este paso se
sistemas binarios en bases de datos como DECHEMA y el Instituto
estudian las alternativas básicas de la integración energética del proceso;
Nacional de Estándares y Tecnología (NIST). Sin embargo, para miles de
El diseño detallado de la integración energética y las redes de
otros sistemas binarios, estos datos no están disponibles. Los parámetros
intercambiadores de calor son el tema de la integración del proceso
de los modelos de equilibrio de fase basados en coeficientes de actividad
acompañada de la simulación de diferentes alternativas.
se pueden calcular mediante métodos de contribución como los
Coeficientes de actividad de grupos funcionales UNIQUAC (UNIFAC); en
Después de la evaluación de las variaciones de la tecnología, incluidos
un diseño de proyecto real, es necesaria la verificación experimental del
los aspectos ambientales, de seguridad y de salud, el diseñador tiene una
modelo de equilibrio de fase.
breve lista de casos que deben simularse mediante un simulador de
procesos para tomar decisiones adicionales.
La calidad de una simulación depende en gran medida de la calidad de El modelado de reactores químicos requiere información sobre la
los datos y los parámetros del modelo utilizado. La calidad y disponibilidad estequiometría de la reacción, la constante de equilibrio y los parámetros
de los datos son dos de los temas más desafiantes en muchas cinéticos de las reacciones químicas. En algunos casos, la conversión de
simulaciones. El software de simulación utilizado en el diseño de la reacción es conocida, por ejemplo, en reacciones muy rápidas, donde
procesos químicos contiene muchas bases de datos de datos de la conversión total se logra inmediatamente.
propiedades de materiales y datos de equilibrio de fase, particularmente Las constantes de equilibrio se pueden calcular minimizando la energía
para componentes convencionales. Sin embargo, en muchos casos, los libre de Gibbs. El software de simulación proporciona estos cálculos; sin
datos experimentales independientes son muy útiles en la verificación de embargo, los valores experimentales de la constante de equilibrio y su
los resultados. Además, para todas las sustancias no convencionales, dependencia con la temperatura brindan mejores resultados y pueden
faltan los datos de propiedades del material y de equilibrio de fase. usarse para la verificación de los datos calculados por la energía libre de
Además, los datos de equilibrio de fase no están disponibles para todos Gibbs.
minimización.
los posibles pares binarios de componentes convencionales. Otro tipo de
datos requeridos son los datos de equilibrio químico y los parámetros La ecuación de velocidad y los parámetros cinéticos de las reacciones
químicas utilizadas permiten calcular la conversión de la reacción y los
cinéticos de las reacciones químicas.
tamaños del reactor cuando se utiliza un modelo de reactor cinético.
Si se trata de componentes que no se califican como convencionales, Una tarea de simulación comienza con la definición del objetivo
como pseudocomponentes, ensayos, mezclas, sólidos no convencionales, (objetivos) en función del requisito de mejora de un proceso existente. El
etc. están presentes, se requiere cierta información sobre sus propiedades siguiente paso es estudiar los esquemas tecnológicos del proceso y la
para su caracterización. documentación para extraer la información requerida para la simulación.
Cuantas más propiedades se conozcan, mejor se podrá lograr la Los esquemas tecnológicos de proceso suelen ser muy detallados y
caracterización de un componente no convencional. contienen
Machine Translated by Google
diferentes tipos de información. Solo se puede utilizar cierta En el caso de una tarea de simulación, el punto de partida es
información para la simulación de procesos. Un ingeniero de el esquema tecnológico del proceso a analizar y simplificar para
procesos tiene que extraer la información necesaria y crear el la definición de los objetivos de la simulación. El PFD requerido
diagrama de flujo del proceso (PFD) basado en los objetivos de la para la simulación del proceso se deriva del esquema tecnológico
simulación y el esquema tecnológico del proceso. En el siguiente mediante la selección de equipos y flujos que pueden afectar el
paso, se deben recopilar datos de la operación de la planta. objetivo de la simulación.
Algunos datos de la planta se pueden utilizar como datos de Los programas de simulación usan módulos de bloques para
entrada al simulador y otros para la comparación de los datos del modelar diferentes tipos de equipos. A veces, el diagrama de flujo
modelo y de la planta real. Además, los datos descritos en la de simulación difiere del PFD real porque el equipo real en un
Sección 1.4 también deben recopilarse para una simulación programa de simulación se puede modelar mediante uno, dos o
adecuada del proceso. incluso más bloques de operaciones unitarias; o viceversa, un
Después de la preparación de un PFD simplificado y la bloque de operación en el software de simulación puede
recopilación de todos los datos necesarios, se puede realizar la representar más de un equipo real.
simulación del proceso con diferentes escenarios. Con base en Se pueden desarrollar dos modos diferentes de un PFD para la
los resultados de la simulación y su comparación con los datos simulación de procesos utilizando un simulador: modo activo y
operativos y el análisis de diferentes escenarios, se puede sugerir modo en espera. Cuando se usa el modo activo, el cálculo de las
la modificación del proceso. unidades se realiza simultáneamente con la creación del PFD
(después de instalar cada bloque de operación de la unidad, se
realiza su cálculo). En el modo en espera, primero se completa el
1.6 Desarrollo de Diagramas de PFD y luego se inicia el cálculo.
Flujo de Proceso
El desarrollo de un PFD para una tarea de diseño generalmente 1.7 Programas de Simulación de Procesos
comienza con diagramas muy simples sin redes de
intercambiadores de calor, modelos cinéticos del reactor, La simulación se define en (1) la simulación de la operación de un
materiales o integración de calor. La Figura 1.3 muestra un proceso o sistema del mundo real a lo largo del tiempo. En (2) se
ejemplo de un diagrama de flujo simplificado creado para el diseño proporciona una definición más detallada de la simulación de
de la producción de acetato de vinilo a partir de ácido acético y procesos: La simulación es un proceso de diseño de un modelo
acetileno. Después de calcular estos esquemas simplificados y operativo de un sistema y realización de experimentos con este
obtener el conocimiento sobre los antecedentes del proceso, los modelo con el propósito de comprender el comportamiento del
esquemas se pueden mejorar al incluir dos intercambiadores de sistema o de evaluar estrategias alternativas para el desarrollo o
calor laterales, modelos cinéticos de reactores y flujos de reciclaje la operación del sistema. Tiene que ser capaz de reproducir
de materiales y energía. En la Figura 1.4 se muestra un ejemplo aspectos seleccionados del comportamiento del sistema modelado
de un PFD más complejo para el mismo proceso (producción de acetatocon
de un grado aceptado de precisión.
vinilo).
C2H2REC
REACTOR ENFRIADOR
CALENTADOR SEP
MEZCLA
S1 S3
S4 S6
C2H2
DESPERDICIAR
PRODUCTO DC1
DC2
S7
S11
S12
S13
Figura 1.3 PFD simplificado para el diseño del proceso de producción de acetato de vinilo
Machine Translated by Google
ysimulació
químicos
procesos
Diseño
8de DIVIDIR
DESPERDICIAR
MIX4 C2H2REC
R2
C2H20 S8
F1
W4 W1
CACL2W
MEZCLA
REACTOR
O2 VT1
MIX3 O1 VT2 VT4
AA0 S1 S21 S22 S3 W8
F2
S4
PAO1
W5 S5
W7 q
1 de septiembre
S6 SEP2
R1 W6 W3 W2
VT3
CACL2W1
VTPARA
W11
VT5 W9
PRODUCTO MEZCLAR5
DC2
W10 RESIDUOS1
BOMBA
DC1
S14
S13 MEZCLA2
S7
S12
S10
S9
S11
Modelos de operación Modelos de En cada nivel, los conjuntos de ecuaciones no lineales se resuelven utilizando
unitaria propiedades físicas
el número de bucles de iteración y el procedimiento de solución interactivo. Más
detalles de los simuladores modulares se presentan en (3).
El primer significado indica el tipo de cálculo (diseño y simulación) como se procesos, incluida la producción de productos químicos, hidrocarburos,
explica en la introducción de este capítulo. El segundo significado es el modelado productos farmacéuticos, sólidos, polímeros, ensayos y mezclas de petróleo y
de un proceso por un simulador. otras aplicaciones. Aspen HYSYS es una herramienta de simulación muy potente
para aplicaciones de hidrocarburos, productos químicos y petróleo. Aspen Plus
La simulación de procesos (modelado) juega un papel crucial en todas las y Aspen HYSYS son partes del paquete de programas AspenOne lanzado por
actividades de ingeniería de procesos, incluida la investigación y el desarrollo, el Aspentech, Inc. Ambos software se aplican en este libro en la simulación de
diseño de procesos y la operación de procesos. diferentes tipos de procesos. Lanzado en diciembre de 2012 AspenOne V8.0, en
Una mayor extensión de la simulación de procesos incluye diferentes actividades 2016 AspenOne V9 y en 2017 AspenOne V10 representaron un salto significativo
basadas en computadora, como la dinámica de fluidos por computadora, pero en comparación con versiones anteriores, no solo en gráficos sino también en
el tema de la simulación en este libro es el diagrama de flujo de los procesos funcionalidad. Los libros anteriores sobre la aplicación Aspen Plus o Aspen
químicos mediante un software llamado simuladores de procesos. HYSYS usaban la versión anterior de AspenOne V7, que tenía gráficos bastante
diferentes a las versiones disponibles actualmente. En este libro, se utilizó
Básicamente, los simuladores de procesos pueden funcionar en dos modos: principalmente AspenOne V9; sin embargo, en algunos ejemplos también se
modo modular secuencial y modo orientado a ecuaciones (consulte la Sección aplicaron las versiones 8.6 y 8.8. Tenga en cuenta que cada software de
1.7.1); sin embargo, funcionan predominantemente en modo modular secuencial, simulación de procesos se ha desarrollado sobre la base de los mismos principios
donde los flujos de salida de un modelo de unidad se evalúan a partir de los de ingeniería química. Si uno gestiona el trabajo con uno de ellos, es fácil
flujos de entrada y los parámetros de diseño deseados. Los modelos de unidades aprender a trabajar con los demás.
Las tareas del nivel de topología del diagrama de flujo son Tanto Aspen HYSYS como Aspen Plus a través de Aspen Dynamic ics
también permiten la simulación dinámica de procesos. Este libro, sin embargo,
secuenciación de los módulos de la trata sólo de procesos en estado estacionario.
unidad, inicialización del diagrama de flujo,
El nivel de modo de propiedad física se ocupa de los flujos, la entrada del módulo cambia y, por lo tanto, el módulo debe
reevaluarse. Mas complejo
cálculo de modelos termodinámicos para equilibrio de fase,
Machine Translated by Google
los sistemas pueden contener bucles de reciclaje en diferentes niveles. número de bucles de reciclaje mucho más rápido que el enfoque
Para cada ciclo de reciclaje, se debe definir un mecanismo de iteración. modular. Sin embargo, se requiere un esfuerzo considerable para
La solución final de todo el sistema se obtiene después de la inicializar una simulación orientada a ecuaciones. Además, la
convergencia de todos los bucles de iteración. Para una convergencia construcción y depuración de una simulación orientada a ecuaciones
exitosa, la selección de un buen flujo lagrimal es muy importante. La requiere más esfuerzo.
ventaja del enfoque modular secuencial es que un problema grande Aspen Plus permite enfoques de modelado secuencial modular y
se descompone en una serie de problemas pequeños, lo que hace que orientado a ecuaciones. Sin embargo, en este libro los ejemplos se
la inicialización de la simulación de procesos sea más sencilla y fácil resuelven usando solo el enfoque modular secuencial.
de usar. Sin embargo, el requisito de buenas corrientes lagrimales
para la convergencia de bucles de reciclado es una desventaja que
hace que el enfoque modular secuencial no sea adecuado para
1.7.2 Iniciar una simulación con Aspen Plus
procesos con un gran número de bucles de reciclado.
A continuación se presenta el inicio de los simuladores Aspen Plus y
En el enfoque orientado a ecuaciones, se resuelve simultáneamente Aspen HYSYS cuando AspenOne V9 está instalado en su
un gran conjunto de ecuaciones que representa el modelo matemático computadora. Para versiones anteriores (AspenOne V8) y versiones
de todo el sistema. El enfoque orientado a la ecuación tiende a hacer posteriores (AspenOne V10), se puede hacer de forma similar.
converger procesos con una gran
Machine Translated by Google
Para iniciar una simulación de estado estable en Aspen Plus, realice Ahora puede iniciar un nuevo caso seleccionando el icono Nuevo o
los siguientes pasos: abrir un caso existente seleccionando el icono Abrir.
En la ventana del entorno de propiedades que se muestra en la
abra Aspen Plus eligiendo Figura 1.10, seleccione el icono Agregar para agregar una nueva lista
de componentes y comenzar su simulación.
1.7.3 Iniciar una simulación con Aspen HYSYS no es significativo y su presencia puede despreciarse. Sin embargo, hay
muchos procesos en los que es imposible definir la composición exacta
Abra Aspen HYSYS seleccionando de los flujos de material basados en componentes convencionales.
Tanto Aspen Plus como Aspen HYSYS contienen enormes como ensayos, mezclas, pseudocomponentes, sólidos no
bibliotecas de productos químicos convencionales y bases de convencionales, polímeros y segmentos, etc.
datos de sus propiedades. En el caso de componentes no Este libro trata no solo de componentes convencionales sino
convencionales, los simuladores utilizan enfoques específicos. también de procesos con componentes no convencionales. Los
capítulos de las Partes II y III están dedicados a
Los componentes no convencionales se clasifican en diferentes categorías
Machine Translated by Google
el diseño y simulación de procesos convencionales. En la Parte IV Los principios básicos del análisis de puntos de pellizco y su
se estudian los procesos con componentes no convencionales. aplicación mediante el análisis de energía de Aspen (AEA) se
analizan en el capítulo 12 de este libro. AEA es un software
separado de Aspen, pero también está integrado tanto en Aspen
Plus como en Aspen HYSYS; y la simulación en estos programas
1.9 Integración de Procesos y Análisis se puede exportar a AEA para el análisis de puntos de pellizco y
Energético el diseño de la red de intercambiadores de calor.
Referencias
1 Banks J, Carson J, Nelson B, Nicol D. Simulación de sistemas 5 Zhu X. Energía y optimización de procesos para las industrias de
de eventos discretos. Nueva York: Prentice Hall; 2001. procesos. Hoboken, Nueva Jersey: John Wiley & Sons, Inc.; 2014.
Disponible en http://onlinelibrary.wiley.com/book/
2 Thome B. Principio y práctica de la ingeniería de sistemas 10.1002/9781118782507.
basados en computadora, serie de Wiley en sistemas basados 6 Towler G, Sinnott R. Diseño de ingeniería química,
en software. Chichester, Reino Unido: John Wiley & Sons, Principios, práctica y economía del diseño de plantas y procesos,
Ltd.; 1992. 2ª ed. Nueva York: Elsevier, 2013.
3 Biegler LT, Grossmann IF, Westerberg W. Sistemática 7 Peters M, Timmerhaus K, West R. Diseño y economía de
Métodos de Diseño de Procesos Químicos. Nueva York: Prentice plantas para ingenieros químicos, 5.ª ed. Nueva York: McGraw
Hall PTR; 1997. Hill; 2004.
4 Demian AC. Diseño Integrado y Simulación de Procesos 8 Turton R, Bailie RC, Whiting WB, Sheiwitz JA. Análisis, Síntesis y
Químicos. Ámsterdam, Países Bajos: Elsevier; 2008. Diseño de Procesos Químicos. Nueva York: PrenticeHall, PTR;
1998.
Machine Translated by Google
15
2.1 Selección de componentes desea cambiar el nombre del componente o eliminarlo y reemplazarlo como
se muestra en la Figura 2.3. Selecciona renombrar.
Compilar una lista de componentes es una de las primeras operaciones en las Compruebe el nombre del componente y la fórmula química.
simulaciones de Aspen Plus y Aspen HYSYS. Una lista de componentes Tenga en cuenta que en los siguientes pasos de la simulación, el componente
puede constar de componentes convencionales puros y también de componentes se identifica solo con el ID del componente.
no convencionales, como análisis de petróleo, componentes hipotéticos, sólidos Utilice la selección de revisión para ver las propiedades escalares básicas
convencionales (sólidos con fórmula química conocida), sólidos no de los componentes puros en la lista.
convencionales, etc. En esta sección del libro, se presenta el método de La tabla de propiedades que se muestra en la Figura 2.4 se puede copiar a
creación de listas de componentes tanto en Aspen Plus como en Aspen HYSYS. una hoja de Excel y usarse para cualquier otro propósito.
análisis de propiedades más detallado. clic en el nombre de los parámetros y mantener presionado el mouse,
aparece el nombre exacto del parámetro, como se muestra en la Figura 2.4.
Los ejemplos a continuación tienen como objetivo mostrar la operación En la tabla 2.1 se presentan la masa molar, el punto de ebullición, la
básica de crear una lista de tipos de componentes convencionales y algunos temperatura crítica, la presión crítica y la entalpía estándar de formación
no convencionales. Algunos otros tipos de componentes no convencionales, del gas ideal a 25 ◦C para todos los componentes.
Solución: Solución:
Abra Aspen Plus siguiendo los pasos descritos en la Sección 1.7.2. Abra Aspen HYSYS siguiendo los pasos de la Sección 1.7.3.
En la tabla de selección de componentes que se muestra en la Figura 2.1, Haga clic en Agregar para agregar una lista de componentes, aparece la página
escriba "etanol" en el ID del componente y haga clic en Intro. de búsqueda de componentes.
En la siguiente línea escriba, "agua" o "H2O" y haga clic en entrar. Escriba el nombre, eventualmente una parte del nombre o la fórmula del
Para nombres más largos, como “acetato de etilo” y “ácido acético”, use la componente en la ubicación de búsqueda, HYSYS encuentra
herramienta de búsqueda o escriba el nombre debajo del Nombre del automáticamente el componente que está buscando.
componente. Utilice el icono (Buscar) para ingresar al entorno de búsqueda
de componentes. Busque acetato de etilo y ácido acético como se muestra Para agregar el componente encontrado a la lista de componentes, haga
en la Figura 2.2 y agréguelos a la lista de componentes. Para una mejor clic en Agregar (consulte la Figura 2.5).
identificación de los componentes en la lista de componentes, la Para los coeficientes de la ecuación de Antoine modificada, haga doble clic
identificación del componente se puede cambiar simplemente reescribiéndola; en el nombre del componente, aparece la tabla que se muestra en la Figura
el software le pregunta si usted 2.6. En esta tabla, puede consultar diferentes
Diseño y simulación de procesos químicos: aplicaciones Aspen Plus y Aspen HYSYS, primera edición. Juma Haydary. © 2019
John Wiley & Sons, Inc. Publicado en 2019 por John Wiley & Sons, Inc.
Sitio web complementario: www.wiley.com/go/Haydary/ChemDesignSimulation Aspen
Machine Translated by Google
Solución:
Tabla 2.1 Algunas propiedades de los componentes del proceso de acetato de etilo
megavatios
kg kmol−1 46.069 18.015 88.106 60.053
seleccione Crear un lote de hipoglucemiantes y especifique el punto de y no todos se pueden encontrar en estas bases de datos. Si conocemos la
ebullición inicial (25 ◦C), el punto de ebullición final (700 ◦C) y el intervalo fórmula estructural de un componente, podemos modelarlo como un
de temperatura entre los componentes de hipoglucemiantes individuales. componente convencional en Aspen Plus y estimar sus parámetros en función
Para crear 30 componentes hipo, se requiere un intervalo de temperatura de su estructura química.
de 25 ◦C; haga clic en Generar hipos para crear el conjunto de componentes
hipo.
Abra Aspen Plus siguiendo los pasos descritos en la Sección 1.7.2. En la
Haga clic en Agregar todo para agregar el conjunto de hipocomponentes tabla de
a la lista de componentes (consulte la Figura 2.8).
selección de componentes que se muestra en la Figura 2.1, elija un ID
para dibenzo(ah)antraceno, por ejemplo,
HAP1
Ejemplo 2.4 El dibenzo(ah)antraceno es un compuesto de hidrocarburo
Desde el panel de navegación principal, use Estructura molecular, luego
aromático policíclico (PAH) con la fórmula estructural que se muestra en la
Editar o haga doble clic en PAH1 (Figura 2.10).
figura 2.9. Modele este compuesto en Aspen Plus y estime sus propiedades
con base en su fórmula estructural.
En la página Estructura molecular, use Estructura y Dibujar/Importar/Editar,
aparece la herramienta de dibujo de la estructura molecular Usando
Solución: Las átomos, enlaces
bases de datos de Aspen Plus contienen una gran cantidad de componentes. y fragmentos, dibuje la estructura molecular como se muestra en la Figura
Sin embargo, se conocen millones de productos químicos diferentes. 2.11
Solución:
cálculo de propiedades físicas de componentes puros y corrientes y módulo de termodinámica IKCAPE deriva de las palabras alemanas
también para la determinación del equilibrio de fase. Un método de para consorcio industrial (Industrie Konsortium) y CAPE (ingeniería de
propiedad en Aspen Plus o Fluid Package en Aspen HYSYS consiste procesos asistida por computadora de cooperación industrial).
en un modelo de equilibrio de fase y diferentes modelos para cálculos
de propiedades físicas. El usuario debe seleccionar un método de Los miembros de IKCAPE son las empresas BASF, Bayer, Hoechst,
¨
propiedad o un paquete de fluidos que proporcione una representación DegussaHuls y Dow. El módulo de termodinámica IKCAPE se creó
suficientemente precisa del sistema. La selección de un método con el objetivo de diseñar un paquete de programas eficiente y completo
adecuado requiere un buen conocimiento de la termodinámica del para cálculos termodinámicos.
sistema y experiencia. A menudo, los diferentes modelos deben
compararse con los datos medidos para seleccionar el modelo más
preciso. Los programas de simulación de procesos permiten ajustar los 2.2.1.4 DECHEMA DETHERM
parámetros del modelo. En algunos casos, puede ser necesario ajustar
DECHEMA Sociedad de Ingeniería Química y Biotecnología fundada
algunos parámetros del modelo para lograr una mejor descripción de
en 1926 es una organización sin ánimo de lucro con sede en Frankfurt.
los datos medidos. Esta organización proporciona la base de datos DETHERM, que
contiene datos de propiedades termofísicas de unos 41 500
compuestos puros y 135 000 mezclas. DETHERM contiene valores de
literatura, junto con información bibliográfica, descriptores y resúmenes.
2.2.1 Fuentes de datos de propiedades físicas
En ese momento, se incluyen 8,54 millones de conjuntos de datos. La
Las propiedades físicas de miles de componentes puros se presentan base de datos se actualiza anualmente y crece continuamente con
en forma de gráficos y tablas en muchos libros de texto, manuales y alrededor del 8% anual. La base de datos DECHEMA recopila datos de
bancos de datos sobre ingeniería química y campos relacionados. Los equilibrio vaporlíquido y líquidolíquido medidos experimentalmente
resultados del trabajo de investigación sobre propiedades físicas se para miles de sistemas binarios y muchos ternarios.
publican en varias revistas de ingeniería. El Journal of Chemical
Engineering Data se especializa en la publicación de datos de
propiedades físicas para el diseño de ingeniería química. Sin embargo,
los bancos de datos de propiedades físicas computarizados son la 2.2.1.5 NIST El
mejor fuente de datos de propiedades físicas. Estos bancos de datos
Instituto Nacional de Estándares y Tecnología (NIST) es un laboratorio
pueden incorporarse al software de simulación de procesos y
de estándares de medición del Departamento de Comercio de EE. UU.
proporcionar datos de propiedades físicas evaluados para el diseño
El sistema de archivo de datos SOURCE del NIST Thermodynamic
de procesos químicos.
Research Center (TRC) actualmente contiene más de 3 millones de
A continuación, se describen algunos de los bancos de datos de
puntos de datos experimentales. Thermo Data Engine como una
propiedades físicas más grandes integrados en Aspen Plus y Aspen HYSYS:
instalación de almacenamiento integral para datos de propiedades
termofísicas y termoquímicas experimentales. La base de datos
2.2.1.1 DIPPR
experimental contiene datos de propiedades sin procesar para una gran
Design Institute for Physical Properties (DIPPR) del Instituto cantidad de componentes (más de 17,000 compuestos).
Estadounidense de Ingenieros Químicos se creó en 1978. Actualmente,
la base de datos DIPPR® 801 contiene recomendaciones para 34 Aspen Plus proporciona una herramienta de correlación, evaluación
propiedades constantes y 15 propiedades dependientes de la y predicción de datos termodinámicos llamada ThermoData Engine
temperatura para más de 2200 compuestos. (TDE). La interfaz Aspen PlusTDE contiene propiedades de valor único
de compuestos moleculares puros, como el punto de ebullición normal,
la presión crítica, la temperatura crítica, la temperatura del punto triple,
2.2.1.2 PPDS El la entalpía de formación y la energía de formación libre de Gibbs. Para
servicio de datos de propiedades físicas (PPDS) proporciona fácil los nuevos compuestos, las propiedades se estiman en función de la
acceso a las propiedades físicas, termodinámicas y de transporte y al estructura molecular utilizando diferentes correlaciones. Para
equilibrio de fase de componentes puros y
Machine Translated by Google
PURO32 Datos de muchas fuentes, incluidas DIPPR®, ASPEN, PCD, API y Banco de datos primario de componentes
Tecnología Aspen puros en los datos del sistema de
NISTTRC Datos del Instituto Nacional de Estándares y Tecnología (NIST), Programa de propiedades físicas de Aspen para una amplia gama de componentes.
Datos de Referencia Estándar (SRDP). Disponible solo con la base de datos
empresarial de propiedades de Aspen.
ACUOSO Parámetros de componentes puros para especies iónicas y moleculares en Cálculos que contienen electrolito
solución acuosa
ASPENPCD Banco de datos proporcionado con Aspen Plus 8.56 Para compatibilidad ascendente
BIODIESEL Parámetros de componentes puros para componentes que se encuentran Procesos de biodiésel
normalmente en los procesos de
COMO producción de biodiésel Parámetros de componentes puros para componentes Cálculos de fase gaseosa a alta
que se encuentran normalmente en productos de temperatura
ELECPURE combustión, incluidos los radicales libres Parámetros de componentes puros para Procesos de amina
algunos componentes que se
ETILENO encuentran comúnmente en los procesos de amina Parámetros de componentes Procesos de etileno
puros para los componentes que se encuentran típicamente en
FACTPCD los procesos de etileno para el método de propiedad SRK Especies FACT Procesos pirometalúrgicos
(componentes a los que se hace referencia en una fase pura o de solución
específica para usar solo con la interfaz Aspen/FACT/Chemapp en Aspen Plus)
HYSYS Parámetros binarios y de componentes puros necesarios para los métodos de Modelos que utilizan métodos de
propiedad de Aspen propiedad Aspen HYSYS
Iniciado HYSYS Parámetros de propiedad y parámetros de velocidad de reacción de Iniciadores poliméricos
descomposición térmica para especies de iniciadores de polímeros.
Disponible en Aspen Polymers y Aspen Properties.
INORGÁNICO Propiedades termoquímicas para componentes inorgánicos en estado de vapor, Aplicaciones de sólidos, electrolitos y
líquido y sólido Parámetro metalurgia
NRTLSAC de componente puro XYZE que contiene representaciones de Cálculos utilizando el método de
segmento de solventes comunes propiedad NRTLSAC
PCSAFT, Propiedades puras y binarias para métodos de propiedad basados en PCSAFT. Hidrocarburos cortos y moléculas pequeñas
POLIPCSF comunes
Tabla 2.3 Submodelos disponibles en Aspen Plus para la determinación de 2.2.2 Modelos de equilibrio de fase
propiedades dependientes de la temperatura de componentes puros (16)
El cálculo del equilibrio de fase es la operación clave realizada
Submodelos disponibles
por un método de cálculo termodinámico seleccionado. La
Propiedad
condición básica de equilibrio de fase es la igualdad de fugacidad
Volumen sólido ASPEN, DIPPR, IKCAPE, NIST de cada componente en cada fase. Si f es la fugacidad del
1
componente i en la fase líquida y
Volumen de líquido ASPEN (Rackett), DIPPR, PPDS,
IKCAPE, NIST f es la fugacidad de este componente en la fase de vapor,
yo G
i
Presión de ASPEN (Antoine ampliado), Wagner,
la siguiente relación es válida en el estado de equilibrio:
vapor líquido BARIN, PPDS, PML, IKCAPE, NIST
l fi = f gramo
(2.1)
Calor de i
ASPEN (Watson), DIPPR, PPDS,
vaporización IKCAPE, NIST En general, se conocen dos métodos para representar las
Capacidad calorífica sólida ASPEN, DIPPR, BARIN, IKCAPE, fugacidades en la relación (2.1): el método de la ecuación de
NIST estado y el método del coeficiente de actividad.
Capacidad de DIPPR, PPDS, BARIN, IKCAPE, NIST En el método de la ecuación de estado:
calor líquido
F gramo
= gramo
yip (2.2)
ASPEN, DIPPR, BARIN, PPDS, i i
Capacidad ideal de
calor de gas IKCAPE, NIST (2.3)
Fi = yo yo
i
chip
Coeficiente del DIPPR
segundo virial
(2.4)
eso
Conductividad ASPEN (Sato–Riedel), DIPPR, PPDS, la fase líquida, yi es la fracción molar del componente i en la fase
térmica líquida IKCAPE, NIST
de vapor, xi es la fracción molar del componente i en la fase
Conductividad ASPEN (Stiel–Thodos), DIPPR, PPDS, líquida y P es la presión total.
térmica de vapor IKCAPE, NIST
El coeficiente de fugacidad i se obtiene a partir de la
superficie liquida ASPEN (HakimSteinbergStiel), ecuación de estado, representada por P en la ecuación (2.5).
tensión DIPPR, PPDS, IKCAPE, NIST
1 ∞
= − RT
Fuente: Aspen Plus R. www.aspentech.com i
RT ∫
V [( P ni ) T,V,nj≠i V ] dV−ln Zm (2.5)
l
donde i es el coeficiente de actividad del componente i y f Este tipo de ecuación de estado se basa en la selección de las
i
es la fugacidad del componente puro a la temperatura de la mezcla. potencias de la expansión:
Con el método de la ecuación de estado, todas las propiedades para B C
+ + (2.9)
ambas fases pueden derivarse de la ecuación de estado. Utilizando
PAG = RT ( 1Vm v2
metro
V3
metro
)
el método del coeficiente de actividad, las propiedades de la fase de
vapor se derivan de la ecuación de estado como en el método de Una descripción detallada de todos los métodos de propiedad utilizados
ecuación de estado. Sin embargo, las propiedades del líquido se en Aspen Plus y Aspen HYSYS está más allá del marco de este libro.
determinan a partir de la suma de las propiedades de los componentes La tabla 2.5 enumera los modelos y los campos de aplicación de los
puros a las que se añade un término de mezcla o un término de exceso. métodos básicos de la ecuación de la propiedad estatal.
Para una relación de equilibrio, se sigue: Debido a las interacciones binarias mejoradas utilizadas en el
paquete de propiedades de PengRobinson en Aspen HYSYS, este
l modelo proporciona resultados ligeramente diferentes a los de la
yo if
ecuación de estado estándar de PengRobinson utilizada en Aspen
a= (2.8)
PAG
Plus. Estas diferencias son visibles en sistemas con cierta no idealidad.
soldado
americano
Peng–Rob Ecuación de estado cúbica estándar de PengRobinson para todas las Aplicaciones de procesamiento de gas, refinería y
propiedades termodinámicas, incluido el coeficiente de fugacidad petroquímica, torres de crudo y plantas de etileno.
de la mezcla de vapor y el coeficiente de fugacidad de la mezcla Mezclas no polares o ligeramente polares, por ejemplo,
líquida, excepto el método API de volumen molar líquido para el hidrocarburos y gases ligeros, como dióxido de carbono,
volumen molar líquido de pseudocomponentes y el sulfuro de hidrógeno e hidrógeno, particularmente adecuadas
modelo de Rackett para componentes reales, temperatura general de en las regiones de alta temperatura y alta presión. Este
componentes puros. modelo de propiedades método de propiedad se puede usar para mezclas químicas
dependientes para la presión de vapor, el calor de vaporización, la polares no ideales, si se usan las funciones alfa y las
capacidad calorífica, la viscosidad, la densidad, la conductividad térmica, reglas de mezcla apropiadas.
la tensión superficial, la entalpía, la entropía y la energía de
Gibbs La ecuación de estado de SoaveRedlich
SRK Kwong para el coeficiente de fugacidad de la mezcla de vapor y el Aplicaciones de procesamiento de gas, refinería y
coeficiente de fugacidad de la mezcla líquida. petroquímica, torres de crudo y plantas de etileno.
Mezclas no polares o ligeramente polares, por ejemplo,
El método PenelouxRauzy para la corrección del volumen hidrocarburos y gases ligeros, como dióxido de carbono,
molar líquido resulta en un volumen molar líquido más preciso. sulfuro de hidrógeno e hidrógeno, particularmente adecuadas
en las regiones de alta temperatura y alta presión. Este
método de propiedad se puede usar para mezclas químicas
Tabla de vapor NBS para calcular las propiedades del agua para una
polares no ideales, si se usan las funciones alfa y las
mayor precisión Reglas
reglas de mezcla apropiadas.
de mezcla de KabadiDanner para sistemas de agua e hidrocarburos
Modelo de
Cuando utilice el método SRK, seleccione
propiedades dependientes de la temperatura de componentes puros STEAMNBS como método de agua libre.
para presión de vapor, calor de vaporización, capacidad calorífica,
viscosidad, densidad, conductividad térmica, tensión superficial,
entalpía , entropía y energía de Gibbs HYSPR en Aspen Plus implementa
HYSP el paquete de propiedades PengRobinson utilizado Para aplicaciones de petróleo, gas o petroquímicas,
(PengRobinson en Aspen HYSYS. La diferencia entre el paquete de propiedades
PR estándar en Aspen Plus (modelo PengRob) y el paquete de – Deshidratación por TEG
Paquete de propiedades
– Deshidratación de TEG con aromáticos
en HYSYS) propiedades PengRobinson en Aspen HYSYS es que el paquete de
– Procesamiento de gases criogénicos
propiedades HYSYS contiene parámetros de interacción binarios
– Separación de aire
mejorados para todos los pares hidrocarburohidrocarburo de la biblioteca
– Torres de crudo Atm
(una combinación de pares de hidrocarburos ajustados). y parámetros de
– Torres de vacío
interacción generados), así como para la mayoría de los binarios de
– Sistemas de alto H2
hidrocarburos y no hidrocarburos.
– Sistemas de depósito
– Inhibición de hidratos
– Sistemas crudos
¨
ELECCIÓN LK La ecuación de estado de LeeKeslerPlocker para calcular LKPLOCK se puede utilizar para aplicaciones de refinería y
(LeeKesler todas las propiedades termodinámicas excepto el volumen molar procesamiento de gas, pero se prefiere el método de
Plocker) ¨ líquido de las mezclas, el método API para el volumen molar líquido de propiedad SRK o PENGROB, para mezclas no polares o
los pseudocomponentes y el ligeramente polares.
Modelo Rackett para componentes reales en mezclas También correlaciones integradas para estimar
Modelo de propiedades dependientes de la temperatura de componentes parámetros binarios de componentes CO, CO2 , N2 , H2 ,
puros para presión de vapor, calor de vaporización, capacidad calorífica, CH4 , alcoholes e hidrocarburos.
viscosidad, densidad, conductividad térmica, tensión superficial, Resultados razonables a todas las temperaturas y
entalpía, entropía y energía de Gibbs presiones; los resultados son menos precisos en la región
cercana al punto crítico de la mezcla
Machine Translated by Google
96
0 0.2 0.4 0.6 0.8 1 Los modelos de coeficientes de actividad presentados
x, y anteriormente no son adecuados para sistemas de electrolitos. En
xHYSPR yHYSPR xPengRob yMoneyRob
los electrolitos, además de las interacciones físicas y químicas
moléculamolécula, se producen reacciones iónicas e interacciones
Figura 2.19 Comparación del diagrama isobárico t–x,y de la mezcla binaria
nheptano/tolueno a 101,3 kPa calculado por el modelo Peng (moléculaion e ionion). El Aspen Property System ofrece dos
Robinson utilizado en Aspen Plus con los utilizados en Aspen modelos básicos para la predicción de electrolitos: electrolito
HYSYS NRTL (ENRTL) y modelo de coeficiente de actividad de Pitzer.
Por supuesto, ENRTL se puede combinar con diferentes
la fase líquida y la ecuación de estado de RedlichKwong para la ecuaciones de estado. Aspen HYSYS proporciona el paquete de
fase de vapor. El método de propiedad marcado solo como NRTL propiedades de electrolitos OLI para modelar sistemas de
utiliza la ecuación de estado del gas ideal para la fase de vapor. electrolitos.
POLYFH, POLYNRTL, POLYSAFT, POLYCSF POLYSK,
En Aspen HYSYS, cuando se selecciona un paquete de
propiedades de coeficiente de actividad, la siguiente opción es la POLYSRK y POLYUF son métodos de propiedades para
selección de la ecuación de estado. sistemas poliméricos (para obtener detalles, consulte la Tabla
El método de propiedad IDEAL en Aspen Plus y el paquete de 16.1 en el Capítulo 16).
propiedad ANTOINE en Aspen HYSYS también se pueden
considerar como una forma especial de modelos de actividad. 2.2.2.3 Otros métodos de propiedad
Estos modelos emplean la ley de Raoult (modelo de coeficiente Tanto Aspen Plus como Aspen HYSYS proporcionan algunos
de actividad ideal para la fase líquida, dondei =1, la ecuación de modelos de equilibrio de fase que estiman los valores de la
estado del gas ideal para la fase de vapor), el modelo de Rackett relación de equilibrio K sin aplicar ecuaciones de estado o modelos
para volumen molar líquido y la ley de Henry para gases no de coeficientes de actividad. Estos modelos suelen ser adecuados
condensables. para uno o más procesos específicos. Por ejemplo, el Braun K10,
que calcula los valores K utilizando correlaciones desarrolladas a
partir de los gráficos K10 (valores K a 10 psia) tanto para
1 componentes reales como para fracciones de aceite.
0.9 Para la simulación de fracciones de petróleo, también se utiliza
a menudo el modelo ChaoSeader, que calcula la fugacidad en la
0.8
fase líquida mediante correlaciones empíricas. El método de
0.7
propiedad CHAOSEA fue desarrollado para sistemas que
y 0.6 contienen hidrocarburos y gases ligeros, como dióxido de carbono
0.5
y sulfuro de hidrógeno, sin embargo, con la excepción del
hidrógeno. Si el sistema contiene hidrógeno, se debe utilizar el
0.4
método de propiedad de GRAYSON.
0.3 Para el proceso Amin (endulzamiento de gas por la solución de
0.2 agua Amin), Aspen PPUS proporciona el método de propiedades
AMIN que emplea el método KentEisenberg (19) para los valores
0.1
de K y la determinación de la entalpía. Aspen HYSY proporciona
0
0 0.2 0.4 0.6 0.8 1 el paquete de propiedades AMIN con modelos termodinámicos
X desarrollados por DB Robinson & Associates basados en el
HYSP ROBO DE DINERO método de KentEisenberg para su simulador de planta de amina
patentado AMSIM.
Figura 2.20 Comparación del diagrama isobárico x,y de la mezcla binaria
nheptano/tolueno a 101,3 kPa calculado por el modelo Peng
Para procesos con solo agua, se pueden usar modelos de
Robinson utilizado en Aspen Plus con los utilizados en Aspen HYSYS mesa STEAM, que proporcionan propiedades del vapor de agua
en diferentes condiciones.
Machine Translated by Google
2.2.3 Selección de un método de propiedad en Aspen Plus análisis (Sección 2.2.6). Como resultado, la herramienta Methods Asis
tance recomienda los métodos NRTLHOC o WilsonNTH (Figura 2.23).
Ejemplo 2.6 Seleccione un método de propiedad apropiado en Aspen Plus
para la simulación del proceso de acetato de etilo.
Figura 2.27 Parámetros de interacción binaria de la ecuación NRTL para pares de componentes en el proceso de acetato de etilo
Machine Translated by Google
Figura 2.28 Parámetros de interacción binaria de la ecuación de estado HOC para pares de componentes en el proceso de acetato de etilo
Solución:
parámetros de la ecuación HOC, seleccione la página Interacción binaria
y HCETA como se muestra en la figura 2.28. Seleccione la herramienta de Asistencia de Métodos como se muestra
en la Figura 2.21, aparece la página de Selección de Paquete de
AsistenciaPropiedad.
En el siguiente paso, elija Especificar tipo de componente y, en la
página siguiente, Sistema de hidrocarburos. En el
2.2.4 Selección de un Paquete de Propiedades en
siguiente paso, la Asistencia de métodos le pregunta si su sistema
Aspen HYSYS
contiene ensayos de petróleo o hipocomponentes, seleccione No.
Ejemplo 2.7 Una mezcla de alcanos C5–C8 debe procesarse en una
columna de destilación. Seleccione un paquete de propiedades HYSYS recomienda utilizar un paquete de propiedades basado en
apropiado en Aspen HYSYS para simular este proceso. ecuaciones de estado, como las ecuaciones de estado de Peng
Robinson, LKP o SRK (consulte la figura 2.29).
Para seleccionar el paquete de propiedades de PengRobinson 2.2.5 Análisis de propiedades de componentes puros
en la hoja de paquetes de fluidos, haga clic en el icono Agregar y
El sistema de propiedades físicas de Aspen permite analizar las
luego elija PengRobinson de la lista de métodos, como se
propiedades de los componentes puros y su dependencia de la
muestra en la Figura 2.30.
temperatura. El análisis de componentes puros junto con el análisis
– Para verificar la ecuación de los parámetros de interacción de
binario, los diagramas ternarios, las curvas de residuos, el análisis
estado, elija la página Coefs binarios (ver Figura 2.31).
de mezclas y la envolvente PT crean una herramienta de análisis en
Aspen Plus. La herramienta de análisis de componentes puros
Se puede extraer información breve sobre el paquete de permite calcular la dependencia de la temperatura de las diferentes
propiedades PengRobinson de la ayuda de Aspen HYSYS (consulte propiedades termodinámicas y de transporte de los componentes
la Figura 2.32). puros.
50,000 2.00E–03
45.000
1.50E–03
40.000 Viscosidad
dinámica
(Pa∙s)
Etanol 1.00E–03
35,000
H2O
vaporización
(kJ∙kmol1)
Calor
de
5.00E–04
30.000 Acetato
Ácido
25,000 0.00E+00
0 20 40 60 80 100
20,000 Temperatura (oC)
0 50 100 150 Agua Etanol Acetato de etilo Ácido acético
Temperatura (°C)
Figura 2.36 Dependencia de la temperatura de la viscosidad dinámica para los
Figura 2.34 Dependencia de la temperatura del calor molar de vaporización componentes del proceso de acetato de etilo
de los componentes del proceso de acetato de etilo
Para ver los resultados en forma de tabla, haga clic en Resultados como
Ejemplo 2.8 Calcular la dependencia de la temperatura de la entalpía de
se muestra en la Figura 2.35.
vaporización en el rango de temperatura de 20–150 ◦C y la viscosidad
Los datos se pueden copiar a una hoja de Excel para su posterior
dinámica en el rango de temperatura de 0–100 ◦C de todos los componentes
análisis.
en el proceso de acetato de etilo en el rango de temperatura de 20–150 ◦C.
Para calcular la dependencia de la temperatura de la viscosidad
dinámica, seleccione Transporte como tipo de propiedad y MU para
Solución: viscosidad.
Los resultados se muestran en la Figura 2.36.
Elija Puro de la herramienta Análisis como se muestra en el paso 1 en
la Figura 2.33.
2.2.6 Análisis binario
Como tipo de propiedad, seleccione Termodinámica, luego DHVL y una
unidad, por ejemplo, kJ kmol−1 . Los resultados de los modelos de equilibrio de fase descritos en la Sección
Seleccione el rango de temperatura y el número de puntos 2.2.3 generalmente se presentan en forma de dependencias entre P, T, xi ,
o incrementos. yi , ki , y los modelos de equilibrio yo _ En muchos casos, las ecuaciones de
De la lista de componentes disponibles, seleccione los componentes de fase contienen parámetros estimados mediante el ajuste de datos
para los que debe calcularse el calor de vaporización. experimentales medidos para sistemas binarios. La no idealidad de una
mezcla puede presentarse mediante parámetros de interacción binaria de
Haga clic en Ejecutar análisis para calcular la propiedad seleccionada. diferentes pares de componentes en la mezcla. En cálculos de diseño y
El gráfico de dependencia de la temperatura de la propiedad seleccionada simulación de muchas unidades de separación, cálculos de etapa de
aparece automáticamente (ver Figura 2.34). equilibrio
Figura 2.35 Visualización de los resultados del análisis de propiedades de componentes puros
Machine Translated by Google
son usados. La capacidad de separación de una etapa de equilibrio en la simulación del proceso que contiene unidades de separación debe
el caso de equilibrio vaporlíquido viene dada por la volatilidad relativa incluir este paso.
definida como El banco de datos NIST integrado en el Aspen Property System
contiene una gran cantidad de conjuntos de datos de equilibrio binario
= A (2.10)
yo
medidos experimentalmente. Otra gran fuente de datos de equilibrio
kj binario es el banco de datos DECHEMA.
En caso de equilibrio líquidolíquido, el factor de separación es la La herramienta de análisis binario Aspen Plus permite calcular los datos
selectividad relativa definida como de equilibrio isobárico T–xy e isotérmico P–xy utilizando diferentes
modelos de equilibrio de fase disponibles en Aspen y comparar los
= Salvado
(2.11) resultados así obtenidos con los datos medidos experimentalmente.
yo
KLj
77.0
76.5
76.0
75,5
Temperatura,
°C
75,0
74.5
74.0
73.5
73.0
72.5
72.0
71,5
0,00 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50 0,55 0,60 0,65 0,70 0,75 0,80 0,85 0,90 0,95 1,00
Fracción de masa de líquido/vapor, Eacetato
Figura 2.38 Diagrama Txy de la mezcla de acetato de etilo y etanol generado por la herramienta de análisis binario de Aspen
Después de ejecutar el análisis, aparece el diagrama Txy para los Aspen permite mostrar diferentes tipos de gráficos.
sistemas binarios de acetato de etiloetanol que se muestra en la Para mostrar cualquier otro tipo de gráfico, por ejemplo, xy, Ki =f(xi ),
figura 2.38. i =f(xi ), use la barra de herramientas Plot como se muestra en la
Para ver los resultados completos en formato de tabla, haga clic en figura 2.40.
Resultados debajo de su Análisis binario como se muestra en la Fig. Transporte los resultados a una hoja de Excel para un análisis más
ura 2.39. detallado o use Merge Plot en Aspen para combinar
Figura 2.39 Visualización de datos completos de equilibrio isobárico calculados por la herramienta Análisis binario
Machine Translated by Google
diferentes parcelas y comparar los resultados de diferentes Los datos experimentales también se pueden transportar a la
modelos. misma hoja de Excel o comparar directamente en Aspen Plus con
Calcule los datos de equilibrio isobárico para el sistema binario los datos del modelo. En la Figura 2.43 se presenta una
acetato de etiloetanol utilizando el modelo NRTL, que a comparación de los datos VLE calculados por NRTL y NRTLHOC
diferencia del NRTLHOC considera la fase de vapor como gas con los datos VLE medidos experimentalmente. En el caso del
ideal. Para calcular los datos de equilibrio utilizando el modelo sistema binario acetato de etiloetanol, los datos calculados por
NRTL, simplemente cambie el método de propiedad en la página NRTL son muy similares a los calculados por NRTL HOC. Los
de entrada de datos de análisis binario (Paso 7 en la figura datos de ambos modelos tienen muy buena coherencia con los
2.37). medidos experimentalmente. Con base en las mediciones y
Transporte los datos calculados a la misma hoja de Excel que modelos experimentales, los datos de esta mezcla binaria crean un
los de NRTLHOC para comparar los resultados de ambos azeótropo con un punto de ebullición de 72 ◦C y una fracción molar
modelos. Los resultados también se pueden comparar de 0,55 de acetato de etilo y 0,45 de etanol.
directamente en Aspen Plus usando Merge Plot para combinar
diferentes gráficos (consulte la Figura 2.40). Para otros sistemas binarios, se puede utilizar el mismo
procedimiento que para el sistema acetato de etiloetanol.
Los datos experimentales para el sistema binario acetato de La figura 2.44 muestra una comparación de datos VLE isobáricos
etiloetanol se pueden extraer de la base de datos del NIST. Para a 101,325 kPa para el sistema binario etanolácido acético. Estos
usar el motor NIST ThermoData, siga los pasos que se muestran diagramas Txy se produjeron y compararon directamente en
en la Figura 2.41. Aspen Plus; Se pueden ver diferencias entre los datos VLE de
El motor NIST ThermoData proporciona una serie de datos etanolácido acético calculados por NRTL y NRTLHOC porque el
diferentes para el sistema binario acetato de etiloetanol, incluidos modelo NRTL no considera la fase de vapor.
varios conjuntos de datos de equilibrio líquidovapor isobárico
(VLE) (Figura 2.42). Algunos de estos datos se midieron a 101,325 Datos experimentales de VLE isobárico medidos a 101 325 Pa
kPa, que es la presión aplicada en nuestro ejemplo, por ejemplo, publicados por Amezaga et al. (21) se puede extraer del motor
datos publicados por Li et al. (20) se puede utilizar. Si es necesario, NIST ThermoData. De la Figura 2.45 resulta que el modelo NRTL
se pueden comparar diferentes conjuntos de datos experimentales. HOC proporciona una mejor descripción de los datos
experimentales.
Machine Translated by Google
en comparación con el modelo NRTL, aunque este modelo azeótropo con punto de ebullición mínimo en torno a los 98 ◦C y
también muestra una desviación, especialmente en los datos T=f(x). fracción molar de agua de 0,85. Sin embargo, este azeótropo no
El último sistema binario analizado en este ejemplo es el fue registrado por experimentos proporcionados por Chang et al.
sistema binario aguaácido acético. Como resulta de la figura (22). Este es un hallazgo muy importante, porque en muchos
2.46, el efecto de asociación de la fase de vapor es más visible cálculos de separación de agua y ácido acético, generalmente se
en este sistema binario. Los datos VLE calculados por el modelo usa el modelo NRTL.
NRTL considerando la fase de vapor ideal muestran desviaciones La descripción de datos VLE por el modelo NRTLHOC es
significativas de los datos experimentales y los datos calculados buena y no muestra ningún azeótropo para esta mezcla binaria.
por NRTLHOC. El modelo NRTL muestra una Como conclusión del ejemplo 2.9, se puede decir
79 79
78 78
77 77
76 76
Temperatura
(°C)
75 Temperatura
(°C) 75
74 74
73 73
72 72
71 71
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
Figura 2.43 Comparación de datos VLE para el sistema binario acetato de etiloetanol calculado por NRTL (b), NRTLHOC (a) y datos experimentales
105
100
Temperatura,
C
95
90
85
80
75
0,00 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50 0,55 0,60 0,65 0,70 0,75 0,80 0,85 0,90 0,95 1,00
Figura 2.44 Comparación de datos VLE de etanolácido acético calculados por NRTL y NRTLHOC directamente en Aspen Plus
125 120
120
115 115
110
105 110
100
Temperatura
(°C)
Temperatura
(°C)
95
105
90
85
100
80
75
0 0.2 0.4 0.6 0.8 1
95
fracción molar, acetato de etilo 0 0.2 0.4 0.6 0.8 1
xNRTL yNRTL xNRTLHOC yNRTLHOC xNIST yNIST Fracción molar, agua x
xNRTL yNRTL NRTLHOC yNRTLHOC xNIST yNIST
que NRTLHOC es un modelo apropiado para la simulación del proceso mezcla ternaria. A partir de la síntesis del mapa ternario, proponer un split
de producción de acetato de etilo. de columnas de destilación para la separación del acetato de etilo.
para mezclas ternarias y determinar la viabilidad de splits para Para ver información sobre los azeótropos, haga clic en Azeotropos
como se muestra en el paso 4 de la Figura 2.48.
columnas de destilación.
La síntesis de destilación de Aspen encontró cuatro azeótropos
diferentes, tres homogéneos y uno heterogéneo: azeótropo binario
Además, los mapas ternarios y las curvas de residuos de Aspen Plus se homogéneo de etanol y agua, azeótropo binario homogéneo de acetato
pueden usar para trazar gráficos ternarios y mapas de curvas de residuos. de etilo y etanol, azeótropo binario heterogéneo de acetato de etilo y agua
y, finalmente, azeótropo ternario homogéneo de agua. acetato de etilo y
Ejemplo 2.10 Usando la síntesis por destilación de Aspen, encuentre etanol. Excepto por un azeótropo ternario, que es un nodo inestable,
todos los azeótropos que forman una mezcla de acetato de etilo, etanol y todos los azeótropos binarios son sillas de montar. La Tabla 2.6 muestra
agua a 101.325 kPa usando el modelo NRTLHOC. información detallada sobre todos los azeótropos. La lista completa de
Especificar el tipo de azeótropos y puntos singulares. Dibuje límites de
destilación y mapas de curvas de residuos para este
los puntos singulares se pueden visualizar seleccionando Puntos triángulo (Figura 2.49). Es muy fácil cambiar la forma del diagrama
singulares. Los resultados se muestran en la Tabla 2.7. usando el ícono apropiado en el lado derecho del diagrama mostrado.
Elija el gráfico Ternario en la página de entrada (Figura 2.48) La síntesis de destilación de Aspen proporciona datos VLE y LLE,
para mostrar el diagrama ternario del sistema acetato de etiloetanol y el gráfico de la curva LLE que especifica el área donde se forman
agua. Un diagrama ternario se puede representar en forma de dos fases líquidas. Posición de los azeótropos y destilación
triángulo rectángulo o equilátero.
Temperatura Número de
( ◦C) Clasificación Tipo componentes EACETATO ETANOL AGUA
Temperatura Número de
( ◦C) Clasificación Tipo componentes EACETATO ETANOL AGUA
78,15 ºC 78,15 ºC
0,9 0,1 0.9
0.6
2
ETANOL 0.4
2
0.6
(100,02 °C) 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 (77,20 °C) 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
EACETATO
también se muestran los límites. Como se muestra en la Figura concentración del área 1. Cuando dicha mezcla se alimenta a una
2.49, los límites de destilación dividen el diagrama en tres áreas columna de destilación basada en el mapa de curvas de residuos
de destilación. Dependiendo de la concentración inicial de como un producto destilado, se obtiene un azeótropo ternario de
alimentación, se pueden esperar diferentes productos de destilación. acetato de etiloetanolagua o su mezcla con el azeótropo acetato
Se puede agregar una curva de residuo (o curva de destilación) etanol, y un se obtiene una mezcla de etanol y agua como producto
al diagrama usando los íconos Agregar una curva o Agregar curva de fondo. Al agregar una cantidad específica de agua al producto
por valor en el lado derecho del diagrama que se muestra. Puede destilado, se obtiene una mezcla ternaria heterogénea, que se
agregar cualquier número de curvas y líneas de unión en el área puede separar en una fase de agua y una fase de acetato. La fase
heterogénea. La figura 2.50 muestra un diagrama ternario con de acetato se destila para producir acetato de etilo prácticamente
curvas de residuos adicionales. puro como producto de fondo. La composición del producto
Los diagramas ternarios con límites de destilación y mapas de destilado está cerca de la composición del azeótropo ternario y
curvas de residuos permiten determinar la configuración de las puede reciclarse al separador líquidolíquido. El
columnas de destilación. Considere una mezcla con el
78,15 ºC
78,15 ºC
0,9 0,9
0.8 0.8
0.7 0.7
0.6 0.6
0.5 0.5 1
71,78 ºC 71,78 ºC
0.4 0.4
0.3 0.3
2
0.2 0.2
70,33 ºC 70,33 ºC
4
3
0.1 0.1
Figura 2.50 Mapa de curvas de residuos del sistema acetato de etilo Figura 2.51 Vía de separación del acetato de etilo de la mezcla ternaria
etanolagua acetato de etiloetanolagua
Machine Translated by Google
Alimentar
Agua Fase
Agua+Etanol de agua
reciclado al
reactor Acetato de etilo
los pasos mencionados se muestran en los mapas de curvas de Ejemplo 2.11 Calcular las relaciones entre la temperatura y la presión
residuos en la Figura 2.51, y el esquema del proceso se propone en del punto de burbujeo/punto de rocío
−1
la Figura 2.52. de una mezcla que contiene 0,51 kmol h 0,21 de etanol,
−1
kmol z h−1 de acetato de etilo y 0,18 kmol h de ácido acético.
Además de las fracciones de vapor 0 y 1, realice el cálculo también
2.2.8 Análisis de envolvente de PT para las fracciones de vapor 0,25, 0,5 y 0,75.
Solución:
El análisis de envolvente de PT permite generar tablas y gráficos
que muestran la relación entre temperatura, presión y fracción de En la barra de herramientas de Análisis, seleccione Envolvente PT.
vapor de mezclas con diferentes composiciones. Con la configuración En la página Envolvente de PT, especifique la composición de
predeterminada, este análisis calcula la temperatura del punto de la mezcla, las fracciones de vapor adicionales y los puntos
rocío (fracción de vapor = 1) y la temperatura del punto de burbuja máximos, como se muestra en la Figura 2.53.
(fracción de vapor = 0) a diferentes presiones. El usuario puede Después de presionar Ejecutar análisis, se muestra la relación
especificar fracciones de vapor adicionales. Este análisis en Aspen gráfica que se muestra en la Figura 2.54.
Plus se limita a las fases líquidovapor; si se forman dos fases Se muestra la relación entre la temperatura y la presión de la
líquidas, se debe utilizar otro tipo de análisis. mezcla ternaria líquidovapor para temperaturas y presiones
desde 49 ◦C y 0,3 bar hasta su
Envolvente PT
80
75
70 VFRAC = 0,0
sesenta y cinco VFRAC = 1,0
60 VFRAC = 0,25
55
VFRAC = 0,75
50
VFRAC = 0,5
45
40
Presión,
barra
35
30
25
20
15
10
5
0
40 60 80 100 120 140 160 180 200 220 240 260 280
Temperatura, °C
Figura 2.54 Resultados del análisis de envolvente PT para todo el rango de temperaturas y presiones
Solución:
Otro paso clave en el diseño y la simulación de procesos químicos es la
Aspen HYSYS permite definir cinco tipos de reacciones:
definición de la química del proceso y las reacciones químicas. El
modelado de reactores químicos se presentará en el Capítulo 5. En este 1. Conversión: use este tipo de reacción si se conoce la conversión.
capítulo, solo se presenta la definición de química de procesos y Tenga en cuenta que las reacciones de conversión se pueden
reacciones químicas en un simulador de procesos. Si se usa Aspen modelar utilizando el modelo de reactor de
HYSYS para la simulación, la reacción química se define completamente conversión 2. Equilibrio: defina la reacción como una reacción de
en el entorno de propiedades antes de ingresar al entorno de simulación. equilibrio si se va a calcular la conversión de equilibrio.
Aspen Plus permite HYSYS puede calcular la constante de equilibrio Ke basada en la
energía libre de Gibbs ΔrG usando la ecuación (2.12)
Machine Translated by Google
ΔrG
ln Ke = (2.12)
RT
′
(conc.)k+ y k– son constantes
r = k+f (conc.) − k−f (2.13) donde
de velocidad de las reacciones directa e inversa, respectivamente;
definida por la ecuación de Arrhenius. Las funciones f y f ' son
productos de concentraciones:
−mi
RT T (2.14)
re k+ = UN mi
−mi′
k− = A d' RT T ' (2.15)
′
F (conc.)
(2.16)
r = k [ f (conc.) − K′ ]
Para crear reacciones, conjuntos de reacciones y agregarlos al
5. Catalítica heterogénea: use este tipo de reacción para modelar una paquete de fluidos en este ejemplo, siga los pasos a continuación:
reacción catalítica heterogénea si la ecuación de velocidad está
dada por el método de LangmuirHinshelwoodHougenWatson:
Abra Aspen HYSYS, cree una nueva lista de componentes y agregue
todos los componentes de las reacciones químicas a la lista (para
r=
k+f (conc.) − k−f ( 1 ′ (conc.) obtener detalles, consulte el Ejemplo 2.1).
(2.17)
+ K1 f1 (conc.) + K2 f2 (conc.)… )n Seleccione un paquete de fluido apropiado como se explica en el
Ejemplo 2.6.
donde K1 , K2 , … Kn representan constantes de adsorción y también Seleccione Reacciones en el menú principal y haga clic en Agregar como
están dadas por la ecuación de Arrhenius.
se muestra en la Figura 2.56; aparece la hoja Información del conjunto de
Las primeras tres reacciones en este ejemplo deben definirse como
reacciones.
el tipo de conversión y la última como el equilibrio. Debido a que se
Después de hacer clic en Agregar reacción (paso 3 en la Figura
utilizan diferentes modelos de reactores para las reacciones de
2.57), aparece el menú Tipo de reacción. El usuario puede elegir
conversión y equilibrio, se deben crear dos conjuntos de reacciones
las propiedades HYSYS seleccionadas por defecto o las
diferentes. Para definir un conjunto de reacciones en Aspen HYSYS y
propiedades Aspen para las reacciones. Deje las propiedades
agregarlo al Paquete de fluidos, se requieren los siguientes pasos:
HYSYS seleccionadas por defecto y elija el tipo de reacción de
Conversión (paso 4 en la Figura 2.57).
1. crear un conjunto de Ingrese a la hoja de reacción haciendo doble clic en el Nombre de
reacciones, 2. agregar reacciones al conjunto la reacción (paso 5 en la Figura 2.57).
de reacciones, 3. definir la estequiometría y otros parámetros de las Defina la estequiometría y la conversión de la reacción como se
reacciones, y 4. muestra en la figura 2.58; HYSYS permite ingresar la dependencia
agregar el conjunto de reacciones al paquete de fluidos (FP). de la temperatura de conversión en forma:
Machine Translated by Google
Figura 2.57 Adición de reacciones a un conjunto de reacciones y el conjunto de reacciones al paquete de fluidos
Conversión (%)=C0 +C1T+C2T Se desconoce 2 ; si la temperatura Agregue un nuevo conjunto de reacciones para la reacción de cambio de
la dependencia de la conversión, ingrese solo C0 como un valor constante agua en el equilibrio (Figura 2.56).
Utilizando los mismos pasos, defina otras dos reacciones de conversión Equilibrio.
Agregue el conjunto de reacciones a FP como se muestra en el paso 6 de Ke como se muestra en la Figura 2.59.
la figura 2.57. Agregue el conjunto de reacción a FP.
Figura 2.59 Definición de la estequiometría y método de cálculo de Ke para una reacción en equilibrio
Ejemplo 2.13 Una solución acuosa de HCl se neutraliza en un 2. Ionización de equilibrio de HCl en agua: H2O + HCl ↔ H3O+
depurador con una solución acuosa de NaOH. +Cl− 3.
Defina la química de las reacciones electrolíticas para este proceso Disociación de NaOH: NaOH ↔ Na+ + OH− 4.
en Aspen Plus. Precipitación de NaCl: NaCl → Na+ +Cl−
Solución:
Para modelar su química en Aspen Plus, abra un caso como se
Si se mezcla una solución acuosa de HCl con una solución acuosa
describe en el Capítulo 1, pero en lugar de una simulación en
de NaOH, se producen las siguientes reacciones electrolíticas: 1.
blanco, seleccione Electrolitos y Electrolitos con valores métricos.
Autoionización
unidades.
de equilibrio del agua: 2H2O ↔ H3O+ Aspen Plus agrega automáticamente agua como componente a
+ OH− la lista de componentes porque el agua es un componente obligatorio
componente en el modelado de sistemas de electrolitos. Además, ELECNRTL Aspen Plus proporciona cuatro métodos para la especificación química:
se selecciona automáticamente como método termodinámico; que puede especificando reacciones o componentes reactivos, seleccionando todos los
consultar en la página Métodos → Especificaciones. componentes como reactivos y especificando componentes inertes. En este
caso, seleccione Especificar reacciones y siga los pasos que se muestran en la
Compruebe los parámetros de interacción binaria y también los parámetros Figura 2.61.
del par de electrolitos (para más detalles, consulte el Capítulo 15). Se pueden definir tres tipos de reacciones electrolíticas: equilibrio, precipitación
Cree una nueva química siguiendo los pasos que se muestran en la Figura 2.60. de sal y reacción de disociación. Para
autoionización del agua e ionización del HCl, use el tipo de coeficientes estequiométricos de iones como se muestra en la
reacción Equi librium; para la disociación de NaOH, utilice el tipo figura 2.63. El conjunto definido de reacciones electrolíticas se
Disociación; y para la precipitación de NaCl, utilice el tipo Sal. presenta en la figura 2.64.
Fase
Agua de agua
Acetato de etilo
C2
Reactor R
CColumna de destilación
Mezclador M
SEPSeparador líquidolíquido
equipos y muestra todos los procesos y flujos de servicios, incluidos los (Ejemplo 2.6) en Entorno de Propiedades. A continuación, cambie al entorno
sistemas de control de procesos. Un PFD más detallado es el diagrama de de simulación. Aparece la Paleta de modelo. Aspen Plus proporciona una
tuberías e instrumentación (P&ID). Sin embargo, estos diagramas son serie de modelos predefinidos de operaciones unitarias comúnmente llamados
demasiado complejos para usarlos en cálculos de simulación y optimización. "Bloques" divididos en diferentes grupos como intercambiadores, reactores,
Por otro lado, no es necesaria la inclusión de todos los componentes del columnas, separadores, cambiadores de presión, etc.
diagrama en todos los tipos de cálculos.
Haga clic en Reactores y seleccione el modelo CSTR como se muestra en
En el caso de la simulación de un proceso existente, en primer lugar, en el la Figura 2.66. En nuestro PFD se incluyen tres columnas de destilación; el
esquema tecnológico del proceso, se especifican las unidades y equipos que modelo apropiado es el modelo RadFrac.
pueden afectar la simulación. El siguiente paso es la creación de un esquema Utilice el modelo RadFrac para una simulación rigurosa de todos los tipos de
de cálculo que contenga solo los modelos que siguen el objetivo del cálculo. operaciones de fraccionamiento líquidovapor de varias etapas, como columnas
En caso de diseñar un nuevo proceso, se utilizan diagramas de bloques muy de destilación, absorbedores, stripers, etc. Para seleccionar el modelo
simples en más alternativas. En el siguiente paso, se desarrollan esquemas RadFrac, haga clic en Columna, seleccione RadFrac, luego haga clic tres
de cálculo más detallados. PFD y P&ID se crean en los siguientes pasos del veces en el campo; después de cada clic, se agrega un modelo RafFrack al
desarrollo del proyecto. campo. Haga clic en la flecha en la esquina superior izquierda de la Paleta de
modelos para dejar de agregar modelos. El orden de los pasos se muestra en
la Figura 2.67.
Para completar el esquema de bloques de proceso, agregue un mezclador y
El ejemplo 2.14 Concepto tecnológico del proceso de acetato de etilo en el un modelo de decantador.
Capítulo 8 describe tres posibles procesos continuos: proceso de reactor de Un PFD consta de bloques y flujos. Aspen Plus proporciona tres tipos de
tanque agitado continuo en fase líquida, proceso de destilación reactiva y flujos: flujos de materiales, flujos de calor y flujos de trabajo. Para continuar
proceso de reactor tubular en fase de vapor. Dibuje esquemas de bloques con la simulación del proceso, se deben dibujar al menos todos los flujos de
simples del proceso de reactor de tanque agitado continuo en fase líquida y material requeridos. Al hacer clic en el Modelo de flujo de materiales, todos
prepare diagramas de flujo de Aspen Plus para este proceso. los flujos de materiales obligatorios se muestran con flechas rojas (Figura
2.68). Las flechas azules representan las corrientes de agua libre, cuyo dibujo
es necesario si hay agua libre presente en el proceso.
Solución: En
El usuario puede cambiar el nombre de todos los flujos y bloques. Ejemplo 2.15 Desarrolle un diagrama de flujo de cálculo simple en Aspen HYSYS que
Para cambiar el nombre de un bloque o flujo, haga doble clic en el nombre del bloque o se usará para calcular
de flujo del acetato de etilo preparado para las siguientes operaciones. 1. requerimiento de calor para precalentar la alimentación a la temperatura de reacción
temperatura,
gases combustibles
H2
tolueno
m ÉL R ÉL ÉL Benceno
100 101 101 102 en 101 en 102 103 101
Hidrógeno
Figura 2.70 Diagrama de flujo simple del proceso de hidrodesalquilación de tolueno para la solución del Ejemplo 2.15
Figura 2.72 Diagrama de flujo de hidrodesalquilación de tolueno basado en los requisitos del Ejemplo 2.15
Si el solucionador está en el modo Activo, HYSYS resuelve automáticamente Modelo de reactor de conversión para el cálculo del balance de materia y
el modelo después de que se conecta y se definen sus parámetros. En el entalpía del reactor, considerando que se conoce la conversión de la
modo En espera, primero se deben crear todas las conexiones y definir los reacción, también Modelo de intercambiador de calor para
parámetros de los modelos y flujos, luego puede comenzar el cálculo al calcular el requerimiento de agua de refrigeración, Modelo de separador
cambiar al modo Activo. La figura 2.71 muestra cómo se seleccionan los para simular la
modelos individuales y se definen sus conexiones. Después de seleccionar separación de fases en equilibrio a alta presión, Modelo de válvula para
el modelo de Model Palette y hacer doble clic en su nombre, aparece la página simular la reducción de presión,
Model. Los flujos de entrada y salida del modelo se seleccionan de los flujos el modelo de separador para simular la separación de fases
ya existentes o se definen directamente en esta página ingresando un nombre en equilibrio a baja presión, el modelo de calentador para simular el
para el flujo. precalentamiento de la
alimentación de la columna de destilación y el subdiagrama de flujo de la
columna de destilación
para simular la separación de benceno del tolueno.
Además del mezclador que se muestra en la Figura 2.71, los siguientes
modelos deben incluirse en este proceso:
El PFD creado en Aspen HYSYS se presenta en la Figura 2.72. También
Modelo de calentador para calcular el trabajo de calor requerido para se muestra la posición de los modelos individuales utilizados.
precalentar la alimentación a la temperatura de reacción,
Referencias
1 Frenkel M, Chirico RD, Diky V, Yan X, † Dong Q, Muzny C. Motor actualizaciones a través de la web. J. Chem. información Modelo.
ThermoData (TDE): implementación de software del concepto 2007;47(4): 1713–1754.
de evaluación dinámica de datos. J. Chem. información Modelo. 3 Joback KG, Reid RC. Estimación de las propiedades de los
2005;45(4):816–838. componentes puros a partir de las contribuciones de los grupos. química Ing.
2 Diky V, Muzny C, Lemmon EW, Chirico RD, Frenkel M. ThermoData común 1987;57(1–6):233–243.
engine (TDE): Implementación de software del 4 Constantinou L, Gani R. Nueva contribución grupal
concepto de evaluación dinámica de datos 2, Ecuaciones de Método para estimar propiedades de compuestos puros.
estado bajo demanda y dinámicas AIChE J. 1994;40(10):1697–1710.
Machine Translated by Google
5 MarreroMorejon J, PardilloFontdevila E. Estimación de propiedades 2001. Ver pág. 9.9 para gas a baja presión y p. 9.35 el modelo
de compuestos puros utilizando contribuciones de interacción de Lucas para alta presión.
grupos. AIChE J. 1999;45(8):615–621. 15 Chung TH, Lee LL, Starling KE. Aplicaciones de teorías cinéticas de
6 Wilson GM, Jasperson LV. Constantes críticas Tc , Pc . gases y correlación multiparamétrica para la predicción de
Estimación basada en métodos de cero, primero y segundo orden. la viscosidad y conductividad térmica de gases diluidos.
En Reunión AIChE, Nueva Orleans, LA; 1996. Ing. Ind. química Fundam. 1984;23(1): 8–13.
7 Ambrose D, Walton J. Presiones de vapor hasta su
temperaturas críticas de alcanos y alcanoles normales. 16 Ayuda de Aspen PlusV9; Bedford, MA: AspenTech; 2016.
aplicación pura química 1989;61(8):1395–1403. 17 Peng DY, Robinson DB. Una nueva ecuación de estado de dos
8 Yamada T, Gunn RD. Volúmenes molares de líquidos saturados. constantes. Ing. Ind. química Fundam. 1976; 15 (1): 59–
La ecuación de Rackett. J. Chem. Ing. Datos 1973;18: 234–236. 64.
18 Fredenslund AA, Jones RL, Prausnitz JM. Estimación de la
9 Riedel L. La densidad del líquido en el estado de saturación. contribución grupal de los coeficientes de actividad en mezclas
Chem.Ing.Tech. 1954; 26:259–264. Modificado en JL líquidas no ideales. AIChE J. 1975;21(6): 1086–1099.
Hales, R. Townsend. Densidades de Líquidos de 293 a 490 K
de Nueve. Hidrocarbonos aromáticos. J. Chem. 19 Kent RL, Eisenberg B. Mejores datos para el tratamiento de aminas.
Termodina. 1972: 4 (15): 763–772. Proceso de Hidrocarburos. 1976;55(2):87–92.
10 Benson SW, Buss J H. Additivity rules for the 20 Li Q, Zhang J, Lei Z, Zhu J, Xing FJ.
estimación de propiedades moleculares. Propiedades equilibrio líquidovapor para acetato de etilo + etanol +
termodinámicas. J. Chem. física 1958;29(3):546–572. tetrafluoroborato de 1etil3metilimidazolio. química
11 Bondi, A. Propiedades Físicas de Cristales Moleculares, Líquidos Cercanamente. Datos 2009;54(2):194197.
y Vidrios. Nueva York: Wiley; 1968. 21 Amezaga SA, Biarge JF. Comportamiento
12 Sastri SRS, Rao KK. Un nuevo método de contribución grupal para Termodinámico de Agunos Sistemas Binarios Ácido Acético
predecir la viscosidad de líquidos orgánicos. química Ing. J. Alcohol en Equilibrio Liquido Vapor. Ene.
1992;50(21):9–25. Quim. 1973;69:569–586.
13 Chung TH, Ajlan M, Lee MM, Starling KE. Correlación multiparamétrica 22 Chang W, Xu X, Li X, Guan G, Yao H, Shi J. Equilibrio líquidovapor
generalizada para propiedades de transporte de fluidos polares y no isobárico para agua + ácido acético + acetato de npentilo
polares. Ing. Ind. química Res. 1988;27(4): 671–679. Huagong Xuebao 2005;56(8), 6, 1379.
23 Turton R, Bailie RC, Whiting WB, Sheiwitz JA.
14 Poling BE, Prausnitz JM, O'Connell JP. Las propiedades de los gases Análisis, Síntesis y Diseño de Procesos Químicos.
y líquidos, 5ª ed. Nueva York: McGrawHill; Nueva York: Prentice Hall, PTR; 1998.
Machine Translated by Google
61
Parte II
63
Intercambiadores de calor
El diseño de intercambiadores de calor y redes de intercambiadores de calor Programas de clasificación (EDR): Shell&Tube, AirCooled o Plate; Este modelo
crea una parte considerable del diseño de ingeniería química. La clasificación de incluye una variedad de tipos de intercambiadores de calor de carcasa y tubos en
los intercambiadores de calor, su aplicación y los métodos de cálculo se dan en comparación con el método detallado que realiza el diseño completo y la
muchos libros de texto de ingeniería química o en libros de texto sobre transferencia simulación de diferentes tipos de intercambiadores de calor, incluidas todas las
de calor e intercambiadores de calor (1–3). combinaciones de transferencia de calor monofásica, de ebullición o de
condensación, con cálculos de caída de presión asociados, vibraciones mecánicas.
Tanto Aspen Plus como Aspen HYSYS permiten el diseño y la simulación de análisis de ruidos y estimación del ensuciamiento máximo.
intercambiadores de calor con diferentes niveles de complejidad. Los modelos de
funcionamiento de la unidad Aspen Plus de los intercambiadores de calor son los
MHeatX: MHeatX permite el modelado de la transferencia de calor entre múltiples
siguientes:
corrientes frías y calientes, como los intercambiadores de calor de GNL (gas
Heater: Heater se usa para modelar intercambiadores de calor si el tema de natural licuado). Asegura un balance general de energía pero no tiene en
estudio o la importancia es solo de un lado. cuenta la geometría del intercambiador.
Heater puede modelar calentadores, enfriadores u otros dispositivos como
tanques, mezcladores y válvulas. El bloque calefactor permite especificar la HxFlux: el modelo HxFlux puede realizar cálculos de transferencia de calor entre
temperatura o el trabajo térmico de la unidad, pero no realiza ecuaciones de un disipador de calor y una fuente de calor utilizando la transferencia de calor
intercambio de calor rigurosas. Se puede especificar cualquier número de por convección. La fuerza impulsora de la transferencia de calor por convección
corrientes de alimentación para el bloque calentador. se calcula en función de la diferencia de temperatura media logarítmica.
de calor y geometría definida del intercambiador de calor. 3.1 Modelos de calentadores y enfriadores
El método riguroso proporciona cálculos de diseño, clasificación o simulación El modelo matemático de los bloques de operación de las unidades de calefacción
basados en el diseño de Aspen y y refrigeración consta de balances de materia y energía.
Diseño y simulación de procesos químicos: aplicaciones Aspen Plus y Aspen HYSYS, primera edición. Juma Haydary.
© 2019 John Wiley & Sons, Inc. Publicado en 2019 por John Wiley & Sons, Inc.
Sitio web complementario: www.wiley.com/go/Haydary/ChemDesignSimulation Aspen
Machine Translated by Google
así como de cálculos termodinámicos de fase. El balance de Aspen Plus utiliza bloques de operación de una sola unidad para
materia es muy simple y viene dado por la ecuación (3.1) si procesos de calefacción y refrigeración. La diferencia viene dada
se consideran los flujos másicos (m): únicamente por el signo del deber de calor. Aspen HYSYS utiliza
j j diferentes bloques de operación de unidades para el calentador y el enfriador.
∑ metro
=∑ metro
= min = boca = m (3.1)
yo, j
yo, en
yo, j
yo fuera
Otra diferencia es que Aspen HYSYS requiere un flujo
de energía definido, mientras que la definición de flujo de
donde i representa flujos, j representa componentes, in
energía en Aspen Plus es opcional.
representa flujos de entrada y out flujos de salida.
−1
Para el balance de entalpía podemos Ejemplo 3.1 5000 kg h de una mezcla binaria equimolar
escribir q = m (hout − hin) + qloss (3.2) de etanol y ácido acético se calientan a 101,3 kPa. La
temperatura inicial de la mezcla es de 20 ◦C y la presión
Creer
de la corriente de alimentación es de 110 kPa. Calcule el
no = cpin dT + Δvhin (3.3) requerimiento de calor usando Aspen Plus si
∫
Golpear 1. Se calienta la mezcla a 40 ◦C.
Todo
2. la mezcla se calienta y se produce vapor saturado.
3 Intercambiadores de calor 65
Abra la hoja de especificaciones del flujo de entrada haciendo clic en panel de navegación a la izquierda o haciendo doble clic en el bloque
SIGUIENTE o seleccionando el nombre del flujo en el menú principal a la CALENTADOR.
se necesita Por supuesto, este cálculo es muy simple y también se Ejecute la simulación y muestre los resultados de la misma manera
puede realizar a mano. La temperatura media del proceso es de 30 que en el caso a. El valor calculado de la carga térmica en el caso b
◦C. es de 1.098,7 kW, que lógicamente es mucho más alto que los 58 kW
La capacidad calorífica específica del etanol a esta temperatura logrados en el caso a. La temperatura de la corriente de salida es de
La temperatura es 2,4 kJ kg−1 −1 , y la capacidad calorífica 104,7 ◦C, que es la temperatura de ebullición de equilibrio de la
K La capacidad calorífica promedio de esta −1 . específica Basada en esto mezcla binaria a 101,3 kPa.
mezcla binaria es 2,1 kJ kg −1 K y el valor de la carga térmica
es Q = 5000 × 2,1 × 20 = 210−1
000 kJ h , que es alrededor de 58 kW.
−1 , 3.2 Modelos de intercambiadores de calor simples
3 Intercambiadores de calor 67
y se deben usar dos chorros de salida. El intercambiador de calor de Seleccione el modelo de intercambiador de calor de la paleta de
carcasa y tubos es el tipo de intercambiador de calor más utilizado en la modelos; si no se muestra la paleta del modelo, muéstrela seleccionando
industria de procesos. HeatX en Aspen Plus y Heat Exchanger en Aspen Paleta del modelo en la barra de menú Ver.
HYSYS son los bloques de operación de unidad correctos para modelar El icono utilizado se puede cambiar haciendo clic con el botón derecho
intercambiadores de calor de dos lados de carcasa y tubos. Si solo se y seleccionando Cambiar icono; los iconos disponibles se mostrarán
requiere el balance de energía, se puede usar el modelo ShortcutDesign como se muestra en la Figura 3.6.
en Aspen Plus o el modelo End Point en Aspen HYSYS. Sin embargo, en Abra la hoja de conexiones y especificaciones del intercambiador de
el tipo de cálculo HeatX, se debe seleccionar Diseño y modelo como calor haciendo doble clic en el icono HE. Aparece la página que se
acceso directo. En el intercambiador de calor en Aspen HYSYS, se debe muestra en la Figura 3.7.
seleccionar el modelo de punto final simple. Otra opción es usar modelos Los flujos de material de entrada y salida se pueden definir y conectar
de calentadores y enfriadores y su interconexión con flujos de energía. a HE directamente en esta página. Escriba los nombres elegidos para
las corrientes de entrada y salida en las posiciones apropiadas.
Seleccione el lado de los tubos para la mezcla de hidrocarburos, que
es la corriente caliente, y el lado de la coraza para el agua de
Ejemplo 3.2 Se enfría una mezcla −1 de un benceno y tolueno enfriamiento, aunque no es importante en este ejemplo porque el
de 6000 kg h que contiene 44 mol% de benceno mediante agua de objetivo es solo adquirir el balance de entalpía del proceso.
enfriamiento desde su temperatura de ebullición a 101 kPa hasta 30 ◦C.
Usando Aspen HYSYS, calcule el flujo de masa de agua requerido si la En la página Parámetros de diseño, se debe seleccionar el modelo de
temperatura inicial del agua es de 15 ◦C y la temperatura aumenta en HE punto final simple y la fuga/pérdida de calor se debe establecer en
3 Intercambiadores de calor 69
dos pueden ser definidos por el usuario, porque el tercero Se debe tener en cuenta la ecuación. En la forma diferencial, el
depende de los otros dos. Si VF=0, la corriente es un líquido calor transferido, dq, es proporcional al área superficial isotérmica
con punto de ebullición; si VF=1, la corriente es vapor saturado. perpendicular a la dirección de transferencia de calor, dA, la
Por la misma razón, el flujo de material de la corriente se puede diferencia de temperatura en la dirección de transferencia de calor
definir solo por un valor, por ejemplo, si se define el flujo molar, entre la corriente caliente, (Th ), y la corriente fría, ( Tc ), (Th −Tc )
no se puede definir el flujo másico o el flujo volumétrico estándar. y el coeficiente global de transferencia de calor local, U.
intercambiadores de calor
ΔT1 − ΔT2
ΔTlm = (3.7)
Si el objetivo del diseño o simulación no es solo el balance de ΔT1 ln ΔT2
energía, la solución de la tasa de transferencia de calor
Machine Translated by Google
Mientras que ΔT1 = Th,in − Tc,in y ΔT2 = Th,out − Tc,out para flujo productos es de 200 ◦C y la presión de esta corriente caliente es de 115
paralelo, ΔT1 = Th,in − Tc,out y ΔT2 = Th,out − Tc,in para flujo a kPa. La temperatura inicial de la mezcla de etanol y ácido acético es de
contracorriente. Para intercambiadores de calor con múltiples pasos por 20 ◦C y la presión de esta corriente es de 110 kPa.
el lado de los tubos y/o por el lado de la carcasa, la ecuación (3.7) se La corriente caliente (productos de reacción) contiene 40% mol de acetato
corrigió utilizando el factor de corrección F: de etilo, 40% mol de agua, 10% mol de etanol y 10% mol de ácido acético.
Usando Aspen Plus, calcule el área de intercambio de calor requerida para
q = UAFΔTlm (3.8)
calentar la mezcla binaria a 80 ◦C. Considere un coeficiente de transferencia
El factor F se puede calcular como de calor general promedio U = 200 W m−2 K temperatura de salida.
−1 . Calcular también la corriente caliente
(1 − S)
√ R2 + 1 ln
(1 − RS)
F= (3.9)
Solución:
2 S ( R + 1 √ R2 + 1 )
(R − 1) en
Comience con las mismas operaciones que en el Ejemplo 3.1; en lugar
2 S ( R + 1 + √ R2 + 1 ) del modelo Heater, use el modelo HeatX, que requiere al menos dos
T S.S − T flujos de material de entrada y dos de salida. Conecte corrientes
en afuera
R= (3.10) calientes y frías en las posiciones apropiadas.
T C −T C
afuera
en
C −T C
Defina ambas corrientes de entrada como se muestra en el Ejemplo
fuera _ en
S= (3.11)
T h −T C 3.1 (Figura 3.2).
en en
Abra la hoja de especificaciones del bloque HE haciendo clic en
Combinando las ecuaciones (3.7)–(3.11) con el balance de energía del SIGUIENTE o seleccione el nombre del bloque en el menú principal de
intercambiador de calor, se puede estimar el valor de UA. Para un valor la izquierda o haciendo doble clic en el bloque HE; aparece la página
conocido de U, se puede calcular el área de intercambio de transferencia que se muestra en la Figura 3.10.
de calor. El bloque de operación de la unidad HeatX en Aspen Plus Seleccione el tipo de modelo Clasificación (paso 1 en la Figura 3.10),
permite calcular UA y A para valores específicos de U, si el tipo de modelo el tipo de cálculo Atajo (paso 2) y la disposición de flujo Contracorriente
Clasificación y el tipo de cálculo Acceso directo están configurados. Para (paso 3).
calcular UA en Aspen HYSYS, se pueden utilizar la operación de la unidad Ajuste la temperatura de salida del chorro frío a 80 ◦C (pasos 4 y 5).
del intercambiador de calor y la calificación simple o los modelos
ponderados simples. Ingrese un valor inicial del área de la superficie del intercambiador de
calor (paso 6).
−1
Ejemplo 3.3 Se calientan 5000 kg h de una mezcla binaria equimolar de Mueva la página de métodos U y seleccione Valor U constante, luego
etanol y ácido acético mediante productos de reacción de esterificación a ingrese el valor de U (Figura 3.11); el método U predeterminado se
101,3 kPa. Temperatura de reacción establece en los valores específicos de la fase
Machine Translated by Google
3 Intercambiadores de calor 71
Figura 3.11 Especificación del método del coeficiente global de transferencia de calor en Aspen Plus
Machine Translated by Google
método, pero en este ejemplo se conoce el valor constante superar las resistencias de la capa límite de fluido caliente, la
promedio de U. Por lo tanto, se debe seleccionar el método resistencia de la pared de separación y la resistencia de la capa
del valor U constante. límite de fluido frío.
Compruebe los resultados; los resultados se proporcionan en las El coeficiente global de transferencia de calor se puede obtener
páginas Resultados térmicos y Resultados de flujo. Para mostrar utilizando la siguiente fórmula:
los resultados requeridos en este ejemplo, abra los resultados
1 1
térmicos: página de resumen y detalles del intercambiador. Como en = =
r
se muestra en la Figura 3.12, la temperatura final de la corriente
. es de 7,9 m2
caliente es de 130 ◦C y el área requerida del intercambiador de calor 1 hacer ln ( hacer di ) 1
+f s + t +
h h yo 2k + Ar ( f ht )
(3.12)
3.4 Diseño Detallado y Simulación de
Intercambiadores de Calor donde r es la resistencia total contra la transferencia de calor
(m2 K W−1 ), hs el coeficiente de película del lado de la coraza
−1
El diseño riguroso, la simulación y la clasificación de los intercambiadores (W m−2 K ), el factor de ensuciamiento
i la eficiencia delsf lado de la carcasa (m2
de la aleta,
de calor incluyen el cálculo de los coeficientes de transferencia de calor de K W−1 ), do el diámetro exterior del tubo (m), din el diámetro interior
ambos lados de un intercambiador de calor, los llamados coeficientes de del tubo (m), k la conductividad térmica de pared del tubo (W m−1 K
−1
película. El flujo de calor, en su camino de fluido caliente a frío, tiene que ), Ar el tubo exterior/interior
Machine Translated by Google
3 Intercambiadores de calor 73
relación de t el factor de ensuciamiento del lado del tubo (m2 K W−1 ), y Tabla 3.2 Correlaciones del coeficiente de transferencia de calor del lado de la carcasa (4)
VARA
No = (3.13) Jensen (17)
Hirviendo
Condensación (9–12)
vertical
adimensionales: l kf lw Re = (3.14)
Laminado Nusselt
ondulado laminar Después
c.p.
Pr = (3.15) Turbulento Labuntsov
kf Corte Rohsenów
dominado
donde Nu, Re y Pr representan los números de Nusselt, Reynolds
Condensación + (13)
y Prandtl, respectivamente, representa la dimensión característica horizontal Kern
del dominio del intercambiador de calor (en el caso de un tubo,
es el diámetro del tubo), (kg m− 3 ) es la densidad del fluido,
−1 ) el simulación rigurosa de intercambiadores de calor si se establecen
(Pa s) la viscosidad dinámica, cp (kJ kg−1 K
capacidad calorífica específica, y kf es la conductividad térmica los parámetros Rigorous Shell and Tube o Dynamic Rating y
del fluido que fluye en el tubo (W m−1 K −1 ). Detail rating.
Diferente Se pueden usar correlaciones empíricas que
Ejemplo 3.4 El enfriamiento de una mezcla de tolueno y benceno
dependen del carácter del flujo y la geometría del intercambiador
de calor para obtener el número de Nusselt. Las correlaciones descrita en el ejemplo 3.2 se realiza en un intercambiador de calor de
carcasa y tubos existente con la geometría que se indica en la tabla
utilizadas en Aspen Plus y Aspen HYSYS se resumen en las
Tablas 3.1 y 3.2. 3.3. El intercambiador de calor se muestra esquemáticamente en la
Figura 3.13. La temperatura inicial del agua de enfriamiento es de 15
El bloque de operación de la unidad HeatX en Aspen Plus −1
◦C y su caudal másico es de 15 000 kg h . Calcule la temperatura de
proporciona una clasificación de diseño rigurosa y una simulación
las corrientes de salida, el lado de la carcasa, el lado de los tubos y el calor general
del intercambiador de calor si se establece el modo Clasificación
detallada o Simulación detallada (18). Además, Aspen EDR
Tabla 3.3 Geometría del intercambiador de calor utilizado en el Ejemplo 3.4
(Exchanger Design and Rating) está integrado en HeatX. Aspen
EDR también está integrado en el modelo de intercambiador de
AEL (AChannel y cubierta removible, E
calor en Aspen HYSYS. El modelo Aspen HYSYS Exchanger permite también
One pass shell,
Tipo de intercambiador de calor LCanal extraíble con tapa plana
Tabla 3.1 Correlaciones del coeficiente de transferencia de calor del lado del tubo (4)
Posición Horizontal
Mecanismo Régimen de flujo Correlación Referencia
Número de deflectores 5
Vistas en la flecha A
A
4235 Total
T1
220 158 170 3155
429
T1
429
T2
T2 S2 S1
447 410
158 700 2100
3180
Longitud de tracción
S1 S2
75 75
112 112
450 450
150 150
2 pernos 2 pernos
fijos deslizantes
Figura 3.13 Detalles del intercambiador de calor existente utilizado en el Ejemplo 3.4
coeficientes de transferencia y caída de presión en los lados de la carcasa y los temperatura de 15 ◦C, presión de 101 kPa y caudal másico de 15 000 kg h
tubos. Vaya a Parámetros de −1 .
Continúe en la solución del Ejemplo 3.2 cambiando la especificación de las 20,3 ◦C. Para conocer los coeficientes de transferencia de calor y la caída de
corrientes de entrada. Elimine todas las especificaciones de las corrientes presión, consulte la página de parámetros de clasificación; el coeficiente de
transferencia de calor de la carcasa es de 265 W m−2 −1 ; tubo de calor
de salida y especifique solo las corrientes de entrada, alimentación de
K el coeficiente de transferencia es de 359 −1 ; presión del lado del tubo
hidrocarburos como VF=0, presión de 101 kPa y agua de refrigeración por
caudal másico de 5.000 kg h −1 , W m−2 K la caída segura es de alrededor de 0,5 kPa y la presión lateral de la carcasa
Machine Translated by Google
la caída es prácticamente nula (alrededor de 0,04 kPa). Para conocer modelo de calificación dinámica está limitando la solución. Haga clic
la media aritmética logarítmica y los perfiles de los intercambiadores de en transferir geometría desde HYSYS y verifique si los datos son
calor, consulte la página Rendimiento; y para el coeficiente de coherentes con los parámetros del intercambiador que se dan en la
transferencia de calor global equivalente, U, consulte la página Dinámica; Tabla 3.3.
su valor es 137 W m−2 −1 . K Después de esta operación, HYSYS calcula automáticamente el
intercambiador de calor utilizando el riguroso modelo de carcasa y
tubos.
3.4.2 Diseño riguroso de intercambiadores de calor de carcasa
Para mostrar los resultados, vaya a Resumen de resultados en
y tubos mediante EDR
Rigorous Shell & Tube y para conocer las condiciones de los flujos
Para obtener una clasificación más detallada del intercambiador de calor de salida, consulte la Hoja de trabajo.
En este ejemplo, usando solo datos transferidos de cálculos previos que el calculado por el modelo de clasificación dinámica.
por el
Machine Translated by Google
3 Intercambiadores de calor 77
Figura 3.17 Selección del modelo riguroso de carcasa y tubos en Aspen HYSYS
Para comprobar el perfil de temperatura calculado por EDR (Figura e instalado. La selección del tipo de intercambiador de calor a utilizar es
3.21), consulte Perfiles en Rigorous Shell & Tube; El perfil de una de las acciones más importantes realizadas por el ingeniero de
temperatura es un perfil típico para intercambiadores de una coraza diseño de procesos. El proceso de selección normalmente incluye una
y de dos pasos de tubos. serie de factores, como los requisitos térmicos e hidráulicos, la
compatibilidad de los materiales, el mantenimiento operativo, las
consideraciones y normativas ambientales, de salud y seguridad, la
3.5 Selección y Costo de disponibilidad y el costo (19).
Intercambiadores de Calor Los principales factores que influyen en el costo de los intercambiadores
de calor son el área de transferencia de calor, el material de construcción,
La evaluación económica del diseño del proceso y la estimación del la longitud del tubo, el diámetro y el grosor del tubo, los fluidos a presión,
costo del proyecto es el tema del Capítulo 12. Aquí, nos enfocamos solo los requisitos del deflector, los requisitos especiales de diseño, etc.
en la selección del tipo de intercambiador de calor para un proceso y la Un resumen de los criterios clave para la selección preliminar del
estimación del costo del equipo comprado. intercambiador de calor y una muy buena presentación gráfica
3 Intercambiadores de calor 79
Figura 3.21 Perfil de temperatura del intercambiador de calor calculado por EDR
del costo de diferentes tipos de intercambiadores de calor fue presentado estimación de costos de equipos. El mapeo predeterminado para todos
por Peters et al. en 20). los equipos de transferencia de calor es el intercambiador de calor TEMA,
Es difícil para los ingenieros de procesos recopilar datos recientes que se puede personalizar para varios tipos diferentes de intercambiadores
sobre los costos de todos los tipos de intercambiadores de calor disponibles de calor.
en el mercado. En cambio, el método más común de estimación preliminar APEA está integrado tanto en Aspen Plus como en Aspen HYSYS. Sin
de costos es el uso de un software comercial de estimación de costos. embargo, está desactivado por defecto. Para iniciar el análisis económico,
Uno de los software más utilizados para la estimación preliminar de costos se puede activar seleccionando Economía en la barra de herramientas y
de procesos es Aspen Process Economic Analyzer (APEA). Esta marcando Económico activo.
herramienta permite mapear los equipos usados en el proceso, su
Ejemplo 3.5 5000 kg h de etanol −1 de una mezcla equimolar de
dimensionamiento en base a los parámetros del proceso y finalmente el
y ácido acético después de calentar a 80 ◦C por los productos de la
costo del proceso incluyendo la estimación de costos de equipos y equipos
reacción (vea el Ejemplo 3.3) se calienta a 250 ◦C a un vapor precalentado
instalados. Si el usuario no especifica exactamente el equipo, APEA
a 101.3 kPa. Seleccione el equipo adecuado para ambos bloques de
asigna el modelo de operación de la unidad utilizado en el diagrama de
operación de la unidad intercambiadora (HE y HEATER) y calcule su
flujo del proceso a un equipo predeterminado.
costo. Compare el costo del calentador para los siguientes dos escenarios:
1. como calentador, se usa un hervidor tipo marmita y 2.
El usuario debe verificar el mapeo predeterminado y seleccionar un tipo
de equipo apropiado para un modelo de operación de unidad dado; de lo como calentador, se usa un horno tipo caja.
contrario, se pueden introducir grandes errores en el
Machine Translated by Google
Solución:
Figura 3.23. Para el bloque HE, el intercambiador de calor TEMA
Continúe en el Ejemplo 3.3 agregando un bloque calentador seleccionado por defecto es el tipo correcto. Al hacer clic en Aceptar
(CALENTADOR) al diagrama de flujo principal. Como corriente de (paso 7), se completa el proceso de mapeo.
entrada, use la corriente de salida fría del bloque HE; defina también El siguiente paso es dimensionar, esta acción se proporciona al
la corriente de salida. hacer clic en Tamaño como se muestra en la Figura 3.24. Seleccione
Ingrese a la página del bloque CALENTADOR y especifique la el nombre del equipo (paso 2 en la Figura 3.24) y verifique los
presión de salida (101.3 kPa) y la temperatura de salida (250 ◦C). parámetros de tamaño calculados por APEA. Estos parámetros se
Ejecute la simulación y verifique los resultados. pueden cambiar simplemente reescribiendo sus valores si es
Si el analizador económico está inactivo, Aspen Plus solicita su necesario
activación después de ejecutar la simulación. Para activar el Seleccione el icono Evaluar de la barra de herramientas (paso 1 en
analizador económico, haga clic en Activar en este cuadro de la Figura 3.25) para proceder con la evaluación económica del
diálogo o marque Activar económico en la barra de herramientas proceso por parte de APEA. Tenga en cuenta que en este ejemplo,
de Economía como se muestra en la Figura 3.22. Esta opción nos enfocamos solo en el costo del equipo; la evaluación
permite un uso muy rápido y eficaz de la APEA en Aspen Plus. Otra económica de otros componentes se discutirá más adelante.
opción es enviar la simulación a APEA usando Send to Economics. Para ver los resultados, haga clic en Ver equipo (paso 2 en la
Si elige esta opción, el software APEA se ejecuta y puede realizar Figura 3.25) y luego en Equipo (paso 3 en la Figura 3.25).
las operaciones necesarias dentro de este software.
Como resultado de la Figura 3.25, los costos básicos calculados
Cuando el analizador económico esté activado, haga clic en Mapa por APEA son 10 200 USD para HE y 28 900 USD para el rehervidor
en la barra de herramientas de Economía; aparece la página Kettle, respectivamente. El costo de instalación de HE es de 62.500
Opciones de mapa; haga clic en Aceptar y continúe. USD y el de Kettle Reboiler es de 114.900 USD.
En la página de vista previa del mapa que aparece, seleccione
CALENTADOR (paso 1 en la Figura 3.23). Como tipo de equipo Para cambiar el tipo de equipo para el bloque CALENTADOR de
predeterminado, se selecciona el intercambiador de calor TEMA hervidor de agua a horno de caja, repita todos los pasos desde el
(DHE TEMA EXCH); para cambiar este tipo de intercambiador a mapeo, pasando por el dimensionamiento y la evaluación. En el
Hervidor de agua con cabeza flotante, siga los pasos en paso de ping del mapa, cambie Hervidor de agua a Horno de caja.
Machine Translated by Google
3 Intercambiadores de calor 81
Como resulta de la Figura 3.25, si el bloque HEATER se es 3,4 veces mayor que el costo del hervidor Kettle. El
mapea a un horno de caja, el precio del equipo es de 279 costo del equipo es unas diez veces mayor que en el caso
700 USD y el costo de instalación es de 389 500 USD, lo que del calderín Kettle.
Referencias
1 Rohsenow WM, Hartnett JP, Cho YI. Manual de transferencia de calor, 6 Steiner D, Taborek J. Flujo de transferencia de calor por ebullición en
3ª ed. Nueva York: McGrawHill; 1998. tubos verticales correlacionados por un modelo asintótico. Ingeniería de
transferencia de calor. 1992;13(2):43–69.
2 Lienhard JH, IV, Linhard JH, V. A Transferencia de calor 7 Sha MM. Una nueva correlación para la transferencia de calor durante
Libro de texto, 3ª ed. Cambridge, MA: Prensa de flogisto; 2008. flujo en ebullición a través de tuberías. ASHRAE Trans. 1976;82(2): 66–86.
3 Verde DW, Perry RH. Manual de ingenieros químicos de Perry, 8ª ed. 8 Sha MM. Correlación de gráficos para transferencia de calor de ebullición
Nueva York: McGrawHill; 2008. saturada: ecuaciones y estudio adicional. ASHRAE Trans. 1981;87(1):185–
4 Ayuda de Aspen HYSYS V9. Bedford, MA: Aspen Tech. Cª.; 196.
2016. 9 Nusselt W. Condensación superficial de vapor de agua. Versión Z.
5 Gnielinski V. Convección Forzada en Ductos. En celo Dtsch, Ing. 1916;60(27):541546.
Manual de diseño de intercambiadores. Nueva York: Hemisphere 10 Kutateladze SS. Fundamentos de la Transferencia de Calor. Nuevo
3 Intercambiadores de calor 83
11 Labuntsov DA. Transferencia de calor en condensación pelicular de 16 Gentry CC. Tecnología de intercambiador de calor RODBaffle.
vapor puro sobre superficies verticales y tubos horizontales. química Ing. prog. 1990;86(7):48–57.
Teploenergetika 1957;4(7):72–80. 17 Jensen MK, Hsu JT. Un estudio paramétrico de la transferencia de calor
12 Rohsenow WM, Webber JH, Ling AT. Efecto de la velocidad del vapor en ebullición en un haz de tubos. En 1987 ASMEJSME
en la condensación de película laminar y turbulenta. Thermal Engineering Joint Conference, Honolulu, HI, 1987, pp. 133–
Trans. ASME 1956;78:1637–1643. 140.
13 Jaster H, Kosky PG. Transferencia de calor por condensación en un 18 Ayuda de Aspen Plus V9. Bedford, MA: Aspen Tech. Cª.;
régimen de flujo mixto. En t. J. Transferencia de masa de calor 2016.
1976; 19:95–99. 19 Kakac¸ S, Liu H, Pramuanjaroenkij A. Intercambiadores de calor:
14 Taborek J. Intercambiadores de calor de carcasa y tubos: Flujo monofásico. selección, clasificación y diseño térmico, 3.ª ed. Boca Ratón, FL:
En Manual de diseño de intercambiadores de calor. Nueva York: CRC Press; 2012.
Hemisphere Publishing Corporation; 1983. 20 Peters M, Timmerhaus K, West R. Diseño y economía de plantas
15 campana KJ. Método de Delaware para el diseño del lado de la carcasa. para ingenieros químicos, 5.ª ed. Nueva York: McGrawHill; 2004.
En Kakac S, Bergles AE, Mayinger F, editores, Heat Exchangers:
ThermalHydraulic Fundamentals. Berlín: Springer; 1985.
Machine Translated by Google
85
4.1 Bombas, turbinas hidráulicas y válvulas El cabezal de succión positivo neto (NPSHA) disponible para un
bomba se define como
Básicamente, una bomba aumenta la presión del líquido entre su entrada
2
y salida en ΔP:
NPSHA = P1 − P ◦ en + + Hs (4.6)
2g
P2 = P1 + ΔP (4.1)
Las relaciones entre el caudal volumétrico de líquido y la cabeza de la 1. flash adiabático para presión de salida específica (presión
bomba, la eficiencia de la bomba, la potencia del fluido o la potencia del cambiador),
freno están representadas por las curvas características de la bomba, las 2. calcular el coeficiente de flujo de la válvula para la presión de salida
denominadas curvas de rendimiento de las bombas. En la Figura 4.1 se especificada (diseño), y 3.
muestra una muestra de la curva de rendimiento de una bomba centrífuga. calcular la presión de salida para la válvula especificada (valor nominal).
Diseño y simulación de procesos químicos: aplicaciones Aspen Plus y Aspen HYSYS, primera edición. Juma Haydary.
© 2019 John Wiley & Sons, Inc. Publicado en 2019 por John Wiley & Sons, Inc.
Sitio web complementario: www.wiley.com/go/Haydary/ChemDesignSimulation Aspen
Machine Translated by Google
5.00
Cabeza
4.00
P.f.
3.00
N.f.
2.00
1.00
0.00
0123456 7 8
q
Para cambiar la presión de una corriente de gas desde la entrada Tabla 4.2 Composición del gas natural utilizado en el Ejemplo 4.2
compresor, hasta la presión de salida, P2 , la presión del
P1 , desarrolla una cabeza dada por Componente CH4 C2H6 C3H8 N2 CO2
P2
Mol% 92 3 2 2.5 0.5
Hcom = VdP (4.8)
∫
P1
donde V es el volumen molar y P1 y P2 son las presiones de ing compresores. Para compresores centrífugos, se puede utilizar
entrada y salida, respectivamente. eficiencia adiabática (isentrópica) o politrópica. Como método
La ecuación (4.8) solo es correcta para gases perfectos. Para el politrópico, se puede seleccionar Schultz (7), Huntington (8) o
cambio de entalpía real por mol de gas, se debe usar el factor de Reference (integración por partes).
eficiencia: h ,
Ejemplo 4.2 Se utiliza un compresor centrífugo para cubrir la caída
Hcom
Δh = (4.9) de presión durante el transporte de gas natural entre dos estaciones
h
compresoras. La temperatura de entrada del gas natural es de 25
En un proceso de compresión politrópica, la relación de ◦C, la presión de entrada es de 5,5 MPa y la relación de presión del
P y V viene dada por compresor es de 1,5. Se procesa un millón de Nm3 por día de gas
natural con la composición dada en la Tabla 4.2. Usando Aspen
PVn = constante (4.10) HYSYS, calcule la temperatura del gas a la salida del compresor,
la carga de fluido politrópico y la potencia requerida, si la eficiencia
donde n es el exponente politrópico. Suponiendo una constante n,
politrópica del compresor es 0.76.
para el cambio de entalpía real de compresión politrópica tenemos
recibir:
Solución:
P1V1
Δh = (4.11)
P2−1
P1
norte ]__
[( )n−1
Abra Aspen HYSYS, seleccione una lista de componentes y el
pags ( norte
norte 1 )
paquete de fluido apropiado (Peng–Robinson en este caso)
donde es la eficiencia politrópica del compresor. como se explica en los Capítulos 1 y 2.
pag
Cambie al entorno de simulación.
Para un proceso isotérmico, n=1, mientras que para un isen donde
Seleccione dos flujos de material de la paleta del modelo como
isotrópicos n= , =cp /cv , capacidad cp y cv son procesos
se muestra en la Figura 4.7; abra uno de ellos haciendo doble
calorífica bárica e isocórica, respectivamente.
clic en él e ingrese la especificación: nombre de la corriente
Para los compresores de una etapa, Aspen Plus proporciona el
(entrada), temperatura, presión y flujo molar (Figura 4.8).
bloque de operación de la unidad Compr, que incluye modelos
isoentrópicos, politrópicos y de desplazamiento positivo. Los
Vaya a la página Composición e ingrese la composición del gas
métodos de cálculo utilizados para tipos de modelos individuales
son los siguientes: como se explica en el Ejemplo 3.2; el color verde en la parte
inferior de la página indica que la especificación de la transmisión
Isentrópico: basado en Mollier (equivalente al gráfico de Mollier), está completa.
GPSA (Asociaciones de proveedores de procesadores de gas) Para la segunda transmisión, especifique solo el nombre de la
(5) y ASME (Sociedad estadounidense de ingenieros mecánicos) transmisión (Outlet).
(6); Politrópico: Seleccione un modelo de compresor y conecte el flujo de entrada
GPSA, ASME, Integración; Desplazamiento a la posición de Entrada y el flujo de salida a la posición de
positivo: GPSA, Integración. Salida en la página DiseñoConexión; defina también un flujo de
energía directamente en esta página.
Compr también se puede utilizar como un modelo de turbina de gas En la hoja de Parámetros de Diseño, seleccione el modo de
seleccionando Turbina en la página Configuración. funcionamiento del compresor como centrífugo; definir la relación
Para simular la compresión de varias etapas en Aspen Plus, se puede de presión; elimine el conjunto predeterminado de eficiencia
utilizar el bloque de operaciones de la unidad MCompr. adiabática y escriba el valor de la eficiencia politrópica; y aplicar
Aspen HYSYS utiliza dos modelos diferentes para compresores el método de Shultz como método politrópico.
y turbinas de gas. El modelo Compressor se usa para simular Los pasos necesarios para hacer se muestran en la Figura 4.9.
compresores de una etapa y Expander se usa para simular turbinas Los resultados requeridos se pueden verificar en la hoja de
de gas y expansores. El compresor permite modelar centrífugos y trabajo y la página de rendimiento; como se muestra en la figura
recíprocos 4.10, la temperatura de salida del gas es de 62.7 ◦C,
Machine Translated by Google
la potencia requerida (potencia consumida) es de 642,5 kW y la Si resolvemos el mismo ejemplo usando Compr en Aspen Plus y
cabeza de fluido politrópico es de 54,15 kJ kg−1 . seleccionamos el método politrópico GPSA, la temperatura de salida
calculada es de 58,2 ◦C y la potencia requerida es de 647 kW. Teniendo
Cambiar el método politrópico solo tiene un efecto insignificante en en cuenta las diferencias entre el paquete de fluidos PengRobinson
los resultados. HYSYS calcula también la eficiencia adiabática, que utilizado en Aspen HYSYS y el método PengRob en Aspen
ronda el 74,8%.
Además del método de propiedad (consulte el Capítulo 2) y la cantidad de fricción y un componente cinético o de aceleración de la caída de presión:
gas procesado, la coherencia entre los resultados de Aspen Plus y HYSYS
es bastante buena. dP dw
+ mw dL (4.12)
dL = mg sen Θ + ( dP dL) f
r
Colebrook (9), Chen (10), etc. La presencia de múltiples fases complica asientos: una válvula de compuerta abierta, una válvula de compuerta
significativamente los cálculos de caída de presión ya que se deben tener semiabierta, dos codos estándar de 90°, un medidor de agua de disco y
en cuenta las propiedades de cada fluido presente. una válvula de globo abierta, que se colocan como se muestra en la figura
4.11. La entrada y la salida de la tubería están en el mismo nivel horizontal.
La presión del agua al comienzo del sistema de tuberías es de 5 bar y su
Aspen HYSYS calcula la pérdida de presión total a lo largo de la tubería temperatura es de 25 ◦C. La temperatura ambiente promedio es de 20 ◦C,
como la suma del gradiente de presión causado por la fuerza gravitacional, y el valor promedio del coeficiente de transferencia de calor es
las pérdidas de presión causadas por el fluido
Tabla 4.3 Correlaciones utilizadas para el cálculo de la caída de presión de la tubería en Aspen HYSYS
(continuado)
Machine Translated by Google
20 W m−2 K −1 . Calcule la caída de presión, la presión y el perfil modelo. En esta simulación, solo se usa agua como componente,
de temperatura a lo largo de la tubería. por lo que se puede usar el paquete de propiedades de Antoine.
Solución:
Después de la conexión de la corriente, seleccione la correlación
Siga los mismos pasos iniciales que en el Ejemplo 4.2, pero en de flujo de la tubería y elija Beggs y Brill (11) para esta simulación.
lugar de Compressor use el segmento Pipe
Disco Válvula
Medidor de globo:
de agua Abierto
50m 50m
50m
Válvula Válvula
de compuerta: de compuerta:
Medio
Abierto 50m 50m
Figura 4.13 Especificación de las condiciones de transferencia de calor desde el sistema de tuberías
Machine Translated by Google
Como se muestra en la Figura 4.12, defina los segmentos de temperatura y coeficiente global de transferencia de calor; seleccione
tubería. Para cada segmento de tubería y elemento de conexión, HTC general e ingrese la temperatura ambiente y HTC general
defina un segmento separado (juntos 11) en el mismo orden en que como se muestra en la Figura 4.13.
se colocan en el sistema de tuberías que se muestra en la Figura 4.11. Aspen HYSYS calcula las condiciones de la corriente de salida, así
Aspen HYSYS ofrece diferentes posibilidades de especificar la como los diferentes perfiles a lo largo del sistema de tuberías. Para
transferencia de calor entre el sistema y el entorno ambiental. En comprobar la presión de salida y la temperatura de salida, consulte
este ejemplo, sabemos ambiente la Hoja de trabajo (Figura 4.14). Salida
la presión de la corriente es de 4,33 bar; esto significa una caída de Tabla 4.4 Rango de operación de algunos tipos de bombas
presión de 0,76 bar en el sistema de tuberías.
Para mostrar el perfil de presión o el perfil de otros parámetros a lo largo Rango de capacidad
de la tubería, use el enlace PerformanceView Profiles. Los perfiles se Tipo (m3 h −1 ) Cabeza (m)
0,05–50 5–60
cambio de presión Diafragma
Admirador
Compresorgas CompresorAire
Admirador
alternativo Centrífugo
Centrífugo
Centrífugo
Soplador Receta.
compensación
Con Con
Horizontal turbina motor
Hélice
Parámetro
Cuando el analizador económico esté activado, haga clic en Mapa Seleccione el ícono Evaluar de la barra de herramientas para
en la barra de herramientas de Economía; aparece la página proceder con la evaluación económica del proceso por parte de
Opciones de mapa (Figura 4.17); haga clic en Aceptar y continúe. APEA.
En la página Vista previa del mapa que aparece, el compresor Con el icono Ver equipo, verifique los resultados.
CentrífugoHorizontal está seleccionado como el tipo de equipo Repita el mapeo, dimensionamiento y evaluación para el
predeterminado; para cambiar este tipo de compresor a otro, compresor alternativo y el compresor alternativo con el motor
siga los pasos de la Figura 4.18 haciendo clic en Aceptar (paso integral.
7), el proceso de mapeo se completa.
APEA proporciona resultados de dimensionamiento y costos
Proporcione el tamaño del compresor haciendo clic en Tamaño; para equipos, costos de servicios públicos y costos operativos
la información más importante de mapeo y dimensionamiento se unitarios. La Tabla 4.6 muestra los costos, el peso y los costos de
muestra en la Tabla 4.5. Para información más detallada, se debe utilidad de los tres tipos de compresores. Como resulta de esta
utilizar el método alternativo de enviar la simulación a APEA y tabla, APEA proporcionó los costos más bajos para el compresor
trabajar en APEA. alternativo con el motor de gas integral.
Referencias
10 Chen NH. Una ecuación explícita para el factor de fricción en la tubería.
1 DW verde, Perry RH. Manual de ingenieros químicos de Perry, 8ª ed.
Nueva York: McGrawHill; 2008. Ing. Ind. química Fundam. 1979;18(3):296–297.
11 Mendicidad HD, Brill JP. Un estudio de flujo de dos fases en
2 ordenadores Hanlon. Manual del compresor. Nueva York: McGraw tuberías inclinadas. J. Mascota. Tecnología 1973;25(5):607–617.
Hill; 2001. 12 Payne GA, Palmer CM. Brill JP, Beggs HD. Evaluación de tubería inclinada,
retención de líquido bifásico y correlación de pérdida de presión
3 Ayuda de Aspen Plus®V9. Burlington, MA; Tecnología Aspen,
utilizando datos experimentales.
Inc.; 2016. www.aspentech.com.
J. Mascota. Technol, 1979;31(9):1198–1208.
4 Ayuda de Aspen HYSYS®V9. Burlington, MA: Aspen Technology,
Inc.; 2016. www.aspentech.com. 13 Duns H, Ros CNJ. Flujo vertical de mezclas de gas y líquido en
5 GPSA (Asociaciones de Proveedores de Procesadores de Gas) pozos. En Actas del Sexto Congreso Mundial del Petróleo, Frankfurt,
Libro de datos de ingeniería. Tulsa, OK: Asociaciones de proveedores Alemania, junio de 1963.
de procesadores de gas; 1979. Capítulo 4. pp. 56 a 510. Sección II, Documento 22PD6, págs. 19–26.
14 Orkiszenolwski J. Predicción de caídas de presión bifásicas en tubería vertical.
J. Mascota. Tecnología 1967;19(6):829–838.
6 ASME (Sociedad Americana de Ingenieros Mecánicos)
Código de prueba de potencia 10. Ciudad de Nueva York: Sociedad 15 Aziz K, Govier GW, Fogarasi M. Caída de presión en pozos productores de
Estadounidense de Ingenieros Mecánicos; 1965, págs. 31 y 32. petróleo y gas. J. Can. Mascota. Tecnología 1972;11(3):38–
48.
7 SchultzJM. El análisis politrópico de compresores centrífugos. J. Ing.
Fuerza. 1962;84(1):69–82. 16 Hagedorn AR, Marrón KE. Estudio experimental de los gradientes de
8 Huntington AR. Evaluación del cálculo politrópico presión que se producen durante un flujo bifásico continuo en
conductos verticales de pequeño diámetro.
Métodos para el rendimiento de turbomaquinaria. J. Ing. Gas Turbine Power
1985;107(4):872–876. J. Mascota. Tecnología 1965;17(4):475–484.
9 Colebrook CF. Flujo turbulento en Tuberías, con particular referencia a la 17 Lockhart RW, Martinelli RC. Correlación propuesta de datos para flujo
isotérmico bifásico de dos componentes en tuberías. química Ing. prog.
región de transición entre las leyes de tuberías lisas y rugosas. J. HIELO.
1949;45(1):39–48.
1939;11(4):133–156.
Machine Translated by Google
18 Mandhane JM, Gregory GA, Aziz K. Evaluación crítica de los métodos 20 Davidson J, Bertele O. Selección de bomba de proceso: A
de predicción de caída de presión por fricción para el flujo Enfoque de sistemas, 2ª ed. Chichester, Reino Unido: John Wiley
de gaslíquido en tuberías horizontales. J. Mascota. and Sons, Ltd.; 2005.
Tecnología 1977;29(10):1348–1358. 21 Towler G, Sinnott R. Diseño de ingeniería química,
19 Peters MS, Timmerhaus KD, West RE. Diseño y economía de Principio, práctica y economía del diseño de plantas y procesos,
plantas para ingenieros químicos, 5.ª ed. Nueva York: McGraw 2ª ed. Ámsterdam, Países Bajos: Elsevier; 2013.
Hill; 2004.
Machine Translated by Google
101
reactores
5.1 Balance de materia y entalpía de un reactor entalpía de la corriente de salida k, ΔrHn reacción entalpía
químico de la reacción n, Q flujo de calor agregado o eliminado del
sistema (si se agrega calor al sistema, Q tiene un signo más
y si se elimina del sistema Q tiene un signo menos ), y W es
El balance de materia del componente i que participa en n reacciones en
la tasa de trabajo agregado al sistema (signo más) o
un reactor con corrientes de entrada j y corrientes de salida k en
realizado por el sistema (signo menos) (en muchas
condiciones de estado estacionario se puede escribir de la siguiente manera:
aplicaciones, este término puede despreciarse).
Se pueden aplicar reglas aditivas para calcular la entalpía y los
=0
parámetros de dependencia de la temperatura, como las
(∑ nuevo)
j − (∑ knik ) + (∑ n en n ) R
en afuera
capacidades caloríficas.
(5.1) El modelado de reactores químicos en Aspen Plus y Aspen
HYSYS requiere una buena formación en ingeniería de
donde nij representa el flujo molar del componente i en la corriente
reacciones químicas. Se encuentran disponibles varios libros
de entrada j, nik es el flujo molar del componente i en la corriente de
de texto de ingeniería de reacciones químicas que estudian los
salida k, in es el coeficiente estequiométrico del componente i en la n
detalles de diferentes tipos de reactores químicos (1, 2).
ésima reacción, y es la extensión de la reacción de la nésima norte
∑ Posada = adentro, adentro − adentro, afuera (5.2) Los modelos de estequiometría (RStoic en Aspen Plus y
norte
Conversion Reactor en Aspen HYSYS) se utilizan básicamente
La conversión de la reacción i en reacciones químicas individuales cuando se conoce la estequiometría de reacción y la
se puede expresar mediante
conversión o extensión molar de cada reacción, y la cinética
de reacción es desconocida o no es importante. RStoic y
= ni, reaccionó = ni,0 − ni,t
(5.3) Conversion Reactor pueden modelar reacciones que ocurren
i
si, 0 si, 0 simultánea o secuencialmente. Además, estos modelos
donde ni,0 es el flujo molar inicial (moles en el caso de reactores
pueden calcular el calor de reacción. Estos cálculos en
discontinuos) del componente i, ni,reaccionado representa los moles
modelos de reactores estequiométricos se basan en las
(flujo molar) del componente i reaccionado y ni,t el flujo molar (moles)
ecuaciones de balance de materia y energía (5.1)–(5.4).
del componente i en la mezcla de reacción en el tiempo t. La conexión de RStoic requiere al menos un flujo de
material de entrada y un flujo de material de salida. Más flujos
Considerando las condiciones de estado estacionario, una energía general de material de entrada, cualquier cantidad de flujos de energía
el saldo de un rector se puede escribir como de entrada y salida y un flujo de salida de agua libre son opcionales.
La conexión del reactor de conversión en Aspen HYSYS
requiere al menos un flujo de material de entrada y dos
(∑ njhj ) − (∑ k nkhk ) flujos de material de salida. Si se modela un reactor no
j en afuera
adiabático, también es obligatoria la conexión de una
+∑ ΔrHn norte +Q+W=0 (5.4) corriente de energía.
norte
RStoic permite la definición de la estequiometría de la reacción
donde nj es el flujo molar de la corriente de entrada j, nk el flujo molar directamente dentro del modelo de operación unitaria, mientras
de la corriente de salida k, hj entalpía de la corriente de entrada j, hk que en el caso del Reactor de conversión de HYSYS las reacciones
Diseño y simulación de procesos químicos: aplicaciones Aspen Plus y Aspen HYSYS, primera edición. Juma Haydary. © 2019
John Wiley & Sons, Inc. Publicado en 2019 por John Wiley & Sons, Inc.
Sitio web complementario: www.wiley.com/go/Haydary/ChemDesignSimulation Aspen
Machine Translated by Google
debe definirse en el entorno Propiedades, agruparse en un conjunto de Cambie al entorno de simulación y prepare el diagrama de flujo del
reacciones y agregarse al paquete de fluidos (consulte el Ejemplo 2.12); proceso con el mismo método que se muestra en el Ejemplo 2.14; solo
dentro del modelo de operación unitaria, el conjunto de reacciones necesitamos un bloque Rstoic, que se puede seleccionar de la paleta
definido simplemente se agrega al modelo del reactor. del modelo Reactors; el diagrama de flujo se completa dibujando dos
flujos de material de entrada y uno de salida (ver Figura 5.1).
Si no se dispone de información sobre las reacciones y su
estequiometría, pero sí la cantidad de componentes individuales Especifique las corrientes de entrada ingresando la temperatura, la
producidos por unidad de masa o unidad de mol de alimentación del presión, la composición y el flujo molar, para el flujo molar de agua use
reactor, use RYield en Aspen Plus o Yield Shift Reactor en Aspen HYSYS. un valor 12 veces mayor que para el flujo molar de óxido de etileno.
RYield proporciona dos opciones principales: especificación de Especifique el bloque Rstoic ingresando la presión y la temperatura
rendimientos de componentes y especificación de ping de asignación (Figura 5.2).
de componentes. Si se selecciona la opción de rendimiento de Defina la reacción química usando la página Reacción; siguiendo los
componentes, los rendimientos de los productos deben especificarse o pasos que se muestran en la Figura 5.3; definir la estoiquiometría de
calcularse en una subrutina Fortran proporcionada por el usuario. RYield la reacción y conversión; repítelo para la segunda reacción.
normaliza los rendimientos para mantener un balance de masa. Los
rendimientos se pueden especificar como moles de un componente por La simulación está preparada para ejecutarse; sin embargo, para
unidad de masa de alimentación total o como masa de un componente obtener la relación requerida de flujo molar de óxido de etileno a agua
por unidad de masa de alimentación total. Los rendimientos de los para cualquier otro flujo molar de óxido de etileno, se puede definir un
componentes no convencionales deben especificarse en masa. bloque de cálculo.
conversión de C2H4O: 95% kmol h Tenga en cuenta que para la corriente mole flow StreamVar
debe seleccionarse como tipo, y para el componente mole flow Mole
HOCH2 CH2OH + C2H4O → (HOCH2 CH2 )2O (R5.2)
Flow debe seleccionarse como tipo.
Defina el flujo molar de la corriente de agua por el mismo método.
conversión de C2H4O: 5%
Vaya a la página de Cálculo y, usando comandos basados en Fortran,
La relación molar de óxido de etileno a agua es 1:12.
El óxido de etileno y el agua entran al reactor a 25 ◦C y 3 MPa. La escriba la condición NH2O =12 NETO (Figura 5.6).
5 Reactores 103
Figura 5.1 Diagrama de flujo del modelo de reactor de conversión en Aspen Plus
5 Reactores 105
Para comprobar la carga térmica del reactor, consulte Resultados Vuelva a ejecutar la simulación y compruebe los resultados en la
en el bloque REACTOR; En este proceso se requieren 2.562,5 página ResultadosReacciones. El calor de reacción de la
kW de calor. hidratación del óxido de etileno es de alrededor de −95 kJ mol−1
Para calcular el calor de reacción, esta opción debe activarse en y el de la reacción lateral de etilenglicol a dietilenglicol es de −96
la página ConfiguraciónCalor de reacción siguiendo los pasos kJ mol−1 (consulte la Figura 5.10). Estos resultados muestran
que se muestran en la Figura 5.9; introduzca la temperatura, la que ambas reacciones son exotérmicas. Sin embargo, el calor de
presión y la fase de referencia dadas. reacción no es suficiente para calentar el
reactivos a la temperatura de reacción; por lo tanto, se la reacción ocurre espontáneamente en la dirección que
debe introducir calor adicional desde el exterior al muestra la ecuación estequiométrica (a la derecha); si ΔrG
proceso. tiene un valor positivo, la reacción ocurre espontáneamente en
la dirección opuesta y si ΔrG =0, el sistema está en el estado
de equilibrio. El cambio de energía libre de Gibbs de un
5.3 Modelos de reactores de equilibrio químico sistema de reacción puede estar dado por
5 Reactores 107
donde ΔrG◦ es la energía libre de Gibbs estándar (referencia), ai La energía libre de Gibbs estándar se puede calcular a partir de
es la actividad del componente i, R es la constante de los gases
las energías libres de Gibbs estándar de formación Δf Gi ◦ de los
y T la temperatura (K). En un estado de equilibrio, ΔrG =0, y de componentes como
la ecuación (5.5), obtenemos ◦ ◦ (5.8)
ΔrG = ∑ iΔf Gi
ΔrG◦ (5.6) i
= −RT ln I
◦ , y estándar
Si el calor de formación estándar, Δf Hi, la
donde Ke es la constante de equilibrio de la reacción química: ◦
entropía absoluta, Si de ,los componentes están disponibles, ΔrG◦
se puede calcular como
(5.7) ◦ ◦ ◦ (5.9)
cuando = ∏
yo
un yo ΔrG = ∑ iΔf hola − T ∑ y Si
i i i
Machine Translated by Google
5.3.1 Modelo REquil de Aspen Plus El reactor de equilibrio calcula el equilibrio químico y de fase
simultáneamente; por lo tanto, se deben conectar al modelo al menos dos
El modelo REquil de Aspen Plus calcula las tasas de flujo de la corriente
flujos de material de salida. Si no se conecta ninguna corriente de energía,
del producto utilizando constantes de equilibrio determinadas a partir de
HYSYS asume un reactor adiabático; si se conecta una corriente de
la energía libre de Gibbs. Las constantes de equilibrio se basan en la
energía al modelo, se deben especificar las condiciones de salida o de
estequiometría de la reacción y la distribución del rendimiento
calor.
especificadas por el usuario. REquil calcula el equilibrio químico y de fase
simultáneamente. REquil se puede utilizar para modelar un reactor cuando
se conoce la estequiometría de la reacción y algunas o todas las
Ejemplo 5.2 Las siguientes dos reacciones ocurren en la producción de
reacciones alcanzan el equilibrio químico.
gas de síntesis a partir de gas natural.
Para cada reacción, puede seleccionar una especificación de equilibrio CH4 + H2O = 3H2 + CO (R5.3)
restringida, ya sea como la extensión molar de la reacción o la aproximación CO + H2O = H2 + CO2 (R5.4)
de la temperatura al equilibrio. Si no proporciona estas especificaciones,
REquil asume la reacción para alcanzar el equilibrio químico. Si La primera reacción es la reacción de reformado y la segunda es la
selecciona el enfoque de temperatura para una reacción, también puede reacción de cambio de agua a gas mediante la cual se produce aún más
proporcionar una estimación de la extensión molar para mejorar la hidrógeno y el CO se convierte en CO2 .
convergencia de los cálculos de equilibrio químico. Modele estas reacciones como reacciones de equilibrio usando Aspen
HYSYS. Para el cálculo de la constante de equilibrio, utilice el método de
REquil realiza cálculos instantáneos monofásicos o bifásicos anidados energía libre de Gibbs.
dentro de un bucle de equilibrio químico. REquil no puede realizar cálculos Para este cálculo, considere el gas natural como metano pobre que se
El reactor de equilibrio proporciona diferentes posibilidades para el temperatura del reactor entre 350 y 900 ◦C.
ecuación (Ke), ln (Ke ), que se supone que es una función únicamente defina también la temperatura de la corriente de salida del gas y se
de la temperatura, se determina a partir de la siguiente ecuación: puede utilizar el valor final del rango de temperatura (900 ◦C).
5 Reactores 109
fracción molar del componente) y Específicos de la variable La página de ReacciónResultados muestra que el valor de la
(componente), luego haga clic en Agregar para agregar el parámetro constante de equilibrio es prácticamente cero y el de la segunda
seleccionado al estudio de caso. Agregue gradualmente las reacción es de 22,6 a 350 ◦C. Sin embargo, a 900 ◦C, la constante de
fracciones molares de todos los componentes al caso de estudio equilibrio de la primera reacción aumenta a 1202, pero la de la
y también la temperatura de la corriente de gas de salida. Tenga segunda reacción es solo de 0,78, lo que indica que la primera
en cuenta que la variable será identificada por HYSYS como reacción tiene lugar a temperaturas más altas y la segunda reacción
Independiente solo si su valor es definido por el usuario. a temperaturas más bajas.
Después de seleccionar todos los parámetros a observar, cierre la Por lo tanto, un modelo más preciso de este proceso consiste en dos
página Navegador de variables. reactores, donde la primera reacción tiene lugar en el primer reactor y
En la página Configuración del estudio de caso, defina el rango de la la segunda en el segundo.
variable independiente y el tamaño del paso como se muestra en la Figura 5.14. reactor.
Ejecute el estudio de caso y verifique los resultados; los resultados
se presentan en forma de tabla y gráficos. La tabla de resultados
se puede copiar fácilmente a EXCEL para su posterior análisis. La 5.3.3 Modelo RGibbs de Aspen Plus y modelo de reactor Gibbs de
figura 5.15 muestra la composición del gas producido a temperaturas Aspen HYSYS
de 350 a 900 ◦C. Ambos modelos consideran la condición de que la energía libre de
Gibbs del sistema reaccionante sea mínima en equilibrio para calcular
la composición de la mezcla de productos. Por lo tanto, estos modelos
no necesitan la estequiometría de la reacción para calcular la
composición del producto y pueden utilizarse si se desconoce la
estequiometría de la reacción y se supone que se alcanza el equilibrio
químico.
RGibbs de Aspen Plus proporciona los siguientes tipos de cálculos:
3. restringir el equilibrio químico: especificar el enfoque de temperatura 5.4 Modelos de reactores cinéticos
o las reacciones, y
4. restrinja el equilibrio químico: especifique el trabajo y la temperatura, Si se utiliza un modelo de reactor cinético, el último término de la
y calcule el enfoque de temperatura. ecuación (5.1) se calcula utilizando las ecuaciones de velocidad de
las reacciones químicas. Básicamente, se utilizan dos tipos de
Una aplicación muy efectiva de RGibbs es el modelado de
expresiones cinéticas para calcular la velocidad de reacción:
combustión o gasificación de sólidos no convencionales como biomasa,
residuos y carbón (consulte el Capítulo 14). 1. expresión de ley de potencia y 2.
El reactor de Gibbs de Aspen HYSYS puede funcionar únicamente modelo generalizado de LangmuirHinshelwoodHougenWatson
como separador, un reactor que minimiza la energía libre de Gibbs sin (LHHW).
un conjunto de reacción adjunto o como un reactor que usa reacciones
La expresión de ley de potencia mostrada por la ecuación (5.13) se
de equilibrio. Cuando se adjunta un conjunto de reacciones, la
usa para modelar modelos de reactores cinéticos en Aspen Plus
estequiometría involucrada en las reacciones se usa en los cálculos
del reactor de Gibbs. (RCSTR, RPlug y RBatch) si la cinética de ley de potencia es
Figura 5.14 Especificación del rango de variable independiente para el caso de estudio
Machine Translated by Google
5 Reactores 111
0.70
0,60
0.50
CO
Fracción
molar
0.40 CO2
H2O
0.30
H2
0.20 CH4
0.10
0.00
300 400 500 600 700 800 900 1000
Temperatura (°C)
seleccionado. La misma expresión es utilizada por Aspen Las constantes de fuerza impulsora k1 y k2 y también la
HYSYS, si el tipo de reacción es seleccionado como Cinético. constante del término de adsorción Ki pueden depender de la
− 1
norte
temperatura. En Aspen Plus, estas dependencias de temperatura
− ( ES )( T1
T0 )∏
yo
_
i (5.13) están dadas por
es = k ( T
T0 )u yo=1
B ln(K) = A + + C ln(T) + DT (5.19)
donde k es el factor preexponencial, E la energía de activación, T T
la temperatura absoluta, T0 la temperatura de referencia, Ci la donde A, B, C y D son constantes.
concentración del componente i y es el exponente i del Para modelar reacciones de velocidad controlada en reactores
componente i.
continuos de tanque agitado (CSTR) en Aspen Plus, use el
Si no se especifica T0 , Aspen Plus utiliza la ley expres bloque de operación de unidad RCSTR; en Aspen HYSYS, utilice
sión en la siguiente forma reducida: el modelo de reactor continuo de tanque agitado. Para modelar
norte
CSTR
MEZCLA
ETANOL
F PAG
ÁCIDO A
dónde escriba e ingrese el valor dado para el volumen del reactor (Figura
5.17).
k = 1,206 × 106 e −E∕RT (5.21) En el caso de modelos de reactores cinéticos (RCSTR, RPlug,
Rbatch), las reacciones químicas primero deben definirse fuera de
E=54.240 kJ kmol−1 (6) y Ke =4,5 es la constante de equilibrio
la página de especificaciones del reactor y luego agregarse al
basada en la molaridad. modelo; para definir un conjunto de reacciones, seleccione
Reacción en el panel de navegación principal y siga los pasos que
b. encontrar la dependencia de la conversión de etanol en el volumen
se muestran en la Figura 5.18 para seleccionar un tipo de conjunto
del reactor en el rango de volumen del reactor de 0,5 a 6 m3
. de reacciones GENERAL.
Cree una reacción haciendo clic en Nuevo, aparece la página de
estequiometría de reacción; seleccione POWERLAW como la clase
Solución: de reacción e ingrese la estequiometría de reacción; seleccione
Abra Aspen Plus, seleccione una lista de componentes y el método esta reacción como una reacción reversible (Figura 5.19).
termodinámico apropiado (NRTLHOC en este ejemplo) como se
explica en los Capítulos 1 y 2. Ingrese los parámetros cinéticos de la expresión de la ley de
Cambie al entorno de simulación y prepare el diagrama de flujo del potencia como se muestra en la Figura 5.20; seleccione la fase
proceso (Figura 5.16). de reacción como Líquido, la molaridad como Base, el volumen del
reactor como Base −1 como Unidad de Tarifa.
Especifique las corrientes de entrada por temperatura, presión,
composición y flujo molar. de velocidad y kmol m−3 s Vaya a la página Equilibrio e ingrese
Especifique el bloque RCSTR; ingrese la presión y la temperatura; el valor de la constante de equilibrio siguiendo los pasos de la Figura 5.21.
para Asalto, seleccione LiquidOnly como la fase válida; use el Vuelva a la página de configuración de CSTR y agregue la reacción
volumen del reactor predeterminado como especificación definida al modelo de reactor como se muestra en la Figura 5.22.
5 Reactores 113
5 reactores 115
En la página de definición de entrada de sensibilidad, defina otras En la página de Fortran, escriba la fórmula para el cálculo de conversión
variables. Para calcular la conversión de etanol, debe conocer el flujo de etanol como se muestra en la Figura 5.25; tenga cuidado de usar los
molar de etanol en las corrientes de alimentación y producto; defina símbolos en la forma exacta en que están definidos en la página Definir.
también un parámetro local y llámelo CON; Para definir el flujo molar de
etanol, se debe especificar el tipo de parámetro, la corriente, la En la página Tabular, especifique la columna de la tabla donde se mostrará
subcorriente, el componente y la unidad, como se muestra en la Figura la conversión de etanol.
5.24. Para parámetros locales, es suficiente especificar solo el tipo como Ejecute la simulación y verifique los resultados del análisis de sensibilidad;
LocalParam. los resultados se presentan en forma de tabla.
Las curvas de resultados se pueden producir directamente en Aspen Plus,
o se pueden exportar a EXCEL para producir curvas de resultados. La velocidad de reacción y los parámetros cinéticos derivados de (7)
La dependencia de la conversión en el volumen del reactor se son
presenta en la Figura 5.26, y la gráfica se produjo en Aspen Plus
1
utilizando Curva de resultado de la barra de herramientas.
( pEB Cuando pST pH2 )
r1 = k1 (5.22)
(1 + KST pST)
r2 = k2 pEB (5.23)
Ejemplo 5.4 Se produce estireno por deshidrogenación catalítica
de etilbenceno (R5.6) a temperaturas a partir de 630 ◦C y presiones r3 = k3 pEB (5.24)
ligeramente superiores a la atmosférica en presencia de vapor de
( −E1
RT )
k1 = A1 e (5.25)
agua. En este proceso tienen lugar dos reacciones secundarias
principales: (R5.7) pirólisis de etilbenceno a benceno y etileno y
A1 =3524,4 kmol m−3 s −1 Pa−1 E1
, =158,6 kJ mol−1 ,
(R5.8) hidrodesalquilación de etilbenceno, donde se producen
tolueno y metano. k1
( −(E1−Ee ) RT )
k1− = (5.26)
Cuando
= ( A1Sí ) e
C6H5 C2H5 ↔ C6H5 C2H3 + H2 (R5.7) (R5.6) ln (Ae )=27,16 Pa, Ee =124,26 kJ mol−1 ,
C6H5 C2H5 → C6H6 + C2H4 C6H5
( −E2
RT )
C2H5 + H2 → C6H5 CH3 + CH4 (R5.8) k2 = A2 mi (5.27)
60
55
50
45
40
Conversión
(%)
35
30
25
20
15
1000 1250 1500 1750 2000 2250 2500 2750 3000 3250 3500 3750 4000 4250 4500 4750 5000 5250 5500 5750 6000
Volumen rector (L)
Figura 5.26 Resultados del análisis de sensibilidad, conversión versus volumen del reactor
Machine Translated by Google
5 Reactores 117
En la pestaña Tasa de reacción, ingrese los parámetros Especifique la corriente de alimentación por temperatura (630
cinéticos para las reacciones directa e inversa, y los parámetros ◦C), presión (137,8 kPa), flujo molar total (165 kmol h y −1 ),
−1
del término de adsorción como se muestra en la Figura 5.27; composición (15 kmol h de etilbenceno y 150 kmol h de
agua). −1
consulte la ayuda cinética para comprender correctamente la
definición de parámetros cinéticos.
Defina la estequiometría y los parámetros cinéticos de ambas
reacciones secundarias de manera similar; sin embargo, no se
consideran reacciones reversibles, por lo que solo debe ingresar
parámetros cinéticos para reacciones directas; recuerde que el
orden del hidrógeno en la ecuación de velocidad, r3 , es
cero.
Agregue un conjunto de reacciones al paquete de fluidos y
muévase al entorno de simulación.
Instale un PFR, el diagrama de flujo del proceso se muestra
en la Figura 5.28; dado que se va a modelar un reactor
adiabático, no conecte la corriente de energía.
Figura 5.30 Dimensionamiento de PFR en Aspen HYSYS Figura 5.31 Selección de un bloque de hoja de cálculo
Machine Translated by Google
5 Reactores 119
En la pestaña Parámetros, en Diseño, seleccione la ecuación de simulación, calculando otros parámetros para ser utilizados en casos
Ergun para el cálculo de la caída de presión. de estudio o exportándolos en cualquier posición en su simulación.
Vaya a la página de reacciones y agregue el conjunto de reacciones
definido al modelo de reactor. Siguiendo los pasos que se muestran en la Figura 5.32, agregue
En la misma página, defina la información de integración y los datos todas las variables como variables de importación: flujo molar de
del catalizador como se muestra en la Figura 5.29, ingrese el diámetro etilbenceno en la alimentación, flujo molar de etilbenceno en la
de las partículas (5 mm), la esfericidad de las partículas (1) y la corriente del producto y flujo molar de estireno en la corriente del producto.
densidad sólida (2500 kg m−3 ) . Vaya a la página Hoja de cálculo (Figura 5.33) y escriba la fórmula
En la pestaña Tamaño, en Clasificación, defina las dimensiones del para el cálculo de selectividad en una celda de selección.
tubo y los parámetros de empaque del tubo como se muestra en la La selectividad de la reacción de etilbenceno a estireno puede
Figura 5.30. definirse como el flujo molar de estireno producido por flujo molar de
Dado que HYSYS no calcula la selectividad directamente, cree un etilbenceno reaccionado.
bloque de hoja de cálculo para calcular la selectividad al estireno. Para mostrar los perfiles de temperatura, presión y composición,
Para hacer esto, seleccione Hoja de cálculo de la paleta modelo utilice la pestaña Rendimiento. Los resultados se presentan en
como se muestra en la Figura 5.31. La hoja de cálculo permite forma de tablas y gráficos. La Figura 5.34 muestra los perfiles de
importar cualquier parámetro de su temperatura y presión en PFR. el reactor
Machine Translated by Google
la temperatura disminuye en condiciones adiabáticas de 630 a 536 Para trazar los perfiles de conversión y selectividad, se debe usar un
◦C, por lo que el proceso es endotérmico. La caída de presión total estudio de caso (consulte el Ejemplo 5.2 para obtener detalles sobre
calculada por la ecuación de Ergun es de alrededor de 16 kPa. el uso de estudios de caso en HYSYS).
Agregue la longitud del reactor como variable independiente, la
La Figura 5.35 muestra el perfil de composición a lo largo del PFR; el conversión real en el bloque PFR y la selectividad definida en el
flujo molar de etilbenceno disminuye de 15 lo que representa la bloque de la hoja de cálculo como parámetros observados; elija los
−1
a 8 kmol h , conversión total de etilbenceno del 53,5%. El flujo bordes de la longitud del reactor y el paso para el estudio de caso.
molar de hidrógeno está ligeramente por debajo del flujo molar de
estireno, porque una porción de hidrógeno reacciona en (R5.8). El Ejecute el estudio de caso y grafique los resultados tanto para la
flujo molar de los subproductos (tolueno, benceno, etileno y metano) conversión de etilbenceno como para la selectividad a estireno (Figura
es generalmente mucho menor que el de los productos principales; 5.36); la conversión total de etilbenceno alcanza valores de alrededor
por lo tanto, podemos esperar una buena selectividad de la reacción del 54%, que es una conversión cercana al equilibrio y el aumento
principal. Como resulta del perfil de composición, la velocidad de la del volumen del reactor no afecta significativamente la conversión; sin
segunda reacción secundaria (R5.8) es más alta que la de (R5.7). embargo, la selectividad al estireno disminuye linealmente con la
longitud del reactor desde 97% a 95,85%.
Temperatura Presión
600,0 138,0
Líquido Presión
590,0
Fluido utilitario 136,0
134,0
580,0
132,0
570,0
Temperatura
(ºC)
Presión
(kPa)
130,0
560,0
128,0
550,0
126,0
540,0 124,0
530,0 122,0
0,000 0,5000 1,000 1,500 2,000 2500 3000 3500 4000 0,000 0,5000 1,000 1,500 2,000 2500 3000 3500 4000
Longitud del reactor (m) Longitud del reactor (m)
Figura 5.34 Perfiles de temperatura y presión en PFR para la producción de estireno, respectivamente
Machine Translated by Google
5 Reactores 121
54 97
52 96.8
50 96,6
Conversión
(%)
48 96.4 Selectividad
estireno
del
46 96.2
44 96
40 95.6
0123456
Longitud del reactor (m)
Machine Translated by Google
5.5 Selección y costeo de reactores Los reactores por lotes se utilizan generalmente para la producción a
Figura 5.37 Lista de agitadores y tanques agitados disponibles en APEA para el mapeo de CSTR
Machine Translated by Google
5 Reactores 123
reactor, reactor de lecho fluidizado de lecho burbujeante, reactor de Costo de los reactores CSTR y, por lo tanto, es el predeterminado para
tanque rociado, etc. (810). el mapeo CSTR. El PFR se puede mapear como una torre empacada,
Para estimar el costo de un reactor por el Aspen Process Economic un intercambiador de calor de un solo tubo, un intercambiador de calor
Analyzer, se debe seleccionar un dispositivo apropiado de la lista de carcasa y tubos, o como un horno tubular.
disponible. Un reactor CSTR se puede mapear como diferentes tipos de
tanques agitados o agitadores. En algunos casos, también se pueden
mapear como batidoras o amasadoras. La figura 5.37 muestra la lista de Ejemplo 5.5 Seleccione un dispositivo apropiado para el mapeo de CSTR
agitadores y tanques agitados disponibles que se pueden usar para en el proceso de acetato de etilo (Ejemplo 5.3).
calcular el costo del reactor CSTR. Tanque agitado, cerrado, encamisado Determine el costo instalado del reactor y el costo de los servicios
es el dispositivo apropiado elegido con más frecuencia para públicos en función del volumen del reactor y la conversión,
respectivamente.
400000 400000
350000 350000
300000 Costo instalado 300000 costo del equipo
250000 costo del equipo 250000
Costo instalado
(US$)
Costo
200000 (US$)
Costo 200000
150000
150000
100000
100000
50000
50000
0 35 55 75 95 0
0123456
Conversión (%)
Volumen del reactor (m3 )
Figura 5.39 Relación entre el volumen del reactor, la conversión y el costo del equipo
Machine Translated by Google
60000 60000
55000 55000
50000 50000
servicios
públicos
(US$)
Costo
los
de
45000 servicios
públicos
(US$)
Costo
los
de
45000
40000 40000
35000 35000
30000 30000
30 50 70 90 0123456
Conversión (%) Volumen del reactor (m3 )
Figura 5.40 Costo de los servicios públicos en función del volumen del reactor y la conversión
Mapee el reactor como cerrado por un tanque agitado, encamisado como se Repita los dos pasos anteriores para volúmenes de reactor de 0,5 a 5 m3 por
Siga los mismos pasos que en los Ejemplos 3.5 y 4.4 para dimensionar y El costo del equipo y el costo instalado para diferentes volúmenes de reactores
evaluar el equipo. (diferentes conversiones) se muestran en la figura 5.39. Tanto el costo del
Como tipo de material predeterminado, se establece acero al carbono (CS). equipo como el costo de instalación muestran un rápido aumento para las
Sin embargo, la presencia de ácido sulfúrico como catalizador y ácido acético conversiones por encima del 65%; la razón es que por encima de estas
como reactivo presenta un potencial de corrosión, por lo que se recomienda el conversiones, el cambio de conversión con el volumen del reactor es mucho
acero inoxidable 304 (SS304) como tipo de material. Cambie el tipo de material más lento (ver Figura 5.26).
Cambie el tamaño y luego vuelva a evaluar el proceso con el nuevo tipo de la conversión del reactor. Con el aumento del volumen del reactor hasta 3,5
material. m3 por encima de este volumen , aumenta el costo de los servicios públicos;
Registre el costo del equipo, el costo instalado y el costo de los servicios del reactor, el coste de los servicios públicos muestra un valor constante.
públicos.
En la pestaña Entrada del bloque CSTR, cambie el valor del volumen del
reactor y vuelva a ejecutar la simulación.
Referencias
1 Levenspiel O. Ingeniería de reacciones químicas, 3.ª ed. deshidrogenación catalítica de etilbenceno en reactores de membrana
Nueva York: John Wiley & Sons; 1999. compuesta de paladio. química Ing. ciencia
2 Missen RW, Mims CA, Saville BA. Introducción a 1999;54(10):1431–1439.
Ingeniería de Reacciones Químicas y Cinética. Nueva York: John 8 Towler G, Sinnott R. Diseño de ingeniería química,
Wiley & Sons, Inc.; 1999. Principio, práctica y economía del diseño de plantas y procesos, 2ª
3 Ayuda de Aspen Plus®V9. Burlington, MA: Aspen ed. Ámsterdam: Elsevier; 2013.
Technology, Inc.; 2016. www.aspentech.com. 9 Peters MS, Timmerhaus KD, West RE. Diseño y economía de
4 Schefflan R. Enséñese usted mismo lo básico de Aspen Plus. plantas para ingenieros químicos, 5.ª ed. Nueva York: McGraw
Hoboken, Nueva Jersey: John Wiley & Sons, Inc.; 2011. Hill; 2004.
5 Ayuda de Aspen HYSYS®V9. Burlington, MA: Aspen 10 Couper JR, Penney WR, Fair JR, Walas SM. Selección y diseño de
Technology, Inc.; 2016. www.aspentech.com. equipos de proceso químico, 3.ª ed.
´ ´
6 Ilavsky J, Valt ´ yni J, Brunovsk a A, Surov y J. ´ Cinética química Ámsterdam: Elsevier; 2010.
´ ´
aplicada´ y teoría de reactores I . Bratislava: Alfa; 1990. 11 Zimmerman CC, York R. Desmetilación térmica de tolueno. Ing.
Ind. química Proceso Des. desarrollo 1964;3(3):254–
¨
7 Dittmeyer R, Hollein V, Quicker P, Emig G, Hausinger G, Schmidt F. 258.
Factores que controlan el desempeño de
Machine Translated by Google
125
Equipos de separación
La mayoría de los equipos de una fábrica de productos químicos se utilizan Muchos modelos de operación de unidades de separación en Aspen Plus,
para la separación. En un caso general, una mezcla heterogénea del reactor como el modelo RadFrac que se usa con mucha frecuencia, permiten
se somete primero a la separación de la fase sólida utilizando equipos cálculos trifásicos (vaporlíquidolíquido). En caso de existir dos fases
como ciclones, centrífugas, filtros, etc. Luego, la mezcla homogénea líquidas en Aspen HYSYS, se deben utilizar diferentes modelos de operación
obtenida se somete a más operaciones de separación para proporcionar unitaria. En este capítulo, nos enfocamos en los modelos de operación de
los productos requeridos. Como se muestra en la Figura 6.1, se pueden unidad de contacto único y contacto múltiple más utilizados implementados
aplicar diferentes métodos para separar una mezcla homogénea; la más tanto en Aspen Plus como en Aspen HYSYS.
utilizada es la separación por creación o adición de una nueva fase. La
destilación, absorción, extracción, destilación extractiva y azeotrópica,
desorción, cristalización, secado, sublimación y evaporación pertenecen a
este grupo de procesos de separación. 6.1 Separación de fase de contacto único
Tanto los modelos de separación de un solo contacto como los de La ecuación (6.1) representa el balance de materia del componente i, la
contacto múltiple que aplican la creación o adición de una nueva fase se ecuación (6.2) el equilibrio líquidovapor y la ecuación (6.4) el balance de
implementan tanto en Aspen Plus como en Aspen HYSYS. entalpía del proceso. El
Diseño y simulación de procesos químicos: aplicaciones Aspen Plus y Aspen HYSYS, primera edición. Juma Haydary. © 2019
John Wiley & Sons, Inc. Publicado en 2019 por John Wiley & Sons, Inc.
Sitio web complementario: www.wiley.com/go/Haydary/ChemDesignSimulation Aspen
Machine Translated by Google
Separación
productos
la constante de equilibrio se calcula mediante el modelo modelo. El FLASH3 de Aspen Plus y el separador trifásico de Aspen
termodinámico elegido (ver Sección 2.2). Las entalpías de las fases HYSYS se utilizan para la separación vaporlíquidolíquido. Para la
de alimentación, líquido y vapor se calculan utilizando relaciones de separación líquidolíquido en una sola etapa, use DECANTER en
correlación para parámetros dependientes de la temperatura de Aspen Plus o ThreePhase Separator en Aspen HYSYS.
componentes puros.
El número total de variables en el esquema del proceso (Figura
6.2) es 3k + 13, donde k representa el número de componentes. El
número total de ecuaciones es 2k + 9, por lo que hay k + 4 grados Ejemplo 6.1 Se precalientan 100 −1
kmol h de una mezcla que
de libertad. Para un problema estándar, se conocen nF , xFi (para contiene 10, 20, 30 y 40 mol% de propano, nbutano, npentano y n
i−1 componentes; para el último componente, se calcula a partir de hexano, respectivamente, antes de ingresar a una torre de
la ecuación de suma), TF y PF ; el número restante de grados de destilación . La fracción molar líquida de la mezcla es, por lo tanto,
libertad es 2. Al definir estos dos parámetros adicionales, se del 80%. Calcule la composición de las fases líquida y vapor y la
especifica el sistema. Se pueden usar diferentes combinaciones de temperatura de la mezcla líquidovapor a 700 kPa. La temperatura
EN
proceso utilizando un modelo de calentador y un modelo de
TV , fotovoltaica , alto voltaje separador; como se muestra en la Figura 6.3. Conecte una
corriente de energía solo al calentador. Si la corriente de energía
del bloque separador no está conectada, funcionará como un
nF , XFi
separador adiabático.
F
Ingrese las condiciones de la corriente de entrada (F0) (25 ◦C.
TF , FP , hF −1 y 100 kmol h 700 kPa, ) y la composición; establezca
también la fracción de vapor (0,2) de la corriente de salida del
q nL , Xi calentador (F) y la presión de la corriente F (700 kPa).
q L
Cuando se ingresan los parámetros dados, el esquema de color
TL , PL , hL
de HYSYS indica que el proceso está calculado y puede verificar
Figura 6.2 Esquema de una separación continua de líquidovapor de los resultados en la página Hoja de trabajo.
una sola etapa
Machine Translated by Google
Nombre F0 F EN L
FLASH3
6.2.1 Método de destilación abreviado
S1 L1 L2 xiD xHKD
Componentes GRAMO
control de calidad
100(N − Nmín) F1
norte + 1 = f ( 100(R −RRmín)
+ 1 ) (6.11)
Qw
La forma gráfica y numérica de esta correlación ha sido publicada
T100
en varias publicaciones de ingeniería química;
Nn Bueno mi Figura 6.6 Diagrama de flujo del modelo de acceso directo en Aspen HYSYS
(6.12)
Nuevo Méjico Nm, min
Solución:
registro nn (6.13)
Nuevo Méjico
= 0.206 log [ nB nD
( xHK xLK ) F ( xLKB
xHKD )2 ] Continúe en el Ejemplo 6.1 y defina un nuevo flujo (F1) con los
mismos parámetros que el flujo F; para hacer esto, ingrese la
En estas ecuaciones, Nmin representa el número mínimo de etapas, secuencia F1 y elija Definir desde otra secuencia; luego seleccione
Rmin relación mínima de reflujo, N número real de etapas, R relación flujo F y OK.
de reflujo externo, = ki∕kj es la volatilidad relativa del
yo componente i al
En la Paleta del modelo de columna, seleccione el modelo de
componente j, y q representa la cantidad de líquido agregado a la destilación de acceso directo.
etapa de alimentación ingresando una unidad de alimentación. Los Conecte la corriente F1 como la corriente de entrada y defina dos
subíndices LK, HK, F, D, B, n y m representan la llave ligera, la llave corrientes de material: D y W, y dos corrientes de energía: servicio
pesada, la alimentación, el destilado, el fondo, la sección de del condensador, Qc, y servicio del rehervidor, Qw. El diagrama
rectificación de columna y la sección de extracción de columna, de flujo del proceso se muestra en la Figura 6.6.
respectivamente. En la pestaña Parámetros, en Diseño, complete los requisitos para
El método de destilación abreviado proporciona solo resultados la distribución de componentes clave y la presión de la columna;
aproximados; calcula con el valor medio de volatilidad relativa en la cuando se especifican estos parámetros, HYSYS calcula la relación
columna. Por esta razón, los resultados son confiables solo para mínima de reflujo. Calcule el valor de la relación de reflujo externo
sistemas ideales. Para calcular las columnas de destilación por el multiplicando este valor por el coeficiente dado e ingréselo como
método abreviado, use el bloque de operación de la unidad DSTWU valor de la relación de reflujo externo (paso 4 en la Figura 6.7).
en Aspen Plus y el modelo de destilación abreviado en Aspen
HYSYS. En la página Rendimiento, verifique los parámetros de la columna
de destilación; en la pestaña Composición, en Hoja de trabajo,
verifique la composición de los productos.
Ejemplo 6.3 Una mezcla precalentada del Ejemplo 6.1 (corriente F)
debe separarse en una torre de destilación.
La fracción molar de npropano en el producto de fondo es 0,05 y la Para la separación de la mezcla de hidrocarburos dada, como el
de nbutano en el destilado también es 0,05. destilado contiene 95% mol de propano y el fondo contiene 95% mol
Utilizando el método de FenskeUnderwoodGilliland, calcule la de nbutano, el valor de la relación mínima de reflujo es 1,54 y el
relación mínima de reflujo y el número mínimo de etapas teóricas número mínimo de platos es 5,65. Si R = 1,5 Rmin, el número real de
necesarias. Suponiendo una relación de reflujo 1,5 veces mayor que bandejas es 10,74, la etapa de alimentación óptima es la séptima
la relación de reflujo mínima, calcule el número real de bandejas, la desde arriba, la temperatura del condensador es de alrededor de 56
etapa de alimentación óptima, la temperatura del condensador y del ◦C y la temperatura del rehervidor es de alrededor de 122 ◦C. El
rehervidor, y la distribución de agua no clave. destilado contiene 35,92, 58,99,
Machine Translated by Google
es la eficiencia de la etapa j.
nF2 Q2 k
2
Sy,j = ∑ yi,j − 1 = 0 (6.21)
nSL2 yo=1
nV3
nL2
nSV3 Ecuación H:
como método predeterminado. Estos modelos de operación unitaria HYSYS. Cuando instala el modelo de diagrama de flujo secundario
también permiten el uso del método NewtonRaphson SC y otros de columna de destilación, HYSYS crea un diagrama de flujo
algoritmos. secundario de columna que contiene todas las operaciones y flujos
Los algoritmos de adentro hacia afuera consisten en dos bucles de asociados con la plantilla de columna que ha elegido. El diagrama
iteración anidados. En el bucle interno, las ecuaciones MESH se de flujo secundario funciona como una operación unitaria en el
resuelven utilizando un conjunto aproximado de parámetros diagrama de flujo principal. Un diagrama de flujo secundario para la
termodinámicos. En el bucle exterior, se emplean modelos columna de destilación proporciona una serie de ventajas:
termodinámicos exactos para actualizar los parámetros de las aislamiento del solucionador de columna, uso opcional de diferentes
ecuaciones empíricas utilizadas en el bucle interior. El algoritmo de paquetes de propiedades, construcción de plantillas personalizadas
adentro hacia afuera se puede emplear en la solución de todas las y la capacidad de resolver varias torres simultáneamente.
operaciones multietapa y multicomponente, como destilación, destilación Continúe en el Ejemplo 6.1; pero en lugar del separador, instale un
azeotrópica, destilación extractiva, absorción, extracción por desorción, diagrama de flujo secundario de la columna de destilación.
etc. Inicie la especificación de la columna haciendo doble clic en el
modelo de operación de la unidad. Ingrese los parámetros en varios
Ejemplo 6.4 Usando el método riguroso de adentro hacia afuera, calcule pasos usando el llamado experto de entrada de columna de
la dependencia entre el número de etapas teóricas y la relación de destilación. Cuando se especifican los parámetros mínimos
reflujo, si la mezcla precalentada del Ejemplo 2.1 (corriente F) se separa requeridos, se activa el botón Siguiente y se pueden ingresar los
en una torre de destilación, el destilado contiene 99 mol% de parámetros requeridos en la siguiente pestaña. En la primera pestaña
componentes ligeros (propano y nbutano) y al mismo tiempo la (Figura 6.10), defina el número de etapas teóricas (como valor inicial,
recuperación común de propano y nbutano es superior al 99%. Utilice puede usar el calculado por el método de destilación abreviado), la
el tipo de condensador de reflujo total y desprecie la caída de presión etapa de alimentación teórica y el tipo de condensador. Conecte el
en la columna y los intercambiadores de calor. flujo de material de entrada F y defina los flujos de material de salida
D y W, así como los flujos de energía Qc y Qw. Después de completar
estos parámetros, se activa el botón Siguiente y puede continuar a
Solución:
la página siguiente haciendo clic en Siguiente.
El diagrama de flujo secundario de la columna de destilación se
utilizará para el cálculo riguroso de la destilación en Aspen
La especificación de la columna por la relación de reflujo y la tasa En la página de especificaciones, ingrese el requisito de pureza del
de destilación no pueden asegurar la pureza y recuperación destilado siguiendo los pasos que se muestran en la Figura 6.14.
requeridas del producto; Para obtener la pureza del producto
requerida y la recuperación de los componentes ligeros, se deben Agregar una nueva especificación para la recuperación de
definir nuevas especificaciones. componentes de columna; defina la especificación de la columna
Se pueden usar dos métodos para definir la especificación de usando los mismos pasos que antes.
la nueva columna: usar la pestaña Especificaciones en Diseño o En la página Monitor, desactive las especificaciones originales
agregar especificaciones directamente desde la página Supervisar. (Relación de reflujo y Tasa de vap de Ovhd) y active la nueva
Para agregar una nueva especificación, seleccione Agregar especificación definida (Fracción de compensación y Recuperación
especificación en la página Monitor y luego elija el tipo de de compensación), HYSYS recalcula la columna.
especificación de la lista (Figura 6.13). Si el cálculo converge, verifique la composición
14
12
R 10
0
5 10 15 20 25 30 35
norte
de los productos utilizando la pestaña Hoja de trabajo o Libro de datos. donde es la volatilidad relativa original del componente A al componente
AB
La Tabla 6.3 muestra las condiciones de balance de masa y la B y del componente A al la volatilidad relativa de com
ABC
composición de las corrientes de alimentación y producto cuando el componente B después de la adición del solvente C.
número de etapas es 15 y se calculó el valor de 4.307 para la relación El punto de ebullición del nuevo solvente es más alto que el de los
de reflujo (Figura 6.15). componentes de la mezcla, y generalmente se alimenta a la parte superior
En la página Rendimiento, analice los resultados del cálculo de la de la columna y se obtiene en el producto de fondo.
columna, incluidos todos los perfiles de temperatura, presión, En caso de destilación azeotrópica, el nuevo solvente debilita los
composición y valor K. enlaces de hidrógeno y crea un nuevo azeótropo homogéneo o
heterogéneo con uno de los componentes de la mezcla original. El
Las nuevas especificaciones no permiten cambiar la composición del
solvente se alimenta generalmente a la parte inferior de la columna y se
producto cambiando el número de etapas; solo se puede cambiar la
obtiene en el destilado.
relación de reflujo y, en un cierto rango, también la tasa de destilado. Por
lo tanto, estas especificaciones se utilizan para analizar la dependencia
Tanto para la destilación extractiva como azeotrópica, se pueden
entre el número de etapas teóricas y la relación de reflujo.
aplicar métodos rigurosos de adentro hacia afuera o SC para la solución
de ecuaciones MESH. El modelo de operación de la unidad RadFrac de
Para obtener la dependencia entre N y R, repita el cálculo para Aspen Plus y el diagrama de flujo secundario de la columna de destilación
diferentes valores de N (cambiando proporcionalmente a la posición de Aspen HYSYS se pueden usar para modelar la destilación extractiva
de la corriente de alimentación) en el rango de Nmin (puede obtenerse y azeotrópica.
mediante el modelo de destilación abreviado) a N∞ correspondiente
a R = ∞ y Rmin, respectivamente.
−1
Ejemplo 6.5 Se separan 3,5 kg s de una mezcla binaria equimolar de n
Lea el valor calculado de la relación de reflujo para cada cálculo en la heptano y tolueno mediante destilación extractiva a presión atmosférica
página Monitor. utilizando n metilpirrolidona (NMP) como disolvente selectivo. Se utiliza
Dibuje el gráfico N versus R (Figura 6.16). una columna con capacidad de separación de 10 etapas teóricas, un
condensador total y un hervidor de agua. La alimentación del punto de
ebullición con una presión de 110 kPa entra en la octava etapa desde
6.3 Destilación azeotrópica y extractiva arriba comenzando por un condensador. El disolvente selectivo elimina el
tolueno y se obtienen en el producto de fondo. Por lo tanto, el nheptano
La destilación azeotrópica y extractiva se utilizan para separar mezclas debe recuperarse en el destilado (la velocidad de flujo del destilado es
con baja volatilidad relativa o mezclas azeotrópicas. En caso de igual al flujo molar de nheptano en la alimentación). La relación de
destilación extractiva, se alimenta a la columna un nuevo disolvente que reflujo es 3.
crea enlaces de hidrógeno con un componente de la mezcla original. La
selección de un solvente adecuado es un tema crucial del proceso de
El tolueno se separa del NMP en una segunda columna de destilación
destilación extractiva. atmosférica con seis etapas teóricas, donde la relación de reflujo es 2,5 y
Se deben considerar diferentes criterios; el más importante de ellos es el la tasa de destilado es igual a la tasa de flujo molar de tolueno en la
factor de selectividad, : alimentación. La alimentación entra en la etapa media de la columna.
C1 C2
NMP TOLUENO
NHEPTAN
W1
W2
el flujo molar es 630 kmol h −1 , su temperatura es de 100 ◦C, y la caída de presión de la columna, pero en este caso no conocemos
presión de 110 kPa, y entra en la etapa óptima; encontrar la etapa de las caídas de presión individuales, por lo que es suficiente ingresar
alimentación óptima para NMP considerando la pureza máxima de n la presión de la etapa superior/condensador.
heptano en el producto destilado de la primera columna. Encuentre la Especifique la columna C2 de manera similar a la columna C1.
dependencia de la pureza de los productos (nheptano y tolueno) y el Ejecute la simulación y verifique los resultados en las pestañas
requisito específico de NMP (kmol de NMP por kmol de alimentación). Results y Stream Results de cada bloque de operación unitaria, C1
y C2, respectivamente. El flujo molar y la composición de todas las
corrientes se muestran en la Tabla 6.4.
Solución:
Defina un bloque de sensibilidad (para obtener detalles sobre la
Inicie Aspen Plus y abra un nuevo caso con Blank Simulation. definición de un bloque de sensibilidad, consulte el Ejemplo 5.3);
como parámetro variable, use FEED STAGE para NMP y deje que
Cree la lista de componentes para esta simulación. Aspen lo cambie de 2 a 8; como parámetro observado, use la
Utilice el método de propiedad NRTL; en la pestaña de parámetros fracción molar de nheptano en la corriente destilada de C1.
de Interacciones binarias de NRTL, deje que Aspen estime los
parámetros que faltan por UNIFAC. Vuelva a ejecutar la simulación y trace los resultados del bloque de
Vaya al entorno de simulación y prepare un diagrama de flujo del sensibilidad; La Figura 6.20 muestra la dependencia de la pureza del
proceso como se muestra en la Figura 6.17 utilizando el bloque de nheptano en la etapa de alimentación de NMP. Si se alimenta cerca
operación de la unidad Radfrac para la columna de extracción (C1) de la parte superior de la columna, una porción del solvente se
y la columna de regeneración (C2). arrastra en el destilado y causa su impureza, por otro lado, si se
Especifique la corriente de alimentación por fracción de vapor (0), alimenta demasiado lejos de la etapa superior, la eficiencia de
−1 separación de las etapas superiores disminuye. y la impureza es
presión (110 kPa), caudal másico (3,3 kg s ) y fracciones molares
(xC7 = 0,5, xT = 0,5). causada por la presencia de tolueno. Para el proceso de destilación
Especifique la corriente de NMP por temperatura (100 ◦C), presión extractiva, es importante encontrar la etapa de alimentación óptima
(110 kPa), flujo molar (630 kmol h fracción molar de −1 ), y para el solvente.
NMP (xNMP = 1). En su bloque de sensibilidad, defina una nueva variable en la página
Especifique la primera columna: C1, Figura 6.18, seleccione Equi Vary (flujo molar de la corriente NMP (Stream VarNMPMole Flow)
librium como método de cálculo, ingrese el número de etapas y establezca su rango de 0 a lo que corresponde con el requisito
(incluidos el condensador y el rehervidor 10 + 2 = 12), y elija el tipo 700 kmol h −1 , específico de NMP de 0 a 5.34; en la pestaña
de condensador total y Hervidor como tipo de rehervidor. Definir, defina también el flujo molar de tolueno en el destilado de la
segunda columna C2.
Introduzca los valores de la tasa de destilado y la relación de
reflujo. Ejecute la simulación con esta nueva definición y represente
En la página Corrientes, especifique las etapas de alimentación de las gráficamente los resultados de la pureza molar tanto de nheptano
corrientes F y NMP como se muestra en la Figura 6.19. Utilice la como de tolueno frente a nNMP/nF (Figura 6.21). El aumento de
convención predeterminada (por encima del escenario). nNMP/nF da como resultado una mayor pureza de nheptano, pero
En la página Presión, ingrese la presión de la columna para la etapa por encima de nNMP/nF = 2, este aumento no es muy rápido y la
superior/presión del condensador (Figura 6.19). Si se conoce la pureza del tolueno comienza a disminuir. Significa que, si se utiliza
caída de presión del condensador o la presión en la etapa 2, se nNMP/nF > 2 para lograr una mejor pureza del nheptano, la columna
puede ingresar; La caída de presión en la columna se puede ingresar de regeneración C2 debería tener más de seis etapas teóricas.
como caída de presión de etapa individual o como caída de presión total.
Machine Translated by Google
C2 C1
Componente fracción molar (nheptano) 0.50000 0.00000 0.99581 0.00041 0.00430 0.00000
Componente fracción molar (tolueno) 0.50000 0.00000 0.00056 0.09413 0.97436 0.00262
Fracción molar componente (NMP) 0.00000 1.00000 0.00363 0.90546 0.02134 0.99738
−1 131.01 630 65.5 695.52 65.5 630.02
Flujo molar (kmol h )
−1
Caudal másico (kg h ) 12,600 62.453,5 6.562,81 68.490,72 6.047,25 62.443,48
0.995
0.99
0.985
heptano
Fracción
molar
de
n
0.98
0.975
0.97
0.965
0,96
0123456789
Etapa de alimentación de NMP
0,95
0.9
heptano
Fracción
molar
de
n
0.85
0.8
0.75
nheptano tolueno
0.7
0123456
Requisito específico de NMP (kmol NMP/kmol alimentación)
Machine Translated by Google
81 85
80.5 80
80
75
79.5
70
temperatura
(°C)
79 temperatura
(°C)
78.5
sesenta y cinco
78 60
77.5 55
77 50
0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1
x,y benceno x,y Acetona
Figura 6.22 Diagramas txy binarios isobáricos de mezclas binarias de benceno/ciclohexano y acetona/ciclohexano
su composición. Utilizando el método termodinámico para obtener benceno puro en el producto de fondo con alta
UNIQUAC, composición y punto de ebullición de la recuperación, incluso
CICLO01(80.78°C)
Figura 6.23 Mapa ternario del sistema de
acetona, benceno y ciclohexano
0.9
0.8
0.7
0.6
0.5
77.54°C
0.4
0.3
53.59°C
0.2
0.1
ACETONA BENCENO
(56,14°C) 0.1 0,2 0,3 0,4 0,5 0,6 0,7 0.8 0.9 (80.13°C)
Machine Translated by Google
Acetona
D Defina dos nuevas especificaciones (fracción molar de benceno en
el producto de fondo y recuperación de benceno en el producto de
fondo) utilizando valores máximos cercanos a 1 tanto para la
Alimentar recuperación como para la pureza.
ESO Reinicie el cálculo de la columna; si no converge en unos segundos
después de hacer clic en Ejecutar, detenga el cálculo y reduzca
el requisito de recuperación y pureza del benceno. Los resultados
EN obtenidos para la pureza del benceno de 0,995 y la recuperación
T100 de benceno de 0,99 se muestran en la Tabla 6.5; la relación de
reflujo para estas especificaciones fue de 6,6 y la tasa de destilado
Figura 6.24 Esquema de una columna de destilación azeotrópica en Aspen fue de 107,64 kmol h. La columna también puede converger −1 . para
HYSYS
valores más altos de pureza del benceno pero con valores altos
se debe agregar más acetona. Use el flujo de mol −1 de acetona de la relación de reflujo y el número de etapas teóricas. Se puede
de 85 kmol h en esta simulación. alcanzar una pureza molar de benceno de 0,999 y una recuperación
Inicie Aspen HYSYS, cree una lista de componentes y seleccione molar de benceno de 0,999 en una columna con 50 etapas teóricas
el paquete de fluidos UNIQAQ. con una relación de reflujo de 5,4, etapa de alimentación de mezcla
Vaya al entorno de simulación y defina dos flujos de material: para azeotrópica 15 y etapa de alimentación de acetona 43.
alimentación azeotrópica y para acetona.
Instale una columna de destilación; utilice el diagrama de flujo secundario
de la columna de destilación.
Defina la columna de destilación usando los mismos pasos 6.4 Destilación reactiva
explicados en el ejemplo 6.4. Sin embargo, en este ejemplo, se
deben considerar dos corrientes de entrada: la alimentación La destilación reactiva combina reacción química y destilación en un
azeotrópica ingresa a una de las etapas superiores y la acetona solo equipo. La destilación reactiva también se puede utilizar para
ingresa cerca del fondo de la columna. Como primera estimación, separar mezclas azeotrópicas o mezclas de punto de ebullición
utilice 22 etapas teóricas; para la mezcla azeotrópica seleccione la cerrado, en las que un nuevo componente añadido reacciona
etapa 7 como la etapa de alimentación y para la alimentación de químicamente con uno de los componentes y crea un componente
acetona use la etapa 18. La figura 6.24 muestra el diagrama de flujo del con
proceso.
diferentes propiedades físicas, que se puede separar fácilmente
Como especificación inicial de la columna, use el valor de para la de la mezcla. Sin embargo, la aplicación más interesante de la
Relación de 100 −1 tasa de destilado y 5 para el reflujo destilación reactiva es en reacciones de equilibrio limitado donde uno
kmol h. o más productos formados pueden eliminarse continuamente para
Ejecute la simulación con las especificaciones iniciales y verifique lograr una alta conversión. Las reacciones químicas suelen tener
los resultados. El destilado contendrá un azeótropo de acetona/ lugar en fase líquida o sobre una superficie sólida de catalizador.
ciclohexano y el resto
Nombre Acetona F D EN
En los procedimientos de cálculo de la destilación reactiva, se debe A la octava etapa teórica entra 4 mol% de agua.
agregar un término de reacción a las ecuaciones de balance de materia Ambas corrientes de entrada ingresan a la columna como líquidos con
y energía de las etapas reactivas, la ecuación Mi,j (ecuación 6.14) se punto de burbujeo a 1 bar. Las etapas reactivas son 810, donde la
transforma así en reacción química
Mi,j = nL j−1 xi,j−1 + nVj+1 yi,j+1 + nFjxFi,j −(nLj + CH3COOH + CH3CH2OH → CH3COOC2H5 + H2O (R6.1)
nSLj)xi,j − (nVj + nSVj)yi,j
Utilice el modelo de operación de la unidad RadFrac de Aspen Plus o y los perfiles de composición y temperatura en la columna.
AGUA
DIC
B1
D1
SP
ETH
D
ORG.
H2O
RD
ÁRBITRO
Y
EN
Figura 6.25 Diagrama de flujo de destilación reactiva del proceso de acetato de etilo
Machine Translated by Google
Definir la columna de destilación reactiva; en la pestaña Configuración Establezca la presión en 1 bar en la columna de la pestaña Presión
en Configuración, defina el tipo de cálculo, el número de etapas, el en Configuración.
tipo de condensador y rehervidor, las fases válidas y el método de En la pestaña trifásica en Configuración, identifique las etapas que
convergencia; si los vapores que salen de la columna se condensan se probarán para dos fases líquidas; como etapa inicial, seleccione
fuera de la columna y se separan en dos fases líquidas, seleccione la etapa 1 y como etapa final, seleccione la etapa 2; como
Ninguno para el condensador; como fases válidas seleccionar componentes clave para identificar la segunda fase líquida, seleccione
VaporLíquidoLíquido. acetato de etilo y agua.
Debido a que la mezcla azeotrópica se destila, use el método Vaya a la pestaña Reacción en la columna de destilación
azeotrópico para la convergencia (Figura 6.26). Especificación y defina la sección reactiva de la columna y la
En la pestaña Corrientes en Configuración, defina la posición de las retención de líquido como se muestra en la Figura 6.27.
corrientes de entrada (etapa 8 para ácido acético, etapa 10 para La reacción química que tiene lugar en las etapas reactivas se
etanol y etapa 1 para la corriente de reflujo). selecciona en el paso 3 (Figura 6.27); sin embargo, no
se definió la reacción química, haga clic en el campo debajo de mezcla de ácido acético y agua con una pequeña cantidad de
ID de reacción, seleccione Nuevo y elija un nombre para la etanol y acetato de etilo; los vapores que salen de la columna (D)
reacción. contienen principalmente la mezcla azeotrópica ternaria de
Vaya a la pestaña de reacción seleccionada en Reacciones, haga acetato de etilo, etanol y agua; el producto destilado final (ETH)
clic en Nuevo y defina la estequiometría de la reacción; como tipo contiene principalmente acetato de etilo con cantidades
de reacción seleccione Equilibrio como se muestra en la Figura considerables de etanol y agua (véanse los números resaltados
6.28. en la Tabla 6.6).
Especificar otros bloques de operación de la unidad incluidos en Para preparar los perfiles de temperatura y composición en la
el diagrama de flujo del proceso (enfriador por la temperatura de columna, desplácese a Perfil en la operación de la unidad de
salida de 25 ◦C y presión de 1 bar, decantador por la temperatura columna RD. La Figura 6.29 muestra el perfil de temperatura y la
de salida de 25 ◦C y presión de 1 bar, divisor por el flujo de −1 Figura 6.30 el perfil de composición en la columna.
mol del Flujo de ETH de 40 kmol h Ejecute la ). Ambos perfiles muestran claramente el efecto de la reacción
simulación de la columna y, en caso de convergencia, verifique química que tiene lugar en las etapas 810.
los resultados. El producto inferior (W) es un
De RD B1 RD DIC DIC SP SP
A RD RD DIC B1 DIC SP RD
Fase: Líquido Líquido Líquido vapor Líquido Líquido Líquido Líquido Líquido Líquido
H2O −1 kmol h 2 9 120 59.66 59.66 12,58 145,56 34.10 6.74 27.36
Flujo de topos −1 kmol h 50 60 120 242.46 242.46 30,00 160,00 202.46 40.00 162.46
−1
Flujo de masa kg h 2.918,55 2.511,66 2.161,83 15.914,42 15.914,42 1.249,84 3.453,22 14.623,03 2.889,06 11.734,04
Caudal volumétrico l min−1 51.55 56.15 38,36 111.599,00 293.31 23.35 60.27 270.14 53.37 216.77
Temperatura ◦C 115.68 77.85 80.00 70.61 25.00 97.41 25.00 25.00 25.00 25.00
100
95
90
85
Temperatura
(°C)
80