TecNM:

Tecnológico
Nacional de México

Instituto  

Tecnológico de
Aguascalientes

Departamento de
Ingeniería Química y
  Bioquímica
Notas sobre el curso
"Simulación de
Procesos  
DCIQ Javier Gerson
Zapiain Salinas 

Agosto de 2017 
Capítulo 3

Simulación comercial.

El uso de un simulador comercial en la solución de problemas de ingeniería de procesos

simpli…ca la solución; sin embargo el simulador no necesariamente podría tener alguna unidad
(operación unitaria o reactor) con las consideraciones particulares requeridas. Se debe de
considerar que un software de este tipo requerirá caracteristicas especí…cas de hardware.
Además de estas desventajas, está el proceso del aprendizaje del manejo de simulador, aún
R R R R
cuando la mayoría de ellos (Aspen Plus , Aspen HYSYS , ProSim , ChemCad , ect.)
presentan una serie de pasos similares en la elaboración de proceso, hay consideraciones
particulares en cada uno de ellos.
R
En el caso particular del simulador Aspen Plus , este presenta las mejores características
comparado con los otros simuladores comerciales; tales como, una base de datos amplia, un
conjunto muy grande de modelos termodinámicos, así como una biblioteca de operaciones
unitarias y reactores que permite elaborar la gran mayoria de los proceso, además de una
seríe de herramientas de análisis de proceso, entre otras La intención de esta Unidad es
mostrar de manera simple la mayoría de estas características y ejempli…car un par de casos
de estudio.

3.1. Manejo de un simulador comercial.

El nombre de Aspen, es un acrónimo de Advanced System for Process Engineering (Sis-
tema Avanzado para Ingeniería de Procesos). Este software fue desarrollado en 1970 en el
MIT.
Aspen Plus es un programa básicamente de simulación de procesos químicos, en el
cual además de simulaciones de diagramas de ‡ujo, se puede realizar: (1) estimación de
propiedades de compuestos (2) análisis de sensibilidad de variables de proceso , (3) obtener
especi…caciones de diseño de proceso, y (4) síntesis y análisis de procesos químicos, entre
otras tareas del diseño de procesos y equipos.

75
76 CAPíTULO 3 SIMULACIÓN COMERCIAL.

3.1.1. Introduccion al Aspen

El simulador Aspen Plus es un simulador tipo hoja de ‡ujo (‡owsheet), esto es un software
de computadora que se utiliza para modelar cuantitativamente una planta de procesamiento
químico, es decir un conjunto de reactores y operaciones unitarias. Este conjunto de procesos
se utiliza para transformar la materia prima en uno o más productos. Así, la simulación de
todo un proceso químico está representada simbólicamente por diferentes iconos donde cada
icono representa una unidad de operación, proceso químico, material de entrada / salida,
corriente, ‡ujo de energía de entrada / salida o señal eléctrica / neumática de entrada / salida.
En términos de la notación de la hoja de ‡ujo de Aspen Plus, habrá un icono de bloque y un
icono de ‡ujo. Los simuladores como Aspen Plus, nos permite predecir el comportamiento
de un proceso.
Los simuladores de hoja de ‡ujo, como Aspen Plus, nos permite predecir el compor-
tamiento de un proceso usando relaciones de ingeniería básicas.
Como ya se reviso en la Unidad I, podemos describir un proceso físico dado utilizando un
conjunto de ecuaciones algebráicas y/o diferenciales linealmente independientes tales que el
número de ecuaciones escritas será igual al número de variables (o cantidades desconocidas).
En general, la escritura de estas ecuaciones se deriva de:
ecuaciones de equilibrio de propiedades termodinámicas extensivas, tales como masa,
mol y energía;
relaciones termodinámicas para el medio reactivo y no reactivo, como las de fase y las
de equilibrio químico;
correlaciones de velocidad para el momento, el calor y la transferencia de masa;
estequiometría de reacción y datos cinéticos;
restricciones físicas impuestas al proceso.
En general, un proceso químico consiste en componentes químicos, o especies diferentes,
que están sujetos a tratamiento físico o químico, o ambos. El objetivo de aplicar dicho
tratamiento es básicamente agregar un valor o convertir el material bruto, barato, en valiosos
productos acabados. Las etapas de tratamiento físico pueden incluir el mezclado, la sepa-
ración (tales como absorción, destilación y extracción) y calentamiento/enfriamiento con o
sin un cambio de fase. Por otro lado, la etapa de tratamiento químico implica un conjunto
de reacciones paralelas, en serie o mixtas, lo que resulta en un cambio de identidad química
de cada una de las especies que reaccionan. Tales pasos de tratamiento se visualizan en el
simulador de diagrama de ‡ujo como componentes que son transportados de una unidad (o
bloque) a otro a través de ‡ujos de proceso.
Podemos traducir un proceso en un modelo de simulación de proceso Aspen Plus real-
izando los siguientes pasos necesarios:

1. Especi…car los componentes químicos en el proceso. Este procedimiento se realiza a)

utilizando las bases de datos de Aspen Plus, o b) de…niendo componentes especí…cos
utilizando la plataforma Aspen Propiertie.

2. Especi…car los modelos termodinámicos para representar las propiedades físicas de los
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 77

componentes y mezclas en el proceso. Estos modelos están incorporados en el simulador.

3. De…na el diagrama de ‡ujo del proceso:

De…nir las operaciones de la unidad en el proceso, una vez de…nidos se seleccionan

los modelos utilizando Aspen Plus Model Library para describir la operación de cada
unidad o síntesis química y colocarlos en el diagrama de ‡ujo del proceso.
De…nir los ‡ujos de proceso que ‡uyen dentro y fuera de las operaciones de la unidad.
Etiquetar cada modelo de operación de la unidad (es decir, bloque) como parte del
diagrama de ‡ujo del proceso y conectar los bloques a través de los ‡ujos de proceso.

4. Especi…car los ‡ujos de los componentes y las condiciones termodinámicas (temper-

atura, Presión y composición) de todas las corrientes de alimentación.

5. Especi…que las condiciones de funcionamiento de los modelos de funcionamiento de la

unidad (es decir, bloques).

3.1.2. Comenzando una simulación

Para describir de forma general los pasos a seguir para la elaboración de una simulación
utilizando el Aspen, se plantea el siguiente ejercicio.

Ejemplo 29 Proceso de obtención del cumeno a partir de benceno. Considere el procesos

que se muestra en la …gura 3.1, se alimenta una corriente de benceno y propileno a un reactor
donde se obtiene cumeno (propilbenceno), al mezcla que sale del reactor se enfría para ser
alimantada a una cámara de ‡ash, ya que hay una cantidad signi…cativa de reactivos que
debe ser retornada al rector, en una corriente de reciclo. La corriente líquida que sale del
‡ash es la del producto deseado.

Figura 3.1: Proceso de obtención del cumeno a partir del benceno

Las características de la corriente de alimentación se muestra en la tabla 3.1, mientras

que las condiciones para cada equipo se muestran en la tabla 3.2.
78 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Tabla 3.1: Características de la corriente de alimentación

Corriente FEED
Temperatura, F 220
Presión, psia 36
Flujo por componente, lbmol /hr
Benceno 40
Tolueno 40

Tabla 3.2: Especi…caciones de los equipos de procesos

Equipo Reactor Enfriador Flash
Temperatura, F Adiabático 130 Adiabático
Presión, psia 3 1 atm
Caida de presión, psi 0 0.1

El procedimiento para la elaboración de este proceso es:

PASO 1.- Se inicializa el simulador Aspen, y se selecciona NEW, esta acción despilega
una pantalla en donde se selecciona el sistema de unidades a emplear en el proceso; el sistema
se elige de acuerdo a los datos del problema, en este caso el sistema ingles de unidades, ver
…gura 3.3 y posteriormente se pulsa la tecla Create.

Tabla 3.3: Selección de plantillas de sistemas de unidades.

PASO 2.- Se deben de capturar todos los componentes que aparecerán en el proceso,
el orden de captura no es importante; sin embargo se debe de escribir un identi…cador
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 79

(component ID), al terminar deberá tener el aspecto que se muestra en la …gura 3.2, a
continuación se pulsa la tecla Next.

Figura 3.2: Formato de captura de componentes.

PASO 3.- Se selecciona el método de propiedades. La elección del método se realizará

como se discutió previamente, para este caso se puede utilizar la ecuación de Peng-Robinson.
Este procedimiento se realiza buscando en el menu deplegable Base Method la ecuación
requerida, ver 3.3. Al terminar pulsar la tecla Next.

Figura 3.3: Selección del método de propiedades.

PASO 4.- Una vez terminado los pasos anteriores se selecciona el modelo de fun-
cionamiento Simulation y se elabora el diagrama del proceso seleccionado primero las
operaciones unitarias requeridas, que en el simulador se les llama Block, a partir de la
80 CAPíTULO 3 SIMULACIÓN COMERCIAL.

biblioteca, y despúes concetar los equipos, cuidando que las conecciones se hagan de forma
adecuado, el resultado …nal se verá como se muestra en la …gura 3.4. Al terminar pulsar la
tecla Next.

Figura 3.4: Hoja de Flujo (Flowsheet) del proceso de obtención de cumeno.

PASO 5.- Como un paso opcional, ya que este paso no es necesario para la resolución
del simulador, en la ventana Simulation se busca la carpeta Setup>Speci…cations; en
esta carpeta se escribirá el nombre de la actividad y el nombre del realizador, como se ve
en la …gura 3.5.a. A continuación se abre la carpetal Setup>Speci…cations; y en ella se
selecciona la opción Fraction basis>Mole, ver …gura 3.5.b, esta acción permite que al
terminar la simulación puedan ser muestradas las fracciones molares de los componentes en
todas las corrientes. Al realizar lo anterior pulsar la tecla Next.

a) b)

Figura 3.5: Pantallas opcionales para las simulaciones. a) Título de la simulación b) infor-
mación ha mostrar de las corrientes de proceso.

PASO 6.- A continuación se deberán a continuación llenar las formas para cada una de
las corrientes de alimentación, en el caso que se está desarrollando solo se tiene una, FEED;
al capturar los datos de la tabla 3.1 tendrá el aspecto que se muestra en la …gura 3.6. En caso
de haber más de una alimentación este paso se realizará por cada alimentación. Al terminar
pulsar la tecla Next.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 81

Figura 3.6: Formato para la captura de los datos de la alimentación.

PASO 7.- El paso siguiente es llenar las formas para cada una de los bloques, que
aparecerán en orden alfabético, en este caso el proceso consta de tres bloques, que fueron
nombrados: COOL, REACTOR y SEP. Los datos de los equipos aparecen en la tabla
3.2. Para el caso del enfriador, bloque COOL en la …gura 3.7 se muestra la información
relaciona. Al terminar pulsar la tecla Next.

Figura 3.7: Datos requeridos para el enfriador, bloque COOL.

Para el caso del reactor, bloque REACTOR, el modelo seleccionado es un reactor es-
tequiométrico, el cual requiere además de las condiciones de operación, ver tabla 3.2, la
estequiometria de la reacción, es decir la reacción química balanceada. En la …gura 3.8a se
muestra el formato para las condiciones de operación, o especi…caciones.Para el caso de la
82 CAPíTULO 3 SIMULACIÓN COMERCIAL.

reacción química, se pulsa la carpeta Reaction, el cual abrirá el administrador de reac-

ciones, de…niendo una neva reacción, NEW, se desplegará el aditor estequimetrico, ambas
pantallas se muestran en la …gura 3.8b. En el editor se introducen los reactivos, los productos
y los coe…cientes estequimétricos de ambos. Así mismo se deberá proporcionar la fracción de
conversión y el reactivo de referencia al terminar pulse la tecla Next.

a) b)

Figura 3.8: a) Condiciones de operación para el reactor estequiométrico, bloque REACTOR.

b) Editor de estequimetría del reactor estequimétrico.

Finalmente, para el caso de la cámara de ‡ash, bloque SEP, las condiciones se muestran
en la …gura 3.9.

Figura 3.9: Datos requeridos para el separador ‡ash, bloque SEP.

PASO 8.- Al realizar la simulación se mostrará la pantalla del Control Panel, donde
se muestra los pasos que realizados por el simulador, ver …gura 3.10, en esta pantalla no
parecen los resultados númericos de las corrientes involucradas. La informacción mostrada
en el panel de control esta relacionada con la secuencia de solución, información particular
para cada bloque, información sobre los métodos de convergencia empleados y, en caso de
existir, mensajes de alertas o errores en solución del simulador.
PASO 9.- Para visualizar los resultados de las corrientes se debe pulsar Results Sumary
>Stream, en donde se mostraran los resultados de todas las propiedades de las corri-
entes del proceso, la cual incluye propiedades termodinámicas, es recomendable seleccionar
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 83

Figura 3.10: Resultados del proceso de simulación, Panel de Control.

Format>GEN_E, esta instrucción presenta solo los datos termodinámicos más relevantes,
los ‡ujos molares y las fracciones molares, ver …gura 3.11.

Figura 3.11: Resultados de la simulación para cada una de las corrientes.

Para que esta tabla aparezca en la hoja de ‡ujo se pulsa el botón Stream Table, con la
cual se genera la salida mostrada en la …gura 3.12.
PASO 10.- Finalmente la simulación deberá ser grabada, utilizando el icono respectivo
o utiizando el el menu File>Save o File>Save As. Todos los ejemplos se grabaran con el
nombre de ActividadXX, donde XX es el número correspondiente de actividad, se debetener
cuidado de veri…car que al momento de salvar el archivo este tenga la extensión apw. Es
recomendable crear una carpeta para guardad todos los archivos de Aspen y tener un acceso
84 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.12: Hoja de ‡ujo con los resultados de las corrientes.

rápido. Este ejemplo desarrollado se grabará como actividad00.apw, como se muestra en

la …gura 3.13.

Figura 3.13: Forma para grabar la simulación.

Los pasos anteriormente descrito se utilizarán para la mayoría de los ejemplos que se
describen a continuación y por lo tanto se hará mención de este procedimiento. En estos
casos, lo único que cambia, son las especi…caciones para cada bloque.

3.1.3. Termodinámica y Método de Propiedades.

Uno de los aspectos principales en la simulación es la selección adecuada de las ecuaciones
para el cálculo de las propiedades termodinámicas que serán utilizadas en la simulación. El
proceso de selección de estas ecuaciones ha generado una gran cantidad de artículos y libros,
sin embargo cada ingeniero tiene sus propios criterios. Una de las estrategias de selección
más utilizada ha sido desarrollada por Eric Carlson (Aspen Plus), la cual se muestra a
continuación en la …gura 3.14.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 85

Figura 3.14: Algoritmo de Carson para la determinación del método de propiedades

En términos generales, en el simulador Aspen, las especies químicas se clasi…can como:

Polares.- Sustancias que presentan una separación de carga dentro de la misma

molécula.

Electrolíticos.- Sustancias que pueden formar iones o que presentan disociación en

caso de haber agua. (Lavado de cenizas, neutralización, producción de ácidos y precipita-
ciones de sal)
No-electrolíticos.- Compuestos orgánicos con enlaces covalentes polares y/o cova-
lentes coordinados. (Aminas, Amidas, Cetonas, Aldehídos, Alcoholes, etc.)

No polares.- Compuestos orgánicos con enlaces covalentes.

86 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Reales.- Componentes puros. (Alcanos, alquenos, alquinos, cíclicos, aromáticos,

etc.)
Pseudocomponentes.- Mezclas no polares en donde algunos de los componentes
no pueden ser identi…cados. (Petróleo, gasolinas, etc.)
En esta sección se van a desarrollar ejemplos de la aplicación de Aspen para la elaboración
de diagramas de mezclas binarias y terciarias, en todos los casos se utilizará el algoritmo de
Carson para seleccionar el método de propiedades.

Ejemplo 30 Diagrama T x y y Diagrama de equilibrio para un sistema ide-

al. Elaborar un diagrama de ebullición y un diagrama de equilibrio para el sistema ben-
ceno/tolueno, la presión de operación es 1 atm. Al terminar salve su simulación con el
nombre de actividad01.apw

Una vez que han sido capturados los componentes benceno y tolueno (Paso 1 y Paso
2), se selecciona el método térmodinámico (Paso 3). Se selecciona Binary Analysis, donde
se selecciona en Analysis Type la opción Txy y se introduce el valor de la presión y se
seleccionan la unidad de atm. En la …gura 3.15 se muestra el aspecto del analisis. Una vez
llenado la forma se pulsa Run Analysis.

Figura 3.15: Plantilla de Binary Analysis. Mezcla Benceno/Tolueno

Como resultado se genera automáticamente el diagrama de ebullición, ver …gura 3.16,a;

utilizando las opciones de Plot, se puede obtener el diagrama de equilibrio, que se muestra
en la …gura 3.16,b. Finalmente se guarda la actividad realizada.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 87

a) b)

Figura 3.16: a) Diagrama de ebullición, mezcla benceno/tolueno; b) Diagrama de equilibrio,

mezcla benceno/ tolueno

Ejemplo 31 Diagrama T x y, diagrama de equilibrio y diagrama de coe…ciente

de actividad para un sistema no-ideal. Elaborar un diagrama de ebullición, un diagrama
de equilibrio y un diagrama de coe…ciente de actividad para el sistema etanol/agua, la presión
será de 1 atm. Al terminar salve su simulación con el nombre de actividad02.apw.

Se seleccionan los compontes etanol y agua (Paso 1 y Paso 2) y el modelo termodinámico

(Paso 3) y se completa el formato de Binary Analysis,ver …gura 3.17, de forma similar al
ejercicio anterior. Despúes se realiza el análisis.

Figura 3.17: Plantilla de Binary Analysis. Mezcla etanol/agua.

88 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Los resultados de esta análisis se muestran en la …gura 3.18.

a) b)

Figura 3.18: a) Diagrama de ebullición, mezcla etanol/agua; b) Diagrama de equilibrio,

mezcla etanol/agua.

Para visualizar los resultados númericos, estos se localizan en la barra de navegación

Analysis>BINARY-1>Results y se muestran en la …gura 3.19a, mientras que para de-
splegar el gra…co del coe…ciente de actividad se selecciona en Home>plot>Activity Coe¤.
y se mostrará el gra…co de la …gura 3.19b. Para terminar se guarda la actividad.

a) b)

Figura 3.19: a) Resultados númericos del análisis binario, sistema etanol/agua. b) Diagrama
de coe…ciente de actividad, sistema etanol/agua.

Ejemplo 32 Digrama triangular. Elaborar un diagrama triangular para el sistema agua/n-

octanol/n-hexano a 1 bar de presión. Al …naliza guardar como actividad03.apw.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 89

Los componentes a considerar, agua, n-octanol y n-hexano, son capturados de la forma

en que se ha realizado anteriormente (Paso 1 y Paso 2) y se selecciona el método térmod-
inámico (Paso 3). En el menú de Home se selecciona Ternary Diag y la forma que se
despliega se deberá llenar de forma adecuada, tal como se muestra en la …gura 3.20. Una vez
que se ha completado el llenado, pulsa Run Analysis.

Figura 3.20: Plantilla de Ternary Diag. Mezcla agua/n-butanol/n-hexano.

De forma automática aparecera el diagrama triangular, con tantas líneas de reparto como
se haya seleccionado en la forma anterior (Number of tie lines), ver …gura 3.21.

Figura 3.21: Diagrama Triangular para el sistema agua/n-butanol/n-hexano a 1 bar de pre-

sión.

Adicionalmente se puede obtener información sobre los azeotropos encontrados para el

sistema al seleccionar Analysis>Terdi-1>Result, que depliega los resultados del sistema.
Si se selecciona la opción Azeotrope se muestran las concentraciones de estos puntos a las
condiciones de presión especi…cadas, ver …gura 3.22. Al concluir guardar la actividad.
90 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.22: Localización de los azeotrópos para el sistema agua/n-butanol/n-hexano a 1

bar de pre- sión.

Aún cuando la base de datos de Aspen es extensa, existen componentes no incluidos en

ella; sin embargo es posible incluir estos componentes conociendo algunas de sus propiedades
físicas, utilizando Aspen Properties. Para conlcuir esta sección se describirá la forma
de incluir especies químicas, se debe de considerar que esta modi…cación en la base de
datos será valida solo en el archivo de trabajo, cada nueva simulación que requiera nuevos
componentes deberán ser dados de alta con este procedimeintos y preferentemente al inicio
de la simulación.

Ejemplo 33 Predicción de propiedades. La vainillina, C8 H8 O3 , es un agente saborizante

en comidas, bebidas y productos farmacéuticos, que esta de…nida en la base de datos de Aspen.
Obtenga sus propiedades utilizando Aspen Properties, si se sabe que esta sustancia tiene
un peso molecular de 152.15 gr/mol, una temperatura normal de ebullición de 284.85 C y
una densidad de 1.056 g/cm3 . La estructura química de la vainillina se muestra en la …gura
3.23. Al terminar salve el archivo con el nombre actividad04.apw.

Figura 3.23: Estructura de la vainillina.

Paso A.- Iniciar una simulación en blanco, y de…nir el sistema de unidades más adecuado
a los datos de las propiedades conocidas (Paso 1). En el formato de componentes, escribir
solo el Component ID de la sustancia a predecir y ver…car que el type sea Conventional.
En este paso solo deberá capturarse los componentes desconocidos, tal como se muestra en
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 91

la …gura 3.24. En la práctica se puede hacer con varios componentes a la vez, tomando en
cuenta el cuidado al capturar las propiedades de cada especie diferente.

Figura 3.24: Formato de componentes.

Paso B.- Se recomienda utilizar el método de propiedades de NRTL. En el menú de

navegación selccionar Methods>Parameters>Binary Interaction>NRTL-1 y veri…car
que este seleccionada Estimate parameters by UNIFAC; ver la …gura 3.25.

Figura 3.25: Formato de Binary Interaction.

Paso C.- Se deberá proporcionar la estructura del componente, para lo cual puede hac-
erse de varias formas a) utilizando un editor de moleculas como ChemO¢ ce, ChemDraw,
etc., los archivos *.mol son compatibles con Aspen Property, b) utilizando una codi…-
cación simple de la estructura y c) utilizando el Molecule Editor, que es la opción que se
mostrará. Para acceder a este editor se pulsa el icono Draw Structutre de menu Home. El
funcionamiento de esta aplicación es muy sencilla, ya que existen estructuras prediseñadas
92 CAPíTULO 3 SIMULACIÓN COMERCIAL.

(Fragments), las cuales facilitan el dibujar la estructura, iconos para enlaces sencillos, dobles
o triples, e iconos para selccionar el tipo de átomo a dibujar. En la …gura 3.26a se muestra
el resultado de la estructura. Finalmente esta guarda con un nombre adecuado, para este
caso vainilla.mol y se cierra el editor. El formato de Molecular Structure se mostrará
como en la …gura 3.26b. Finalmente se deberá pulsar el botón Calculate Bonds, para que
se obtenga la estructura desarrollada de la molecula.

a) B

Figura 3.26: a) Formato Molecular Editor y b) Formato Structure.

Paso D.- Una vez de…nida la estrucutura se deberán capturar las propiedades conocidas
de la sustancia, para esto es necesario acceder a Methodos>Parameters>Pure Compo-
nents y pulsar New. En la ventana que aperecerá, New Pure Component Parameters,
dar el nombre del nuevo componente, en este caso VAINILLA, el uso de mayúsculas o minús-
culas es indistinto y no afecta la predicción, ver la …gura 3.27a. Al pulsar el boton Ok, se
desplagará una forma de captura. En la columna Parameters, se seleccionará la propiedad
a capturar. Mientras que en la columna Component se seleccionará el nombre del compo-
nente. Una vez seleccionado la propiedad a capturar, el valor númerico de esta propiedad
estará debajo del nombre del componente y en la casilla correspondiente capturar el valor
de la propiedad. En la columna Units se pueden seleccionar las unidades adecuadas para
el valor conocido, al …nializar la captura de datos de este ejercicio la forma será igual a la
mostrada en la …gura 3.27b.

Paso E.- Para terminar con este proceso se pulsa el icono RUN, para que el sistema
genere las propiedades. Los valores estimados se muestran en Methodos>Parameters>Pure
Components>PCES-1, ver la …gura 3.28. En la mayoría de las veces aparecerán mensajes
de precaución, pero no son de importancia. Finalmente guardar la simulación.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 93

a) b)

Figura 3.27: a) Forma Pure Components y b) Forma Input

Figura 3.28: Propiedades estimadas de la vainillina.

3.1.4. Mezcladores y Divisores de ‡ujo.

Los elementos más simples en Aspen son los mezcladores de corriente (Mixer) y los
divisores de ‡ujo (FSplit). Los iconos de estas operaciones se muestran en la …gura 3.29.
Los mezcladores combinaran las corrientes de entrada, realizando los balances de materia
y energía correspondientes. Se pueden utilizar tantas alimentaciones como sea requerido; para
el caso del dividor de ‡ujo sólo fraccionará las salidas, pero no cambiarán las composiciones
de las corrientes de salida, en este caso se deberá especi…car el porcentaje de salida para
m 1 corrientes de salida, donde m es el número de salidas del divisor.
94 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.29: Mezclador y divisor de ‡ujo. Paleta de Modelos.

Ejemplo 34 Mezclador de corrientes. Considere el sistema que se muestra en la …gura

3.30, el cual consiste de un mezclador (Mixer) en donde se alimentan tres corrientes. Los
componentes en este proceso son propano (C3), n-butano(nC4), n-pentano (nC5) y n-hexano
(nC6). Las especi…caciones de las corrientes se muestran en las tablas 3.4 y 3.5. Al terminar
salve el archivo con el nombre actividad05.apw.

Figura 3.30: Mezclador con tres alimentaciones

Tabla 3.4: Condiciones de las corrientes de alimentación del mezclador

Corriente 1 2 3
Temperatura ( F) 200 250
Presión (psi) 300 400
Fracción vapor 200 0.5

En los mezcladores solo se especi…ca la presión, en caso de no hacerlo, como en este

ejercicio, la presión del equipo y de la corriente de salida será la presión más baja de todas
las alimentaciones. Antes de iniciar la simulación se selecciona, utilizando el algoritmo de
Carlson, el método termodinámico adecuado. Para este ejemplo se repiten los Pasos del 1
al 6 discutidos con anterioridad, se omite el Paso 7.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 95

Tabla 3.5: Flujos molares por componente para las corrientes de alimentación para el mez-
clador.
Corriente Flujo (lbmol /hr)
C3 nC4 nC5 nC6
1 20 30 30 20
2 30 30 20 20
3 50 0 30 20

Finalmente se concluyen con los ultimos Pasos, del 8 al 10. Los resultados de la simu-
lación se muestra en la …gura 3.31.

Figura 3.31: Resultados del mezclador de corrientes.

3.1.5. Manipuladores de Presión.

Se disponen en Aspen de un conjunto de equipos que permiten varíar la presión en las
corrientes de un proceso, ya sea para incrementarla, como las bombas (Pump) y compresores
(Compr), reductores de presión como las válvulas (Valve) y los conjuntos complejo de
tuberías que incluyen distintos accesorios (Pipe), entre otros. Los iconos de estos equipos,
o cambiadores de presión (Pressure Changers) se muestran en la …gura 3.32.
Todos estos modelos permiten interconectar equipos de reacción o separación; normal-
mente se utilizan en combinación con los intercambiadores de calor para adecuar las condi-
ciones de la presión y la temperatura de las corrientes. En la simulación rigurosa de una
planta de procesos deberían incluirse estos equipos, aunque es posible realizar simulaciones
sencillas sin incluir este tipo de equipos.
96 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.32: Cambiadores de presión. Paleta de Modelos.

Ejemplo 35 Sistema de bombeo de agua.- Se desea simular el sistema de bombeo de

agua, mostrado en la …gura 3.33, en este hay dos tanques, que pueden ser representados con
dos mezclador, cuatro válvulas, una bomba y dos secciones iguales de tuberías. Se bombean
10,000 kg/hr de agua, que se encuentra a 25 C y 2 atm de presión. Los sistemas de tuberías
PIPE-1 y PIPE-2 son idénticos, son secciones de acero de cédula 40S, 1 km de longuitud,
de diámetro nominal de 2.5 in y que operan de forma adiabática, los accesorios de estos
sistemas, que están unidos por tornillos, se muestan en la tabla 3.6. Las condiciones de las
válvulas se muestran en las tablas 3.7 , para el caso de la bomba sólo se requiere que la presión
de descarga sea de 3 bares. Al término de la simulación grabar como actividad06.apw.

Figura 3.33: Sistema de bombeo de agua.

Tabla 3.6: Accesorios del sistemas de tubería del sistema de bombeo.

Válvulas de compuerta 2
Válvulas de mariposa 1
Codos largos de 90 6
T rectas 10

Cuando uno de los componentes es agua, y este no interviene como un reactivo, si no que
se utiliza en el sistema de enfriamiento o calentamiento, es más adecuado utilizar tablas de
vapor como método de propiedades, la cual aparece como STEAM-TA, como se muestra en la
…gura 3.34. En el simulador Aspen es posible seleccionar más de un método termodinámico.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 97

Tabla 3.7: Parámetros de las válvulas del sistema de bombeo.

Válvulas 1 2 3 4
Condiciones de operación
Tipo de cálculo A partir de especi…caciones de la válvula
% de apertura 50
Parámetros de la válvula
Tipo de válvula Mariposa Bola Globo Globo
Manufactura Neles - Jamesbury
Serie/estilo Ansi-300 Metal-Seated V810- % igual V810-Lineal
Tamaño 3 in 3 in 3 in 3 in

Se realizan los Pasos del 1 al 6 discutidos con anterioridad, para el caso del Paso 7 se
deberá capturar la información para cada uno de los equipos, con excepción de los tanque
(Mixer); para el caso de los sistemas de tubería (Pipe-1 y Pipe-2) que son iguales, las
formas de captura quedarán como se muestra en la …gura 3.35.
Para la bomba (Pump) la forma correspondiente tendrá el aspecto que se muestra en la
…gura 3.36.
Para cada una de las válvulas las que las condiciones de operación son las mismas, sin
embargo en el sistema hay tres tipos de válvulas diferentes, una tipo mariposa (Valve-1),
una de compuerta (Valve-2) y dos de globo (Valve-3 y Valve-4); en las …guras 3.37 y
3.38 se muestran las formas para cada una de ellas.
Finalmente se concluyen con los ultimos Pasos, del 8 al 10. En este sistema, los valores
de interés son los de la presión entre cada equipo, para mostrar los valores de la presión
en la hoja de ‡ujo se selecciona la casilla de Pressure del menú Modify, los resultados se
muestran en la …gura 3.31.

Figura 3.34: Selección para usar Tablas de Vapor como método de propiedades.
98 CAPíTULO 3 SIMULACIÓN COMERCIAL.

a) b)

Figura 3.35: a) Parametros de la tubería, b) Accesorios de la tubería.

Figura 3.36: Parámetros de la bomba.

a) b)

Figura 3.37: Parámetros de las válvulas a) Mariposa y b) Bola

3.1 MANEJO DE UN SIMULADOR COMERCIAL. 99

a) b)

Figura 3.38: Parámetros de las válvulas a) Globo (proporcional) y b) Globo (lineal)

Figura 3.39: Resultados de la presión para la simulación de la planta de bombeo.

3.1.6. Intercambiadores de Calor.

En este simulador se disponen de varios intercambiadores de calor, siendo los más comunes
el Heater y el HeaterX, lo cuales no requieren ninguna información sobre los parámetros
del equipo. El primero sólo permite una corriente, que puede ser modi…cada en su tem-
peratura con o sin cambio de fase; mientras que el segundo permite utilizar dos corrientes,
en este módulo es posible obtener el área de transferencia de calor, el coe…ciente global de
transferencia y las curvas de transferencia de calor. Los iconos de estos modelos se muestran
en la …gura 3.40.

Figura 3.40: Intercambiadores de calor. Paleta de Modelos.

A continuación se presentan tres ejemplos que permiten conocer el uso de estos equipos,
en el primero se desarrolla un calentador para una corriente de freón, en el segundo un
100 CAPíTULO 3 SIMULACIÓN COMERCIAL.

sistema más elaborado en donde una corriente de freón enfría a una corriente de etilenglicol
y …nalmente un acomparación entre los modulos de Heater y HeaterX

Ejemplo 36 Calentador de corriente. Considere una corriente de Freon-12, dicloro di-

‡uoro metano, de 12,000 kg/hr que se encuentra a 250 K y 8 atmósferas de presión, que
se alimenta a un intercambiador de calor. Se espera aumentar en 50 la temperatura de la
corriente. En la …gura 3.41 se muestra el proceso de calentamiento. Considere una caída
de presión de 0.1 bar. Estime la cantidad de calor que se requerida. Al terminar grabe la
simulación como actividad07.apw.

Figura 3.41: Calentador de una sola corriente, HEATER.

Para resolver esta simulación se deben considerar lo siguiente, que el freón-12 está en
la bases de datos como R-12 y para de…nir una caida de presión se debe de expresar con
un signo negativo, para este caso la caída se expresa como -0.1 bar. Se realizan los Pasos
del 1 al 6 discutidos con anterioridad, considerando en el Paso 2 que el nombre con que se
encuentra el Freón-12; para el caso del Paso 7 la forma para el intercambiador aparece en
la …gura 3.42, considerando la manera en que se captura la caída de presión.

Figura 3.42: Formato de Especi…caciones para el calentador HEATER.

Finalmente se concluyen con los ultimos Pasos, del 8 al 10. La simulación de este proble-
ma debe de proporcionar el calor removido, por lo cual se despliega los resultados del bloque,
Blocks>HEATER>Results, ver …gura 3.43. Por lo que el calor removido es igual 1800.49
kW.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 101

Figura 3.43: Resultados del calentador de freón.

Ejemplo 37 Intercambiador de calor.- Considere dos corrientes que entran a un inter-

cambiador; la primera, la caliente, es un ‡ujo de 5,000 kg/hr de etilenglicol a 370 K y 3
atmósferas de presión; la segunda, la fría, es una corriente de Freon-12 a 250 K y 8 atmós-
feras de presión. El diagrama del proceso se muesta en la …gura 3.44. se desea saber el valor
de la temperatura de la corriente de salida del etilenglicol (EG-OUT) si a la corriente fria
(FR12-IN) le es transferida 1800.049 kW de calor. Para este caso se utilizará un mode-
lo HEATX, con las parametros que se muestran en la tabla 3.8. Al términar grabar como
actividad08.apw.

Figura 3.44: Intercambiador de dos corrientes.

Tabla 3.8: Parámetros del intercambiador de calor.

Especi…cación
Modelo Shortcut
Dirección del ‡ujo Countercurrent
Modo de cálculo Design
Especi…cación Exchanger duty

La corriente fría es la misma que la del ejemplo anterior, por lo que el componente del
Freón-12 como ya se menciono anteriormente stá en la bases de datos como R-12. Se realizan
102 CAPíTULO 3 SIMULACIÓN COMERCIAL.

los Pasos del 1 al 6 discutidos con anterioridad, para el caso del Paso 7 la forma para este
tipo de intercambiador aparece en la …gura 3.45.

Figura 3.45: Formato de Especi…caciones del intercambiador HeatX.

A continuación se concluyen los Pasos, del 8 al 10. Los resultados de esta simulación
se localizan en Blocks>INTERCAM>Thermal Results, ver …gura 3.46a. Este bloque
también realizar los cálculos de área y coe…ciente, ver …gura 3.46b.

a) b)

Figura 3.46: Resultados del intercambiador a) temperatura de todas corrientes b) parámetros

de diseño.

Ejemplo 38 Comparación entre dos Heat y un HeatX. Se puede utilizar dos módulos
Heat para replicar el comportamineto de un HeatX, para este caso se utilizará una corriente
de agua para enfriar una corriente de hidrocarburos. El diagrama del proceso se muestra en
la …gura 3.47. La corriente de hidrocarburos se encuentra a 200 C y 4 bar de presión, con
un ‡ujo másico de 10000 kg/hr; la composición en porcentaje en peso de esta corrientes es
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 103

de 50 % de benceno, 20 % estireneo, 20 % etilbenceno y 10 % de agua. Mientras que la cor-

riente de agua de enfriamiento se encuentra a 20 C y una presión de 10 bar, en este caso
el ‡ujo másico es de 60000 kg/hr En este proceso no hay caidas de presión y la corriente
de hidrocarburo se matiene en fase vapor. Considere que el sistema de hidrocarburos puede
generar tres fases (dos líquidas y una de vapor). Se recomienda por las características del sis-
tema utilizar los métodos termodinámicos de NRTL-RK, para la corriente de hidrocarburos,
junto con las tablas de vapor, para el agua. Al términar la simulación se debe guarda como
actividad09.apw.

Figura 3.47: Comparación entre los módulos Heater vs HeatX.

Como se observa en la …gura 3.47 se utilizará nuevamente el módulo DUPL, una para
cada corriente de alimentación, se deberá considerar el uso de dos modelos térmodinamicos,
uno para la corriente de hidrocarburos y otro para la corriente de agua (se recomienda
ampliamente las tablas de vapor). A partir de los anterior se procede con los Pasos del 1
al 6 de la metodología propuesta, teniendo en cuenta que en el Paso 4 se deberá usar un
nuevo tipo de línea, que es la tipo Heat para unir los dos módulos HEATER. Este tipo de
línea indica una transferencia de calor, por lo que tiene una dirección, para esta simulación
el calor ‡uirá del Heater2 al Heater1; esta situación de…ne además que módulo requiere
más variables al momento de de…nir los equipos.
En del Paso 7 se deberán llenar los formatos para cada uno de los intercambiadores. Para
el caso de los HEATER, sólo es necesario de…nir tres variables entre los equipos; en este
caso la caída de presión (en ambos intercambiadores) y la condición de fase para el líquido
caliente, ver …guras 3.48a y 3.49a, lo anterior debido a que el calor ‡uye de la corriente
caliente a la fría.
104 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Para el intercambiador Heater donde ‡uye el líquido frio (el que transporta agua) deberá
especi…carse el método termodinámico para seleccionar las tablas de vapor, para lo cual se
selecciona Heater1>Block Options, ver …gura 3.48b.

a) b)

Figura 3.48: Formatos para el intercambiador Heater1. a) Especi…caiones b) Seleccion de

método termodinámico.

Para el caso del intercambiador de la corriente caliente (el que transporta los hidrocar-
buros), es necesario considerar las tres fases, por lo que se deberá seleccionar esa opcion
en la forma de Speci…cation, ver …gura 3.49a. Para esta corriente se selecciona el método
termodinamico NRTL-RK, seleccionando Heater2>Block Options, ver …gura 3.49b.

a) b)

Figura 3.49: Formatos para el intercambiador Heater2, a) Especi…caciones b) Seleccion de

método termodinámico.

Para el intercambiador HeatX, se deben de…nir todas las variables antes mostradas para
el sistema de dos intercambiadores simples, para esto se comienza de…nido las condiciones
generales de la simulación, de hecho este modelo permite incluir aspectos más complejo que
los discutido en este ejemplo, como la geometría y tipo del intercambiador, distribución de
tubos, la presencia de ba‡es, el tipo de material, etc. Para este ejercicio se utilizará el modo
más simple, que el Shorcut, y el modo de cálculo será: Desing. Las demás especi…caciones
se muestran en la …gura 3.50.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 105

Figura 3.50: Formato Especi…caciones para el para el intercambiador SHeatX.

Como en el caso anterior es necesario de…nir las fases presentes para cada una de las
corrientes involucradas en el intercambiador, para este equipos se hace según se muestra en
la …gura 3.50a, para acceder a esta forma se pulsa SHEATX>Options. Para seleccionar
los métodos termodinámicos para cada corriente, se pulsa SHEATX>Block Options y se
seleccionan de forma adecuada para cada corriente como se muestra en la …gura 3.50b.

a) b)

Figura 3.51: Formatos para el intercambiador SHeatX, a) especi…caciones de fase b) Seleccion

de método termodinámico..

A continuación se concluyen los Pasos, del 8 al 10. Los resultados de la simulación del sis-
tema de intercambiadores se localizan en HEATER1>Results y HEATER2>Results1,
los cuales se muestran en las …guras 3.52a.y 3.52b respectivamente. Mientras que para el
intercabiador HEATX, los resultados se localizan en Blocks>SHEATX>Thermal Re-
sults, ver …gura 3.53a. Como se comento anteriormente en este bloque también se realizan
los cálculos de área y coe…ciente, los cuales se muestran en la …gura 3.53b, pulsando la forma
Exchanger Details.
106 CAPíTULO 3 SIMULACIÓN COMERCIAL.

a) b)

Figura 3.52: Resultados para los intercambiadores de una corriente a) Heater1 y b) Heater2.

a) b)

Figura 3.53: Resultados para el intercambiador de dos corrientes a) temperaturas de salida

y b) detalles de diseño.

3.1.7. Reactores.
Uno de los modelos más importantes de una planta química son los reactores, en el
simulador Aspen se desarrollaron un conjunto muy amplio de reactores, los cuales se muestran
en la …gura 3.54. De acuerdo a la información requerida para con…gurar cada modelo, se
pueden clasi…car como:

1. Reactores de cinética simple.- En este caso se agrupan los modelos RYield, RGibbs,
REquil y Rstoic.

2. Reactores de cinética compleja.- Se encuentran los modelos RCSTR, RPlug y Rbatch.

A continuación se explica de forma amplia cada uno de estos tipos.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 107

Figura 3.54: Reactores. Paleta de Modelos

Reactores con cinética simple.

En este tipo de reactores no es necesario conocer la cinética del sistema reactivo, solo en
algunos de ellos la estequimetría de la reacción. En el RYield se basa solo en el redimiento
de los componentes. En el RGibbs se estima el equilibrio de fases a partir del criterio de la
minimización de la energía libre de Gibbs. En el reactor REquil se calcula el equilibrio de
fases pero en este caso hay que de…nir la reacción química (estequimetría). Finalmente en el
Rstoic requiere la estequiometría del sistema reactivo y además la conversión o avance de
la reacción. En el siguiente ejemplo se muestran las caracteristicas de estos modelos.

Ejemplo 39 Para mostrar las características de cada uno de los reactores de cinética simple
se puede considerar el caso de la reacción para la producción de hidrógeno a partir de metanol,
para la cual considere la siguiente reacción química CH4 O ! 2H2 +CO. Considere un ‡ujo
de metanol puro de 100 kmol/hr, las condiciones en todo el proceso de temperatura y presión
son de 150 C y 1 atm, respectivamente. En la …gura 3.59 se muestra el sistema de los cuatro
reactores, un RYield (PRODUCC), un RGibbs (GIBBS), un REquil (EQUILIB) y un
RStoic (ESTEQU); en esta con…guracion se utiliza un bloque DUPL. En la tabla 3.9 se
muestran la información requerida en los reactor. Al términar la simulación se debe guarda
como actividad10.apw.

Figura 3.55: Sistemas de reactores con cinética simple.

108 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Tabla 3.9: Parámetros requerido para los reactores de cinética simple

Reactor de redimiento.
Componente Producción molar
Metanol 0.61388
Hidrógeno 7.45262
Dioxido de Carbono 3.72631
Reactor estequiométrico.
Conversión con respecto al metanol 0.862

En este caso el uso del modulo DUPL es solo para no hacer cuatro corrientes de al-
imentación. En la solución de este sistema se siguen los Pasos del 1 al 6 discutidos con
anterioridad, para el caso del Paso 7 en las …guras se muestra los formatos para los reac-
tores.
En el caso de reactor de producción, RYield, la producción, en términos molares y en
términos másicos, se de…ne como:

moles del componente i en la salida

producción molar del componente i = (3.1)
moles totales en la alimentación

masa del componente i en la salida

producción másica del componente i = (3.2)
masa totales en la alimentación
Para este módulo, se llene al formato con la información proporcionada en los datos
previos, que en este caso la producción es molar, como se muestra en la …gura 3.56.

Figura 3.56: Párametros para el reactor RYiled.

Para el caso del reactor RGibbs, solo es necesario especi…car la presión y la temperatura
del reactor, ver …gura 3.57.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 109

Figura 3.57: Especi…caciones para el reactor RGibbs.

Mientras que para el reactor REquil, o reactor en equilibrio, se deben llenar dos formas,
la primera para especi…car las condiciones de temperatura y presión y las fases válidas, como
se muestra en la …gura 3.58a; para este modulo es necesario de…nir la estequimetría de la
reacción, al pulsar el botón Next se desplegara el administador de reacciones y al seleccionart
New se activa Edit Stoichiommetry, que es la forma donde se escribrirá la reacción
química. En general todas las reacciones se de…nen en términos de reactivos y productos,
primero se escriben los nombre de las especies y despúes el coe…ciente estequimétrico, en el
caso de los reactivos el valor deberá escribirse con un signo negativo y para los productos
tendrá que ser positivo, para este caso la forma del editor de estequimetría quedará como se
muestra en la …gura 3.58b.

a) b)

Figura 3.58: Especi…caciones para el reactor REquil, a) Condiciones de opreraciones y b)

Editor estequiométrico.

Para el rector estequimétrico, RStoic, el cual es similar al reactor anterior, se deben llenar
dos formatos, el de las condiciones de operación, …gura 3.59a y el del editor estequimétrico,
que en este caso deberá especi…carse ya sea la conversión o el avance de la reacción, en este
caso con los datos conocidos esta forma deberá tener el aspecto que se muestra en la …gura
3.59b.
110 CAPíTULO 3 SIMULACIÓN COMERCIAL.

a) b)

Figura 3.59: Especi…caciones para el reactor RStoic, a) Condiciones de opreraciones y b)

Editor estequiométrico..

A continuación se concluyen los Pasos, del 8 al 10. Los resultados de la simulación del
sistemas de los cuatro reactores se presentan en la …gura 3.60.

Figura 3.60: Resultados de los reactores de cinética simple.

Reactores con cinética compleja.

Para este tipo de reactores es necesario conocer en su totalidad la cinética química de la

reacción, en el simulador es posible realizar reacciones catalíticas y no-catalíticas, en forma
general los primeros utilizan una cinética tipo Langmuir, Hinshelwood, Hougen y Watson
(LHHW) y los segundo tipo ley de potencia (POWERLAW).
Las cinéticas menos complejas, que son las tipo ley de la potencias, se discutirán inicial-
mente; los del modelo LHHW o catalíticos se ejempli…carán posteriormente.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 111

Reactores con formas cinéticas tipo Ley de la Potencia. En estos modelos es nece-
sario conoce la ley de velocidad de una reacción. Considere, para ejempli…car la forma en
que se interpretan los datos cinéticos en el simulador, la reacción de descomposición de la
acetona, de acuerdo a la reacción:

CH3 COCH3 ! CH2 CO + CH4

Se sabe que esta reacción es de primer orden con respecto a la acetona, por lo que su ley
de velocidad o Ley de la potencia estará dada por:

r = k [CH3 COCH3 ]
en donde k es la constante de velocidad especí…ca, de acuerdo a la Ley de Arrhenius está
dada por:
E 1 1
k = k0 e R T T0
(3.3)
en esta expresión, k0 es el factor preexponencial; E, la energía de activación y T0 es la
temperatura de referencia.
Dado que las expresiones cinéticas son bastante complejas, en el simulador las reacciones
se capturan por separado. Los reactores de este tipo, sólo requieren las especi…caciones de
las condiciones de operación y las de diseño, como los son el volumen del reactor o el tiempo
de residencia (para los reactores tipo tanque), longitud y diámetro de los tubos, número de
tubos, etc. (para los tubulares).
Otro aspecto a cuidar en este formato es la base de la concentración, es decir en términos
de que está expresada la ley de velocidad; normalmente en términos de concentracion molar
(mol/lt) para las fases líquidas y presión parcial para los fases vapor o gas. En el simulador
se pueden utilizar las siguientes opciones, que se muestran en la tabla.3.10.

Tabla 3.10: Formatos para expresar la base de la concentración en las reacciones con cinética
compleja
[Ci] basis Término de concentración
Molarity Concentración molar de componente (kmol/m3 )
Molality Molalidad de componente (mol/kg de H2 O)
Mole fraction Fracción molar de componente
Mass fraction Fracción masa de componente
Partial pressure Presión parcial de componente (N/m2 )
Mass concentration Concentración masa de componente (kg/m3 )

El siguiente ejemplo permite conocer la forma en que esta información cinética se adapta
al simulador.
112 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Ejemplo 40 Una de las reacciones industriales más estudiada es la de síntesis del acetato
de etilo a partir de etanol y ácido acetico. La cual se lleva a cabo mediante la reacción
reversible:
CH3 CH2 OH + CH3 COOH CH3 COOCH2 CH3 + H2 O (3.4)
Estas reacciones ocurren a 70 C y 1 atmósfera de presión, y bajo estas condiciones se conocen
las siguientes leyes de velocidad para la reacciones directa y la indirecta, en fase líquida:
5;95 107
rD = 1;9 108 e RT [CH3 CH2 OH] [CH3 COOH] (3.5)
5;95 107
rI = 5 107 e RT [CH3 COOCH2 CH3 ] [H2 O] (3.6)
Las unidades de la energía de activiación son J/kmol. Se utilizarán dos reactores cinéticos
para simular la reacción. En la …gura 3.61 se muestra el diagrama de los reactores, el primero
es un CSTR y el segundo es un reactor FTR. Las parámetros de diseño de cada reactor se
muestra en la tabla 3.11. Considerese que la corriente de la alimentación esta compuesta por
8.892 kmol/hr de agua, 186.59 kmol/hr de etanol y de 192.6 kmol/hr de ácido acético. En este
sistema reactivo se recomienda utilizar el método termodinámico NRT-RK. En este ejercicio
se debe permitir la presencia de las fases vapor y líquido, aún cuando la reacción ocurre en
fase líquida. Considere que en ambas leyes de velocidad está en términos de conctraciones
molares. Al términar grabe la simulación como actividad11.apw.

Figura 3.61: Sistemas de reactores con cinética compleja.

Se reliza la simulación como se ha discutido anteriormente, Pasos del 1 al 6, cuidando

que se asigne el modelo termodinámico indicado, NRTL-RK.
Para el Paso 7, comenzaremos con el reactor tipo tanque, TANQUE. Como se muestra
en la …gura 3.62, para este reactor solo es necesario llenar una forma, en donde se especi…cará
las condiciones de operación (temperatura y presión), se deberá especi…car las fases válidas
en la reacción y …nalmente el parámetro de diseño, que para este caso es el volumen del
reactor.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 113

Tabla 3.11: Parametros de diseño de los reactores de cinética compleja.

CSTR
Volumen del tanque 0.14 m3
PFR
Longitud de tubo 2m
Diametro del tubo 0.3 m

Figura 3.62: Formato de especi…caciones del RCSTR.

Para el caso del otro reactor, TUBULAR, se requieren dos formas de datos, la primera
require especi…car las condiciones relacionadas con la temperatura, que pueden ser: operación
adiabática, con per…l de temperatura o con un sólo dato de temperatura conocida, ya que
es un tubo la caída de presión será calculada a travéz del reactor; en el caso del ejemplo
se especi…cará la temperatura a la entrada del reactor, que será la misma que la de la
corriente de la alimentación, ver …gura 3.63a. En una segunda forma (Con…guration) se
deberá capturar la información sobre la geometría del reactor tubular, y las fases válidas en
el mismo, como se observa en la 3.63b.
En los formatos de los reactores, existe una carpeta llamada Reactions, que es la que
permite relacionar los sistemas reactivos con los equipos, ya que como se comenta con ante-
rioridad se de…nen de forma independiente a los reactores.
Para las reacciones se pulsa Reaction desde el panel de navegación, el cual desplegará el
administrador de reacciones; ahí se pulsará New, y se la dará un nombre al sistema reacción,
R-1, y se selecciona el tipo de reacción, en este caso POWERLAW, como se muestra en
la …gura 3.64.
114 CAPíTULO 3 SIMULACIÓN COMERCIAL.

a) b)

Figura 3.63: Formatos para el RPlug a) Especi…caciones y b) Con…guración.

Figura 3.64: Formato de Editor de reacciones complejas.

Al cerrar este formato y pulsar Next se abrirá el editor de reacciones, en él se deberá

escribir las especies involucradas, su coe…ciente estequiomético y el orden en la ley de ve-
locidad para cada una de las reacciones. Para este sistema reactivo la información obtenida
a partir de la reacción química, ecuación 3.4 y las leyes de velocidad, ecuaciónes 3.5 y 3.6.
Los datos generados a partir de esta ecuación se muestran en la tabla 3.12.
Por lo que con esta información y los valores númericos de las leyes de velocidad se
de…nirán dos reacciones, en cada una de ellas se trabajará de forma similar, primero se
de…nirá la reacción directa como se muestra en la …gura 3.65a, en donde se observa que en
los productos no es necesario indicar que el orden parcial, o exponente, es cero; de la misma
forma se muestra el formato para la reacción indirecta, ver …gura 3.65b. La última etapa en
la de…nición de los sistema reactivos es llenar los formatos relacionados con los valores de la
Ley de Arrhenius, para cada reacción, como se muestra en la …gura 3.65c y d; en la primer
…gura se observa el especto que deberá tener el formato para la reacción directa, teniendo en
cuenta que la reacción ocurre en fase líquida, mientrasque en la segunda es la información
correspondiente a la reacción inversa.
En este caso hay que tener en cuenta que no es lo mismo la cinética química de origen
que la operación del reactor, para este ejemplo los reactores validan un sistema bifásico
(liquido-vapor) pero los datos cinéticos fueron obtenidos en una sola fase (líquida).
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 115

Tabla 3.12: Coe…cientes (estequioméricos) y exponentes (orden de velocidad) para el sistema

reactivo de cinética compleja.
Reacción directa
CH3 CH2 OH CH3 COOH CH3 COOCH2 CH3 H2 O
Coe…ciente 1 1 1 1
Orden parcial 1 1 0 0
Reacción indirecta
Coe…ciente 1 1 1 1
Orden parcial 0 0 1 1

a) b)

c) d)

Figura 3.65: Formatos para el editor estequimétrico a) Reacción directa y b) Reacción indi-
recta. Formatos para los datos cinéticos c) Reacción directa y d) Reacción indirecta.
116 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Finalmente se ejecuta la simulación y se obtiene los resultados que se muestran en la

…gura 3.66.

Figura 3.66: Resultados de los reactores de cinética compleja.

Reactores con formas cinéticas tipo LHHW. Una teoría ampliamente aceptada para
los sistemas catalítico heterogéneos, donde el catalizador es un sólido poroso, es el propuesto
por Langmuir y Hinshelwood; en este modelo se considera que para que exista una reacción
química al menos uno de los reactivos necesariamente deberá ser adsorbido por el super…cie
del catalizador y que es aquí, en la super…cie catalítica, donde ocurre la generación de los
productos, que aparecerán adsorbidos a la super…cie catalítica; …nalmente los productos
deberán ser desorbidos del catalizador. Los procesos de adsorción y desorción son también
reacciones química, por lo que la expresión de velocidad siempre es mucho más compleja que
las consideradas con expresiones basadas en la Ley de la Potencia.
El simulador Aspen considera este modelo en base a tres términos

(factor cinético) (expresión de fuerza gobernante)

r= (3.7)
(términos adsorción)

1. Factor cinético.- Este el término relacionado con la Ley de Arrhenis, ver ecuación
3.3. La única situación a considerar es que las unidades del factor preexponencial sean
las del sistema internacional (SI).

2. Fuerza gobernante.- Esta relacionado con los términos de concentración de la reac-

ción y genéricamente tiene la forma:

Y
N Y
N
(expresion de fuerza gobernante) = k1 Ci i k2 Cj j (3.8)
i=1 j=1
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 117

donde las k, son constantes experimentales y la Ci i son las expresiones de concentración

elevadas a las potencias experimentales de ajuste.

3. Términos de adsorción.- Son los términos asociados con la adsorción de reactivos

y productos en el mecanismo de reacción, está dado por la expresión general:
!!m
XM YN
(términos adsorción) = Ki Civi (3.9)
i=1 i=1

donde, las Ki son las constantes de adsorción y las Civi son las expresiones de concen-
tración elevadas a los coe…eintes estequiométricos en las reacciones propuestas en el
mecanismo de la reacción.

Todas las constantes, independientemente del término al que se re…era, se expresará en

el simulador como:
B
k =A+ + C ln(T ) + DT (3.10)
T
donde A, B, C y D son los valores que serán capturadas en Aspen, por lo que siempre es
necesario adecuar los valores experimentales para poder ser usadas.
Para ejempli…car este modelo y hacer más comentarios para el uso de esta cinética se
desarrollará el siguiente proceso.

Ejemplo 41 Producción catalitica de metanol a partir de dioxido de carbono. En la síntesis

de metanol, utilizando Oxido de Zinc - Cobre soportado el alumina, ocurren las reacciones
en fase gas:
CO2 + 3H2 CH3 OH + H2 O
CO2 + H2 CO + H2 O
Bussche y Froment has determinado las expresiones de velocidad LHHW (con unidades de
kmol
kgcat s
/) para estas reacciones:

PC H 3 O H PH 2 O
KR1 KD1 PCO 2 PH 2 KI1 PH2
2
r1 = 3
P
1+ K1 PH22 O + K2 PH0;5
2
+ K3 P H 2 O
H 2
PC O PH 2 O
KR2 KD2 PCO 2 KI2 PH 2
r2 = PH 2 O
1 + K1 PH2
+ K2 PH0;5
2
+ K3 P H 2 O
2

donde los valores de las constantes de las ecuaciones se muestran en la tabla 3.13. Se desea
hacer una simulación de este sistema, para lo cual se alimentará una corriente a 250 C y 110
bares de presión, la cual contiene 4,066 kmol/hr de monóxido de carbono, 28,920 kmol/hr
118 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Tabla 3.13: Valores de las constantes asociadas al mecanismo de velocidad para la síntesis
de metanol.
4413;76
KR1 1;07 10 13 e( T )

KD1 1
7059;726
KI1 3;9084 1020 e( T )
11398;244
KR2 122e( T )
4773;259
KI2 0;009354e( T )
K1 3453;38
2068;44
K2 1;578 10 3 e( T )
14928;915
K3 6;62 10 16 e( T )

de hidrógeno y 3,976 kmol/hr de dióxido de carbono; utilizando un PFR multitubular de

8000 unidades, de 12.2 m de largo y con un diámetro de 3.75 cm, considere un per…l de
temperatura a lo largo de los tubos, en la entrada del reactor la temperatura se de 150 C,
al 90 % de la distancia de los tubos la temperatura presentrá un máximo de 280 C y a la
salida disminuirá a 267 C. El proceso se muestra en la …gura 3.67. Al términar grabar como
actividad12.apw.

Figura 3.67: Reactor tubular con modelo cinético LHHW.

Se inicia esta simulación de la forma que se ha discutido anteriormente, Pasos del 1 al

6, se recomienda para este tipo de sistemas reactivos con catalizadores utilizar los métodos
termodinámicos SR-POLAR, SRK o PSRK. En este caso se utilizará SRK.
En este caso la con…guración del reactor tubular deberá se requieren tres formas de
datos, la primera require especi…car las condiciones del per…l de temperatura a partir de
la infromación dada, como se observa en la 3.68a, la segunda que esta relacionada con la
geometría del reactor y la fase válida en la simulación, ver …gura 3.68b; la tercera, que esta
relacionada por la presencia del catalizador, en la …gura 3.69 se muestra el aspecto de esta
forma.
Para de…nir el sistema reactivos, se procede como en el ejercicio anterior, se selecciona
Reaction desde el panel de navegación, el cual desplegará el administrador de reacciones; ahí
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 119

a) b)

Figura 3.68: Formatos para el RPlug a) Especi…caciones y b) Con…guración.

Figura 3.69: Formato para el RPlug para de…nir las características del catalizador .

se pulsará New, y se la dará un nombre al sistema reacción, R-1, y se selecciona el tipo de

reacción, en este caso LHHW. Como se tiene dos reacciones, deberán ser capturadas de for-
ma independiente; para este modelo sólo se deberá ser de…nidos coe…cientes estequimétricos
de las reacciones involucradas, como se muestra en la …gura 3.70.

a) b)

Figura 3.70: Editor de reacciones a) reacción principal y b) reacción secundaria.

Para de…nir los términos del moldelo LHHW, se explicará la forma en que deben de ser
adaptados los datos experimentales a las formas de Aspen Plus. Inicialmente se simpli…carán
120 CAPíTULO 3 SIMULACIÓN COMERCIAL.

las leyes de velocidad, ecuaciones 3.5 y 3.6, que pueden tomar las siguientes formas:

PC H 3 O H PH 2 O
K1 KD1 PCO 2 PH 2 K1 KI1 PH2
2
r1 = 3 (3.11)
P
1+ K1 PH22 O + K2 PH0;5
2
+ K3 P H 2 O
H 2

PC O PH 2 O
K2 KD2 PCO 2 K2 KI2 PH 2
r2 = PH 2 O
(3.12)
1 + K1 PH2
+ K2 PH0;5
2
+ K3 P H 2 O
2

A partir de estas nuevas expresiones de velocidad se puede escribir los términos:

1. Factor cinético.- Al realizar la simpli…cación en las ecuaciones, los valores de K0 y E

para cada factor cinético se reducen a su expresión más sencillas, para las ecuaciones
3.11 y 3.12, los valores para este término se muestran en la …gura 3.71.

a) b)

Figura 3.71: Formatos para el factor cinéticos a) reacción principal y b) reacción secundaria.

2. Fuerza gobernante.- Para determinar la información requerida, se deberá analizará

cada ecuación cinética; en general ambas ecuaciones tiene dos términos y dos cons-
tantes. Para el caso de la primera, ecuación 3.11, se tendrá que los valores asociados
a las constantes en términos de los requisitos del simulador, se muestran en la tabla
3.14.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 121

Tabla 3.14: Valores de las constantes (formato requerido en Aspen) para el término de fuerza
gobernante para la reacción principal.
B
k =A+ T + C ln(T ) + DT
Constante valor en la ecuación A B C D
4413;76
K1 KD1 1;07 10 13 e( T ) -29.866 4413.76 0 0
2645;966
K1 KI1 4;182 107 e( T ) 17.549 -2645.966 0 0

Mientras que para los términos de las concentraciones de la misma rección, que están
expresadas en función de la presión parcial, se muestran en la tabla 3.15.

Tabla 3.15: Valores para los términos de concentraciones para la fuerza gobernante de la
reacción principal.
Término 1 Término 2
Componente Reactivos Productos Reactivos Productos
CO2 1 0
H2 1 -2
CH3 OH 0 1
H2 O 0 1

Con esta información se llena el formato para la primera reacción, para el Driving
Force, como se ve en la …gura 3.72.

Utilizando la misma metodología para la segunda reacción, ecuación 3.12, la informa-

ción …nal requerida en para los formatos de Driving Force, valores de las constantes
y los exponentes de cada uno de los términos, se muestran en la …gura 3.73

3. Términos de adsorción.- Para este término se deben expresar las constantes de

adsorción en la forma de requerida por el simulador y la relación de presiones parciales
para cada una de las especies. En este caso en las reacciones 3.11 y 3.12, tiene la misma
expresión, solo el exponente es diferente; para este caso los valores de las constantes
serán las mostradas en la tabla 3.16.

Para determinar los componentes asociados a cada una de las constantes se analizan
cada uno de ellos, se deberá considerar que el primer término del polinomio, el valor
1, también se incluye, los resultados se muestran en la tabla 3.17.
122 CAPíTULO 3 SIMULACIÓN COMERCIAL.

a) b)

Figura 3.72: Formas para la fuerza gobernante para la reacción principal a) reacción directa
(Term1) y b) reacción inversa (term 2).

a) b)

Figura 3.73: Formas para la fuerza gobernante para la reacción secundaria a) reacción directa
(Term1) y b) reacción inversa (term 2).

Para determinar los componentes asociados a cada una de las constantes se analizan
cada uno de ellos, se deberá considerar que el primer término del polinomio, el valor
1, también se incluye, los resultados se muestran en la tabla 3.17.

Con la información anterior se llena la forma Adsorption Expression, que para la

primera reacción, ecuación 3.11 se muestra en la …gura 3.74a, considere que además de
esta información deberá incluirse el exponente al cual está elevado el polinomio, que
en este caso es tres.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 123

Tabla 3.16: Valores de las constantes (formato requerido en Aspen) para el término de
adsorción para la reacción principal y secundaria.
B
k =A+ T + C ln(T ) + DT
Constante valor en la ecuación A B C D
K1 3453;38 8.14711 0.0 0 0
2068;44
K2 1;578 10 3 e( T ) -6.45416 2068.44 0 0
14928;915
K3 6;62 10 16 e( T ) -34.9513 14928.9 0 0

Tabla 3.17: Valores para los términos de concentraciones para el factor de adsorción para la
reacción principal y secundaria.
Componente Término 1 Término 2 Término 3 Término 4
H2 0 -1 0.5 0
H2 O 0 1 0 1

a) b)

Figura 3.74: Forma para el término de adsorción a) para la reacción principal b) para la
reacción secundaria.

Como se discutió anteriormente, la expresión de adsorción es la misma para la segunda

reacción, ecuación 3.12, con la única diferencia que el exponente de la expresión es uno, en
la …gura 3.74b se muestra la forma para los términos de adsorción de esta reacción.
Terminada esta captura se puede resolver la simulación, los resultados de las corrientes
invlocradas en el reactor catalítico se muestran en la …gura 3.75.
124 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.75: Resultados del reactor de ‡ujo con cinética tipo LHHW.

3.1.8. Separadores.
Aún cuando los reactores son considerados como el corazón de una planta química, los
separadores son las etapas más importantes en los procesos, ya que casí nunca los sistemas
reactivos tienen conversiones cercanas al 100 %, por lo que es necesario puri…car las corrientes
que salen de los reactores, para separar los productos, subproductos y materia prima sin
reaccionar. El simulador Aspen tiene un conjunto muy amplio de este tipo de procesos, que
pueden ser clasi…cados en dos tipos: a)separadores elementales y de una sola etapa; y b)
separadores multietapas.
Los primeros involucran:
1. Separadores genéricos, basados unicamente en los balances de materia, y que no re-
quieren mayor información de diseño.
2. Destiladores ‡ash, de dos y tres fases, como se muestra en la …gura 3.76.

Figura 3.76: Separadores. Paleta de diseño.

Los segundos representan a columnas o torres de multiples etapas, con ellos se pueden
simular columnas de destilación, torres de adsorción, torres de extracción, entre otras opera-
ciones, en la …gura 3.77 se muestran los modelos disponibles en el simulador.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 125

Figura 3.77: Columna. Paleta de diseño.

Separadores elementales.
El modelo más simple que presenta el simulador para este tipo de procesos son los bloques
Sep y Sep2, en los que unicamente se requiere especi…car la cantidad, másica o molar, o
fracción de los componentes en las corrientes de salida. La diferencia entre estos modelos,
consiste en el número de corrientes de salida, en el Sep se pueden incluir todas las corrientes
necesarias, mientras que Sep2 sólo pueden tener dos.

Ejemplo 42 En destilador que opera a 1 atmósfera, se alimentan dos corrientes, como se

muestra en la …gura 3.78. El ‡ujo de la corriente FEED1, el cual es un líquido saturado a
14.696 psi, es de 50 lbmol /hr de metanol y 100 lbmol /hr de agua. Para la corriente FEED2,
se tiene un ‡ujo de 150 lbmol /hr de n-hexano puro, a 120 C y 14.696 psi. Las especi…caciones
de los componentes a la salida son los siguientes: a) La mitad del metanol se obtiene en la
corriente superior; b) 99 % del agua en la alimentación en la corriente superior y c) 99 %
del n-hexano de la alimentación se recupera en el fondo de la columna. Al términar grabar
como actividad13.apw.

Figura 3.78: Sistema de separación utilizando Sep2.

Se realiza la simulación como se ha discutido anteriormente, Pasos del 1 al 6. Para el

Paso 7, es necesario considerar lo siguiente; ya que este modelo no requiere las condiciones de
operación (presión y temperatura del separador) ya que la solución se calcula utilizando los
balances de materia, y se considera que las corrientes de salida tienen las mismas consiciones
que las de la alimentación. Por lo que solo es necesario de…nir las especi…caciones para cada
componente. En este ejercicio hay condiciones para la corriente superior e inferior, para
la superioir se especi…ca la fracción requerida y la unidad adecuada, como se muestra en
126 CAPíTULO 3 SIMULACIÓN COMERCIAL.

la …gura 3.79a, mientras que para la unica especi…cación de la corriente inferior se deberá
llenar la forma como se muestra en la …gura 3.79a.

a)

b)

Figura 3.79: Formatos de especi…caciones para el Sep2 para a) corrienteDESTIL y b) corriente

BOT.

Al terminar este ejercicio se encontrán los resultados que se presentan en la …gura 3.80.

Figura 3.80: Resultados para el sistema de separación utilizando Sep2.

Destilación Flash.
Como se comento previamente, una de las operaciones iniciales de separacion son las
separaciones ‡ash, las cuales se emplean para separar componentes que presentan valores
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 127

muy distintos de volatilidad, expresado de otra forma, que presentan valores de presión de
vapor muy diferentes; el simulador presenta dos tipos de cámaras de ‡ash, la que generan dos
fases (vapor - líquido) Flash2 y la que genera tres fases (vapor - líquido -líquido) Flash3, y
en general se con…guran de igual manera.

Ejemplo 43 A una cámara de ‡ash se alimentan dos corrientes, tal como se muestra en
la …gura 3.81. Las alimentaciones se encuentran a a 70 C y 760 mm de Hg de presión. La
cámara opera en condiciones adibáticas a 760 mm de Hg de presión. La primera corriente
FEED1, tiene un ‡ujo de 40 lbmol /hr de etanol y 10 lbmol /hr de acetona, mientras que
la corriente FEED2 esta constituida solo por 50 lbmol /hr de agua. Realice la simulación
utilizando un modelo Flash2. Al términar grabe sus resultados como Actividad14.apw.

Figura 3.81: Sistema de separación, utilizando Flash2.

Se realiza la simulación, como en las ocasiones anteriores, para el caso del Flash2, se
deberá llenar la forma que se muestra en la …gura 3.82. La operación adiabática se especi…ca
utilizando la opción Duty, que se re…ere a la cantidad de calor que se suministra si se le
asigna un cero se tendrá la condición indicada.

Figura 3.82: Formato para el Flash2.

Al resolver la simulación, se encuentran los resultados que se muestran en la …gura 3.83.

128 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.83: Resultados del sistema de separación utilizando Flash2.

Destilación multietapas.
La simulación de columnas de destilación es una de las herramientas más amplias de-
sarrolladas en el simulador, ya que a diferencia de la mayoría de los procesos simulados en
Aspen, en estos módulos se puede hacer el diseño inicial del proceso, como el módulo DST-
WU, y la simulación rigurosa, módulo RadFrac, en la que se involucra la e…ciencia de los
platos, la con…guración, tipo, entre otras más.
El primer módulo, el DSTWU, diseña la columna de destilación a partir del método
corto de Fenske-Underwood-Gilliland (FUG), en el Ápendice B se presenta un resumen
del método; como resultados de esta simulación se obtiene el número de etapas, el re‡ujo
de operacion y la ubicación del plato de alimentación. Uno de los concepto fundamentales
en el que se basa el metodo de FUG, es del componente clave, que es quién determina
donde se realiza la separación, como se tiene dos corrientes de salida deberá especi…carse
dos componentes clave. De estos dos, al más volátil se le da el nombre de componente clave
ligero (LK) y al menos volátil, componente clave pesado (HK).
La cantidad esperada de cada uno de ellos en la corriente de destilado o fondos deberá
especi…carse inicialmente, este valor se expresa como una fracción (fLK o fHK ), a partir de
las siguientes expresiones:

‡ujo molar de LK en el destilado

fLK = (3.13)
‡ujo molar de LK en la alimentación

‡ujo molar de HK en el destilado

fHK = (3.14)
‡ujo molar de HK en la alimentación

Ejemplo 44 Se desea diseñar una columna de destilación para separar una corriente de
hidrocarburos, en la tabla 3.18 se muestra los componentes de la mezcla y sus ‡ujo molares.
Utilice una columna DSTWU para este proceso, como se muestra en la …gura 3.84; las
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 129

condiciones de la alimentación son 78.4 psia y 243.13 F, mientras que la columna opera a
70.4 psia en la cámara de destilado, y tiene una caída de presión en la columna de 10 psi.
Considere que el re‡ujo de operación es 1.2 veces el re‡ujo mínimo. Se desea recuperar 442
lbmol /hr de n-butano (componente clave ligero), en la corriente de destilado, y 23 lbmol /hr de
isopentano (componente clave pesado) en l corriente de fondos. Si la columna está equipada
con un condensador total, determine el número de etapas y el valor de re‡ujo de operación.
Al términar grabe sus resultados como Actividad15.awp.

Tabla 3.18: Composición de la corriente de alimentación al sistema de separacion DSTWU.

Componente ‡ujo, lbmol /hr
isobutano 12.0
n-butano 448.0
isopentano 36.0
n-pentano 15.0
n-hexano 23.0
n-heptano 39.1
n-octano 272.2
n-nonano 31.0

Figura 3.84: Sistema de separación DSTUW.

Se elabora el simulador de la forma conocida, para la llenar la forma del modelo DSTWU,
se deben conocer las fracciones de separación para cada uno de los componentes clave,
ecuaciones 3.13 y 3.14, por lo que se puede estimar:
442
fLK = = 0;9866
448
130 CAPíTULO 3 SIMULACIÓN COMERCIAL.

36 23
fLK = = 0;6389
36
Para el caso de la relación de re‡ujo, en el simuladar se debe indicar con un signo negativo,
la realación entre en re‡ujo de operación y el re‡ujo mínimo, para este caso como la relación
es de 1.2 veces, se deberá escribir como -1.2. Los demás valores se capturan en la forma, la
cual deberá tener el aspecto que se muestra en la …gura 3.85.

Figura 3.85: Formato Especi…caciones para la columna DSTWU.

Al concluir la simulación se pueden encontrar los resultados de las corrientes de salida,

como se muestra en la …gura 3.86a. Mientras que los resultados de diseño se muestran en la
…gura 3.86b. En este caso la columna que realiza la separación deberá tener aproximadamente
15.9 etapas, el plato de alimentación localizado en la etapa 4.88 y un valor de 0.622 como
razón de re‡ujo de operación.

a) b)

Figura 3.86: Resultados del sistema de separación DSTWU a) ‡ujos de salidas y b) valores
de diseño.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 131

Otro de los módulos de separación que puede ser utilizado para distintos tipos de sepa-
ración (columnas de destilación, torres de absorción, torres de extracción, entre otras) es el
RadFrac. Este modelo resuleve las ecuaciones del balance de materia, de energía y de equi-
librio, conocidas como ecuaciones MESH; lo cual permite conocer los per…les temperatura,
presión, composición, ‡ujos internos y equilibrio a los largo del equipo de separación, para
utilizar este modulo es necesario conocer el número de etapas, la ubicación de las alimenta-
ciones y los valores de la presión en la parte superior e inferior del separador, a continuación
se desarrollarán cuatro simulaciones que permitirán observar como se utiliza este modulo.

Ejemplo 45 Considere una columna de destilación que está con…gurada con 17 etapas teóri-
cas, una alimentación en la etapa 8, que posee un condesador y un rehervidor parcial. A esta
columna se alimenta la corriente FEED, que se encuentra a 213.9 F y 200 psia de presión,
los ‡ujos molares de los componentes en la corriente de alimentación se muestra en la tabla
3.19. La columna opera de tal forma que se produce un ‡ujo molar de la corriente de destila-
do de 23 lbmol /hr con una velocidad de re‡ujo de 150 lbmol /hr. Si la columna se encuentra a
una presión constante de 250 psia, cuales serán las composiciones de las corrientes de salida
y cómo serán los per…les de composición, temperatura y ‡ujos internos. La con…guración
de la columna se muestra en la …gura 3.87, observe la posición de la corriente DESTIL. Al
terminar guarde sus resultados como Actividad16.apw.

Tabla 3.19: Composición de la corriente de alimentación al sistema de separacion RadFrac.

Componente ‡ujo, lbmol /hr
etano 3.0
propano 20.0
n-butano 37.0
n-pentano 35.0
n-hexano 5.0

Al desarrollar este simulador con los pasos ya comentados con anterioridad, se debe de
tener encuenta que al indicar que el condensador es parcial se debe de utilizar la conección
adecuada; en la …gura 3.88 se muestran las opciones para este modelo; la conección A
corresponde a un condensador parcial, es decir que la corriente del destilado puede ser de
vapor saturado o una mezcla en equilibrio de líquido y vapor; la B corresponde a una
salida de agua del condensador; y C es la opción que debe tomarse para el caso de utilizar
un condensador total, el destilado solo esta constituido por líquido saturado. Solo puede
selccionarse una de las tres y deberá corresponder al diseño deseado. Si se va a utilizar un
condensador parcial, como en el caso del ejercicio, se deberá seleccionar desde el momento
de elaborar el diagrama de ‡ujo.
132 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.87: Sistema de separacion RadFrac.

Figura 3.88: Con…guración del condensador de la columna RadFrac.

Una vez que se ha seleccionado el tipo del condensador de forma adecuada, se comienza
con el llenado de las formas; los datos requeridos por el modelo se divide en tres formatos;
el primero, ver …gura 3.89, corresponde a la con…guración general de la columna, número de
etapas, tipo de condensador y rehervidor, y dos especi…caciones de diseño, el ‡ujo de destilado
y la razón del re‡ujo. En el caso de tener un condensador parcial se deberá especi…car si solo
hay vapor (Partial-Vapor) o tambien hay líquido (Partial-Vapor-Liquido).
El segundo formato está relacionado con la ubicación de las corrientes de la columna,
tanto las alimentaciones como las corrientes de salida; teniendo en cuenta que el simulador
cuenta las etapas de arriba hacia abajo, siendo la número uno, el condensador parcial (si lo
hay) y la última (etapa N) el rehervidor. En la …gura 3.90a se muestra como debe quedar
este formato para la simulación del ejemplo.
El ultimo formato corresponde a las condiciones de presión dentro de la columna, ya que
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 133

Figura 3.89: Formato para la con…guración de la columna RadFrac.

se puede de…nir una operación constante o un per…l general o especí…co de presión. Si solo
se especi…ca una sola presión (Top stage) la columna operará a presión constante, como
es el caso del ejemplo; si se especi…ca además la presión en la etapa 2 (Stage 2 pressure) se
indicará que la columna y el condensador operan a presiones distintas, la primera será la del
condensador y la segunda la de la columna; …nalmente, si se especifíca alguna de las caídas
de presión (Pressure drop for rest of column), ya sea por plato o por toda la columna, se
establecerá una presión en la columna, y un per…l de presiones desde el domo de la columna
hasta sus fondos. Para el ejemplo de trabajo, lo forma quedará como se muestra en la …gura
3.90b.

a) b)

Figura 3.90: Formato para la columna RadFrac a) Corrientes y b) Presión.

Los ‡ujos de salida de la columna al resolver la simulación se muestran en la …gura

3.91. En la …gura se observa las concentraciones de salida de las dos corrientes y se observa
que la corriente de destilado contiene principalmente etano y propano, generalmente a los
productos de la corriente de destilación se les suele llamar ligeros (por el método de diseño
anteriormente discutido), mientras que el la corriente de fondos, o pesados, se encuentran
los restantes componentes.
Sin embargo lo que es de interés para una columna de destilación, resuelta por un método
riguroso, es poder conocer los per…les de las distintas variables, lo cual permitirá determinar
134 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.91: Resultados del sistema de separación RadFrac.

el diametro, la altura de la columna, y el diseño hidráulico de la misma. Al seleccionar el

bloque de la RadFrac, es posible generar de forma sencillas los diferentes per…les. En la
…gura 3.92 se muestran los per…les de composición molar de la fase líquida para todos los
componentes en cada una de las etapas de la columna.

Figura 3.92: Per…l de composición molar para el sistema de separación RadFrac.

En la …gura 3.93 se muestran los per…les de la temperatura y los ‡ujos molares internos;
estos valores en los métodos de diseño más simpli…cados de diseño, como el McCabe-Thiele se
consideran constantes, por lo que una solución rigurosa proporciona información más precisa
en el diseño de las columnas.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 135

a) b)

Figura 3.93: Per…les para el sistema de separación RadFrac a) Temperatura y b) Flujo

molares.

Funciones adicionales del modulo RadFrac. El módulo de RedFrac incluye varias

rutinas adicionales a la simulación multietapas, ver …gura 3.94, entre ellas podemos consid-
erar:

Figura 3.94: Funciones adicionales de la columna RadFrac.

1. Especi…cación de variables de diseño..- Es decir que se puede especi…car uno o

más valores en las corrientes de salida, como lo son el ‡ujo, la composición, el calor
adicionado o removido, entre otras (Desing Speci…cations - Vary).
2. Especi…cación de la e…ciencia de la etapa teórica.- Se pueden especi…car esta
característica, ya sea para toda la columna o por secciones (E¢ ciencies).
136 CAPíTULO 3 SIMULACIÓN COMERCIAL.

3. Especi…cación del tipo de plato o empaque.- Se puede selaccionar el tipo de plato

(cachuchas, valvulas, plato perforado) o el empaque (estructurado o al azar), en toda
la columna o en una sección (Tray Sizing).

Estas funciones solo proporcionarán una mayor información sobre la columna, se re-
comienda resolver inicialmente la columna sin utilizar alguna de estas rutinas; a excepción
de que se especi…que lo contrario, por variables de diseño, las composiicones de las corrientes
de salida no cambiarán, el unico efecto sobre la solución original es un cambio en la cantidad
de energía involucrada, ya sea en el condensador o en el rehervidor. Para demostrar estas
características se plantea el siguiente ejercicio, la columna es la resuelta en el caso anterior
y solo se revisarán las funciones adicionales.

Ejemplo 46 Considere la columna RadFrac que se resolvió con anterioridad, se desea re-
solver nuevamente bajo las siguientes características:

1. Que la corriente de salida del destilado tenga una composición molar 0.85 para el
propano. Al …nalizar este apartado grabe sus resultados como Actividad17.apw.

2. Evaluar el desempeño térmico de la columna si, junto con la condición de la salida, se

considera además que la e…ciencia de la etapa es del 75 % en los platos de la columna
y 90 % en los equipos de transferencia de calor.Al …nalizar este apartado grabe sus
resultados como Actividad18.apw.

3. Determinar el per…l del diámtero de la columna, resuelta en las condiciones anteriores,

si la columna está equipa con platos con cachuchas. Al …nalizar este apartado grabe sus
resultados como Actividad19.apw.

Apartado 1.
Para poder especi…car la fracción molar de un componente, es necesario seleccionar la
opción Desing Speci…cations de la carpeta Speci…cations. Se abrirá un administrador
y se seleccionará la opción New, con lo que se desplegará el formato correspondiente, que
tiene tres secciones. En la primera, …gura 3.95, se de…ne la especi…cación, para el caso del
ejemplo se selcciona Mole purity en Type, y se establece el valor requerido en Target.
En la segunda forma y tercera forma, se deberá especi…carse cual es el componente que
se desea separar, ver …gura 3.96a, y cual es la corriente en la que se encuentra, ver …gura
3.96b.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 137

Figura 3.95: Formato Especi…cación de la herramienta Desing Speci…catios, columna Rad-

Frac.

a) b)

Figura 3.96: Formatos de la herramienta Desing Speci…cations, columna RadFrac a) Com-

ponentes y b) Corrientes.

Cuando se utiliza un Desing Speci…cations necesariamente deberá estar ligado a un

Vary, ver 3.94, que es la forma en la que se indica cual de las variables de la columna
deberá ser modi…cada para cumplir con la especi…cación. Se selcciona la carpeta Vary, y
en el administrador que se desplegará pulsar New. Se desplegara el unico formato de esta
opción, para seleccionar la variable a modi…car (Type) es necesario conocer el efecto que
puede provocar en las demás variables de la columna, en general las dos variables a utilizar
son la velocidad del re‡ujo o el ‡ujo de destilado; en este caso se modi…cará la velocidad
del re‡ujo, para mantener constante el ‡ujo de destilado, sin embargo se esperaría que el
diamtero de la columna aumente y que los gastos de calor cambien. En la …gura 3.97 se
muestra este formato, en este caso se especi…ca un rango amplio de valores para variar el
valor de la razón de ‡ujo.
Al realizar estas modi…caciones en el problema anterior, el simulador presentará una
leyenda indicando que los datos han sido modi…cados, por lo que es posible volver a resolver
138 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.97: Formato Especi…cación de la herramienta Vary, columna RadFrac.

el sistema. Al hacerlo, se encuentran los resultados que se muestran en la …gura 3.98, donde
se observa que se cumplió con la especi…cación requerida.

Figura 3.98: Resultados de la herramienta Desing-Vary, columna RadFrac.

Apartado 2.
Como ya se comento con anterioridad, en el simulador Aspen, el conteo de las etapas se
hace de arriba a bajo, la etapa uno corresponde al condensador, siempre y cuandos sea un
condensador de tipo parcial, mientras que la última la corresponde al rehervidor.
Para especi…car la e…ciencia de platos se selecciona la carpeta de E¢ ciencies, la cual
despliega la forma que se muestra en la …gura 3.99a. Se deberá seleccionar la e…ciencia de
Murphree y la posibilidad de de…nir esta característica por secciones. Al realizar esta acción
se deberá llenar el formato de Vapor-Liquid, en ella, la columna se dividirá en secciones
de arriba hacía abajo, se deberá especi…car la etapa inicial, la …nal y la e…ciencia de cada
sección, como se muestra en la …gura 3.99b.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 139

a) b)

Figura 3.99: Formato de la herramienta E……ciencies, columna RadFrac a) Opciones y b)

Vapor-Líquido

Al resolver esta columna, los resultados de las concentraciones resultantes y los ‡ujos son
iguales a los anteriores; los efectos de este cambio, los valores de los consumos de energía se
muestran a continuación en la tabla 3.20.

Tabla 3.20: Consumo de energía en el condensador y rehervidor para un sistema de destilación

resuelto con diferentes herramientas.
Diseño Condensador, BTU/lb 10 6 Rehervidor, BTU/lb 10 6

Original 0;8585 1;0346

Con especi…cación 0;8837 1;0559
Con e…ciencia 2;1612 2;3441

Apartado 3.
Finalmente para determinar el diámetro de la columna, se selecciona la carpeta Tray
Sizing, con lo que despleglará la forma respectiva, en donde solo hay que especi…car la
sección de la columna y el tipo de etapa, ver …gura 3.100.
Al resolver esta simulación se encuentran los per…les de las características asociadas al
diseño hidráulico de la columna, el diametro el área activa, etc., que se muestran en la …gura
3.101a. El simulador puede generar las grá…cas de los per…les de esta sección, por lo que
el per…l del diámetro se muestra en la …gura 3.101b, en donde se observa que en la parte
inferior de la columna se encuentra el diámetro máximo.
140 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.100: Formato Especi…caciones de la herramienta Tray Sizing, columna RadFrac.

a) b)

Figura 3.101: Resultados de la herramienta Tray Sizing, columna RadFrac a) Per…les de las
variables y b) Per…l del diámetro de la columna.

Destilación reactiva. Una de las usos más importantes del modelo RadFrac, es que puede
ser utilizada para simular columnas reactivas, La destilación reactiva es un proceso que
integra la reacción y la separación en un mismo equipo. Esta integración tiene diversas
ventajas potenciales, como son: reducción en los costos, mejoras en la conversión de los
reactivos y en la selectividad de los productos, así como ahorros en el uso de energía y
reactivos entre otros. En el siguiente ejemplo se discutirán los aspectos principales de este
proceso.
Ejemplo 47 El proceso de producción de MTBE (Metil TerButil Éter) es una de los primeros
que ha sido adaptado a las columnas reactivas, como se muestra en la …gura 3.102. En este
caso se debe tener información sobre a cinética y sobre las carcaterísticas de la columna. La
reacción química en fase líquida que ocurre dentro de la columna es:
i-buteno + metanol MTBE
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 141

y ocurre de acuerdo a las siguientes expresiones de velocidad:

92440 xiB
rD = 3;67 1012 e
RT xM eOH

134454 xM T BE
rI = 2;67 1017 e
RT x2M eOH
Se dispone de dos corrientes de alimentación que entraran a la columna de destilación, como
se muestra en la …gura 3.102, la corriente METANOL tiene un ‡ujo de 711.3 kmol/hr y
se encuentra a 47.4 C y 11.652 bar de presión, en esta corriente únicamente se suministra
metanol; en la corriente BUTENOS se alimentan 1965.8 kmol/hr de una mezcla del 0.36
porciento molar de isobuteno y el resto de n-buteno, esta corriente se encientra a 180.4
C y 11.885 bar de presión. La columna que se empleará esta equipada con 17 etapas de
equilibrio, un condensador total y un rehervidor, el MTBE se espera que se recupere por la
corriente inferior del equipo, por lo que se espera un ‡ujo de 640 kmol/hr de la corriente
BOTTOMS, además de que la columna operarará con una razón de re‡ujo igual a 7. La
corriente de metanol se alimentará en el plato 10 y la de butenos en el 11. La columna estará
empacada desde la etapa 2 hasta la 16, utilizando en empaque Mellapack Sulzer 250Y, el
di{ametro de la columna es de 6 m y tiene una HETP = 1 m. La zona reactiva de la columna
estará comprendidad entre los platos 4 al 11, utiizando 1000 kilogramos de catalizador en
cada una de estas etapas. Se desea simular la columna utilizando un modulo RadFrac y
conocer las composiciones de las corrientes de salida, así como los per…les de composición.

Figura 3.102: Columna Reactiva, sistema de separación RadFrac.

La simulación de este proceso se lleva de la misma forma en que se resuelve la columna

de RadFrac, como el ejemplo anterior; sin embargo en la forma de Con…guration, deberá
considerarse seleccionar Strongly non-ideal liquid en la opción de convergencia, ver …gura
3.103a, los datos de la ubicación de las alimentaciones se deberá introducir en el formato de
Stream, como se muestra en la …gura 3.103b.
142 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.103: Formatos para la columna reactiva RadFrac. a) Con…guración y b) Corrientes.

Para terminar la con…guración básica de la columna, se deberá establecer la presión de

la columna, y en este caso la caída de presión a travéz de ella, como se muestra en la …gura
3.102.

Figura 3.104: Formato Presión para la columna reactiva RadFrac.

Hasta aquí termina la con…guración básica de la columna, para continuar es necesario

abrir la carpeta Speci…cations>Reactions, ver …gura 3.94. Con esto se muestra la forma
que se muestra en la …gura 3.105, en la cual en la primera parte se deberá especi…car las
etapas en que se va a llevar a cabo la reacción y se de…ne, además, la reacción que estará
ocurriendo, 3.105a. En la segunda de estas formas se deberá especi…car las etapas en donde
está presente el catalizador, no necesariamente la misma que la anterior, y la cantidad total
de catalizador, ver …gura 3.105b.
3.1 MANEJO DE UN SIMULADOR COMERCIAL. 143

a) b)

Figura 3.105: Formatos de la herramienta Reaction, columna rectiva RadFrac. a) Especi…-

caciones y b) Catalizador.

Para terminar de con…gurar la columna, se selecciona la carpeta Sizing and Rating>Pac-

king Rating, que es la forma en al que se puede con…gurar el empaque de una columna, con
esta cción aparecerá un administrador en el cual se deberá pulsar New para que aparezca la
foma que se muestra en la …gura 3.106. En ella se tiene que especi…car la sección empacada,
el tipo de empaque y su demás carcaterísticas.

Figura 3.106: Formato Especi…cación de la herramienta Packing Rating columna rectiva

RadFrac.

Finalmente se de…ne el sistema reactivo, utilizando el procedimiento para una cinética

tipo PowerLaw, los valores de las leyes de velocidad para cada reacción se muestra en la
…gura 3.107, se debe considerar que la reacción ocurre en fase líquida y que la ley de velocidad
está en términos de fracciones molares.

Los valores de las corrientes de salida, se muestran en la …gura 3.108a, mientras que los
per…les de composición de la fase líquida a travéz de la columna están en la …gura 3.108b.
144 CAPíTULO 3 SIMULACIÓN COMERCIAL.

a) b)

Figura 3.107: Formatos de los datos para las reacciones para el MTBE. a) reacción directa
y b) reacción indirecta.

a) b)

Figura 3.108: Resultados para la columna reactiva. a) Valores de las corrientes y b) Per…l de
composición dentro de la columna.

3.2. Solución de casos de estudio de diseño, análisis y

optimización de procesos en ingeniería química.
El desarrollo de la simulación de plantas de procesos se ejempli…ca a continuación,se
utiliza un ejemplo desarrollado por Luyben en el texto "Principles and case studies of Si-
multaneus desing", para la producción de cumeno, el cual se transcribe a continuación y
posteriormente se establecen las estrategias para resolver el ejemplo.
3.2 SOLUCIÓN DE CASOS DE ESTUDIO DE DISEÑO, ANÁLISIS Y OPTIMIZACIÓN D

3.2.1. Desing and control of the cumene process.

"... The chemistry of the cumene process features the desired reaction of benzene with
propylene to form cumene and the undesirable reaction of cumene with propylene to form
p-diisopropyl benzene (PDIB). Both reactions are irreversible. Since the second has a higher
activation energy than the …rst, low reactor temperatures improve selectivity of cumene.
However, low reactor temperatures result in low conversion of propylene for a given reactor
size or require a large reactor for a given conversion. In addition, selectivity can be im-
proved by using an excess of benzene to keep cumene and propylene concentrations low,
but this increases separation costs. Therefore, the process provides an interesting example
of plantwide economic design optimization in which there are many classical engineering
trade-o¤s: reactor size versus temperature, selectivity versus recycle ‡owrate, and reactor
size versus recycle ‡owrate. Design optimization variables a¤ect both energy costs and cap-
ital investment. They also a¤ect the amount of reactants required to produce a speci…ed
amount of cumene product. The economic e¤ect of reactant consumption is very large, an
order of magnitude greater than the impact of energy or capital.
The process is presented in the design book by Turton et al.1 and consists of a cooled
tubular reactor and two distillation columns. The liquid fresh feeds and the benzene recycle
stream are vaporized, preheated, and fed into the vapor-phase reactor, which is cooled by
generating steam. Reactor e- uent is cooled and fed to the …rst column, which produces a
distillate stream of benzene that is recycled back to the reactor. The second column separates
the desired cumene product from the undesired p-diisopropyl benzene.

Introduction.

The cumene process illustrates some interesting design optimization features and princi-
ples. The two dominant design optimization variables are reactor size and benzene recycle.
Increasing either one reduces the amount of undesirable by-product that is produced, but
increasing either of these variables raises capital and/or energy costs.
Of more importance is the e¤ect of these design optimization variables on conversion and
selectivity, which translates into changes in the amounts of the reactants required to produce
a …xed amount of cumene. For example, if a bigger reactor is used, reactor temperature can
be lower, which improves selectivity while achieving the same conversion and also reduces
the amount of raw materials. If the reactor is inexpensive (cheap catalyst), the optimum
design is a large reactor. However, if the reactor is expensive, the optimum design is a small
reactor.
As Douglas pointed out (Douglas Doctrine) two decades ago, the costs of raw materials
and products are usually much larger than the costs of energy or capital in a typical chemical
process. Therefore, the process must be designed (investing capital and paying for energy)
so as to not waste feedstocks or lose products (particularly in the form of undesirable prod-
ucts). Process economics dictate that the conversion of reactant must be quite high. These
principles are clearly illustrated in the cumene process in several ways.
146 CAPíTULO 3 SIMULACIÓN COMERCIAL.

1. The fresh propylene feed stream contains some propane impurity, which is inert in
the reactor and must have a place to get out of the process. Since the separation of
propylene and propane is di¢ cult, the economics strongly favor designing the reactor
for a very high conversion of propylene. The propane and any unreacted propylene are
‡ashed o¤ and burned, so they only have fuel value. High propylene conversion can be
achieved by either running at high temperatures or increasing the size of the reactor.
The former increases the production of undesirable by-product. The latter increases
capital cost.
2. The undesirable by-product is also burned, so it only has value as fuel. Since it takes
reactants to produce this product instead of producing cumene, there is a strong in-
centive to keep its production rate small.
The economic analysis presented in the chapter will show that raw material fresh feed
costs are an order of magnitude larger than energy or annual capital costs.

Process studied.
The cumene process is described in Turton et al. in a fair amount of detail. The basic
‡owsheet is presented, and the kinetics of the reactions are given. This latter information
is quite useful since it is often di¢ cult to …nd kinetic information in the literature. The
authors provide a good description of the basic units in the process with equipment sizes and
operating conditions. Their intent is not to determine the optimum design but to illustrate a
typical multiunit chemical process that features a reaction section coupled with a separation
section and connected by a recycle stream.

Reaction Kinetics The production of cumene (isopropyl benzene) involves the reaction
of benzene with propylene in a high-temperature, high-pressure, gas-phase reactor.

C6 H6 + C3 H6 ! C9 H12 (3.15)
There is also a sequential reaction of cumene and propylene to form PDIB.

C9 H12 + C3 H6 ! C12 H18 (3.16)

Table 3.21 gives the reaction kinetics provided by Turton et al.1 All reaction rates have
units of kmol s 1 m 3 . Concentration units are molarity. The reactions occur in the vapor
phase in the presence of a solid catalyst (assumed to have 0.5 void fraction and a 2000 kg/m3
solid density). The reactor is run at high pressure (25 bar) since the moles of reactants are
more than the moles of product (Le Châtelier’s principle).
Notice that the activation energy of the undesirable reaction is larger than that of the
desirable reaction. Therefore low reactor temperatures improve selectivity. In addition, selec-
tivity is improved by keeping the concentration of cumene and propylene low in the reactor.
This can be achieved by using a large excess of benzene, but the excess must be recovered
and recycled.
3.2 SOLUCIÓN DE CASOS DE ESTUDIO DE DISEÑO, ANÁLISIS Y OPTIMIZACIÓN D

Tabla 3.21: Parámetros cinéticos para las reacciones de la obtención del cumeno.
Diseño <1 <2
k 2;8 107 2;32 109
E (kJ/kmol) 104; 174 146,742
Concentration terms(kmol/m3 ) CP CB CC CP

Phase Equilibrium Assuming that a separation between propylene and propane is uneco-
nomical, the only separations required are achieved fairly easily by distillation. The normal
boiling points of benzene, cumene, and PDIB are 80.1, 152.4, and 209.8 8 C, respective-
ly. The NRTL physical property package is used in the Aspen simulations used in this
chapter. Figure 3.109 gives the Txy diagrams for the benzene/cumene system and for the
cumene/PDIB system at atmospheric pressure. The curves are fairly fat, so the required
number of trays and the re‡ux ratios (RR) in the distillation columns are fairly low.

a) b)

Figura 3.109: Diagramas de equilibrio a) sistema benzeno/cumeno y b) sistema cu-

meno/PIDB.

Flowsheet Figure 3.110 shows the ‡owsheet of the process with the equipment sizes and
conditions used by Turton et al. The fresh feed streams of benzene and mixed C3 (propylene
and propane) enter the process as liquids. The fresh feed ‡owrate of the C3 feed is set at
110 kmol/h. The composition of this feed is 95 mol % propylene and 5 mol % propane. Since
propane does not react, the (110)(0.05) 41 5.5 kmol/h of propane must leave the process
somewhere. It is vented in a gas stream from the ‡ash drum.
The fresh feed of benzene is 104.2 kmol/h. Note that if there were only the cumene re-
action [Eq. (3.15)] and if all the propylene were converted and none were lost, the fresh
feed ‡owrate of the pure benzene would have to satisfy the stoichiometry of the reac-
tion:(110)(0.95) 14 94.5 kmol/h. But some of the propylene is consumed in the second reaction
148 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.110: Diagrama del proceso para la obtención del cumeno, Luyben.

and some is lost. The result is a benzene fresh feed ‡owrate of 104.2 kmol/h, which shows
that over 10 kmol/h of benzene is being undesirably consumed or lost. The production rate
of cumene is 92.86 kmol/h in the Turton et al. design.
The liquid fresh feeds are combined with a liquid recycle stream (benzene recycle D1, the
distillate from the …rst column C1) and fed into a vaporizer. The total benzene fed to the
reactor (sum of fresh plus recycle) is 207 kmol/h. The saturated gas leaves the vaporizer at
209 C and 25 bar. It is preheated in two heat exchangers. The …rst recovers heat from the
hot reactor e- uent at 427 C. The second adds additional heat to bring the reactor inlet
temperature up to 360 C.
Reaction Section. The reactor is a cooled tubular reactor that generates high-pressure
steam from the exothermic reactions. There are 342 tubes, 0.0763 m in diameter and 6 m in
length. They are …lled with a solid catalyst with a void fraction of 0.5 and a solid density of
2000 kg/m3. The temperature on the steam side of the reactor is 360 8C. An overall heat
transfer coe¢ cient of 0.065 kW m 2 K 1 is used.We assume that the reactor inlet temperature
is set to the same value as the coolant (steam) temperature in the reactor.
Reactor e- uent leaves at 427 C, is cooled to 322 C in the feed-e- uent heat exchanger
(FEHE) and sent to a condenser in which it is cooled to 90 C using cooling water. The
twophase stream from the heat exchanger is fed to a ‡ash tank. The gas from the tank is
used as fuel. The liquid is fed into the …rst distillation column C1.
It is important to note that this gas stream in the Turton et al. design is very large
(21.77 kmol/h), and it contains a large concentration of unreacted propylene (31.9 mol %
3.2 SOLUCIÓN DE CASOS DE ESTUDIO DE DISEÑO, ANÁLISIS Y OPTIMIZACIÓN D

propylene). Clearly this operation would be quite uneconomical because of the excessively
large losses of the raw material propylene.
As we will show in the optimization studies in Section 9.3, the propylene conversion-
mustbe drastically increased. There are 8.9 kmol/h of propylene in the reactor e- uent in the
Turton et al. design, which has a reactor inlet temperature and a steam temperature of 360
C. It is a simple matter to increase the propylene conversion by increasing these two tem-
peratures. However, the higher temperature will increase the production of the undesirable
product. So it may be more economical to increase the size of the reactor (use more tubes),
which would give a higher propylene conversion but at a lower temperature. This is one of
the important design trade-o¤s in the cumene process.
Benzene Recycle Column C1. This column has 15 stages and is fed on stage 6, which
is the optimum feed stage to minimize reboiler heat input. The operating pressure is 1.75
bar,which gives a re‡ux-drum temperature of 60 C, so cooling water can be used in the
condenser.The RR is small (0.44). The distillate is mostly benzene and is recycled back to
the reactor. Its composition is 95.6 mol % benzene with small amounts of propylene and
propane thatare in the liquid from the ‡ash drum.
The design speci…cation is to keep benzene from dropping out the bottom and a¤ecting
the purity of the cumene product leaving in the distillate of the downstream column. Since
the speci…ed cumene purity is 99.9 mol %, a very small benzene composition in the bottoms
(0.05 mol %) is required. The column diameter is 1.36 m, and the reboiler heat input is 1.64
Gcal/h. Both of these increase if more benzene recycle is used.
Cumene Product Column C2. This column has 20 stages and is fed on stage 12. The
operating pressure is 1 bar, which gives a re‡ux-drum temperature of 152 8C. The RR is
low (0.63). Energy consumption 1.2 Gcal/h, and the column diameter is 1.26 m.
The design speci…cation is to attain high-purity cumene in the distillate and minimize
theloss of cumene in the bottoms. The bottoms composition is set at 0.1 mol % cumene. The
distillate composition is 99.9 mol % cumene using the 0.63 RR.
Notice that the bottoms B2 ‡owrate is 1.952 kmol/h in this ‡owsheet. This represents a
loss of reactants, and the stream only has fuel value. The production of the undesired product
can be reduced by lowering reactor temperatures (which would require a larger reactor) or
increasing benzene recycle (which would require a larger distillation column C1 andwould
increase energy consumption in the vaporizer and in the C1 reboiler)... "
Para resolver esta planta de procesos se debe de seguir los siguientes pasos:

1. Identi…car los equipos presentes en el diagrama de ‡ujo, ver …gura 3.110, y encontrar
el módulo más adecuado de Aspen que represente a cada uno de ellos, para este caso se
tienen las siguientes equivalencias que se muestran en la tabla 3.22, se añade el nombre
con que se identi…cará en el simulador.

2. Identi…car las especies que se utilizarán, ya se como reactivos, inertes, servicios, etc.,
y si no está indicado en el documento, proponer el método termodinámico a utilizar.
En este caso se indica utilizar NRTL.
150 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Tabla 3.22: Equipos y nomenclatura para el proceso de obtención del cumeno.

Equipo Modelo Aspen Identi…cador
Mezclador 1 Mixer MIX1
Vaporizador Heat VAPORIZA
Intercambiador HeatX FEHE
Calentador Heat HX01
Reactor tubular Rplug REACTOR
Enfriador Heat HX01
Tanque de Flash Flash2 FLASH
Columna de destilación RadFrac C1
Mezclador 1 Mixer (MIX2) MIX2
Columna de destilación RadFrac C2

3. Identi…car las cinética involucradas, el tipo de reacción y los parámetros requeridos.

Para este caso existen dos reacciones, para las que se utilizarán el modelo de POWER-
LAW, la reaccion ocurre en fase vapor y las unidades empleadas en la ley de velocidad
es de concentración.
4. Identi…car las condiciones, presión y temperatura, de las corrientes que entran y salen
de los equipos, así como los parámetros de los equipos y las especi…caciones impuestas;
para este caso se encontro que:

a) Los mezcladores (MIX1 y MIX2) no requerien ninguna especi…cación.

b) El vaporizador (VAPORIZA) opera a 34.32 bar y la salida será un vapor saturado
c) Intercambiador (FEHE) opera sin caída de presión y la temperatura de la corriente
fría de salida es de 330 C.
d) El calentador (HX01) opera sin caída de presión y la temperatura de salida es de
360 C.
e) El reactor tubular (REACTOR) es un reactor multitubos, que opera con una
temperatura de salida de 427 C, la con…guración del reactor tiene 342 tubes, con
un díametro de 0.0763 m y una longitud de 6 m. Se utiliza un catalizador sólido
con una fracción vacío de 0.5 y una densidad de sólido de 2000 kg/m3 .
f) El enfriador (HX02) opera a 1.75 bares de presión y la temperatura de salida es
de 90 C.
g) La cámara de ‡ash (FLASH) opera a 90 C y un ‡ujo de la corriente de vapor de
21.77 kmol/hr.
3.2 SOLUCIÓN DE CASOS DE ESTUDIO DE DISEÑO, ANÁLISIS Y OPTIMIZACIÓN D

h) La primera columna de destilación (C1), que separa los reactivos sin reaccionar,
esta equipada con 15 etapas y opera con un ‡ujo de destilado 102.4 kmol/hr y
una razón de re‡ujo de 0.44, la alimentación se ubica en la etapa 8 y la columna
opera a 1.75 bar.
i) La segunda columna de destilación (C2), que separa en el destilado al producto
principal y por los fondos los productos de la segunda reacción, la columna está
equipado con 20 etapas y opera con un un ‡ujo de destilado 92.86 kmol/hr y una
razón de re‡ujo de 0.63, la alimentación se ubica en la etapa 12 y la columna
opera a 1 bar.

5. Con la información se elabora la hoja de ‡ujo que se muestra en la …gura 3.111, se

debe de observar que se añadió una bomba en la corriente que une la salida del líquido
de la cámara de ‡ash a la primera columna, debido a que la presión de la corriente de
entrada en menor que la presión de operación.

Figura 3.111: Hoja de ‡ujo del proceso para la obtención del cumeno.

La solución de este proceso se muestra en la …gura 3.112, en ella solo se muestran las
corrientes principales, la entrada (FEED1) y salida del reactor (RXN01), la corriente
vapor de la cámara de ‡ash (VAP), la corriente de reciclo (BZNRECYC) y las salidas de
la segunda columna (PRUDUCTO y B2)
En la …gura que e muestran a ontinuación se muestran las grá…cas que se obtienen de
los distintos equipos; en la …gura 3.113 se muestra el per…l de composición a lo largo de la
longitud del reactor; mientras que en las …guras 3.114 y 3.115, se muestran los per…les de
composición de las columnas C1 y C2 respectivamente.
152 CAPíTULO 3 SIMULACIÓN COMERCIAL.

Figura 3.112: Resultados del proceso de obtención de cumeno.

Figura 3.113: Per…les de concentraciones del reactor, proceso de obtención del cumeno.

3.3. Interpretación de resultados.

La ultima etapa de la simulación es la comparación de los resultados de la simulación, ya
sea con los datos reales de la planta o con datos reportados de otra simulación, para poder
evaluar el desempeño de la simulación. Es necesario ser capaz de explicar la razón de las
diferencias, que siempre existirán, entre la simulacion y los datos. Las principales razones
pudieran ser:

1. Diferente modelo termodinámico seleccionado, ya sea por consideraciones particulares o

por una mala selección; en muchos de los métodos termodinámicos existen variaciones,
ya sea por reglas de mezclado o por los valores de los parámetros binarios.
2. Modelos cinéticos diferentes, consideraciones de sobre el uso o no de los datos de
catalizadores, per…les de temperatura o condiciones especiales en los reactores.
3. Especi…caciones de diseño en las columnas de destilación, selección de parámetros o
consideraciones diferentes al simulador.
3.3 INTERPRETACIÓN DE RESULTADOS. 153

Figura 3.114: Per…les de composición en la columna C1, proceso de obtención del cumeno.

Figura 3.115: Per…les de composición en la columna C2, proceso de obtención del cumeno.

Es necesario hacer una revisón en cada una los equipos para analizar si se observan
variaciones signi…cativas en los ‡ujos, en las composiciones o condiciones de temper-
atura y presión. Los punto más críticos a revisar son los reactores, ya que al cambiar
las cantidades de productos generados modi…ca a los equipos de separación, sobre este
aspecto también es necesario veri…car la cinética empleada en eñ simulador, si no se
simpli…co demasiado el modelo o se omitieron algunas reacciones.