PRÁCTICA 1.

SIMULACIÓN DE UNA PLANTA DE PRODUCCIÓN DE AMONÍACO

Roberto Fabio Díaz Bermúdez, María Alejandra González Bedoya, Raul Esteban
Mina Riascos, Aura Cristina Taborda Acosta, Miller Damian Zarate Riveros.
Estudiantes de noveno semestre de Ingeniería Química.
Universidad Pontificia Bolivariana, Cq. 1 #70-01, of. 11-259, Medellín, Colombia.
Resumen
Se simuló el proceso de producción de amoníaco en Aspen Hysys con el objetivo de comparar los resultados
con los obtenidos por otro autor en el simulador Chemcad, y, además, fortalecer los conocimientos básicos
en la simulación de equipos empleados en un proceso de ingeniería como reactores, separadores,
intercambiadores de calor, entre otros. Se tomaron de referencia de la simulación en Chemcad los valores
de entrada del proceso, la temperatura a la salida del condensador, la relación de recirculación y la caída de
presión en algunos equipos. En Aspen hysys se realizó la simulación de la reacción en un reactor de
equilibrio, a diferencia de la simulación en Chemcad donde se utilizaron varios reactores empacados en
serie con intercambiadores de calor. Debido a esto, la conversión obtenida en Hysys (36.83%) es mayor a
la obtenida en Chemcad (35.31%), porque el reactor de equilibrio simula la máxima conversión posible a
las condiciones dadas. También se obtuvo una producción similar para ambas simulaciones con una
diferencia de 16tn/día más en Hysys, la cual se debe principalmente a la mayor conversión obtenida y a la
tolerancia de este simulador al integrar un reciclo. El modelo termodinámico y las reglas de mezclado
utilizadas en hysys (Peng – Robinson) también hacen que existan diferencias en las entalpías de las
corrientes y por ende en los calores generados, ya que en Chemcad se utiliza un modelo termodinámico
diferente el cual no está especificado.

Palabras clave: Simulación, procesos, conversión, reactor, Aspen Hysys.

1. Introducción
El amoníaco anhidro es un gas incoloro, de olor irritante y tóxico, es muy soluble en agua, se hidrata
formando NH4OH que se ioniza, generando soluciones de fuerte carácter básico. Por su importancia como
nutriente es el segundo producto químico que más se produce industrialmente a escala mundial después del
ácido sulfúrico. Es la materia prima para la fabricación del ácido nítrico, del nitrato amónico y otros nitratos
inorgánicos, así como de la urea, todos ellos de empleo masivo como fertilizantes (Rodríguez, et al, 2007).
La fabricación de amoníaco sintético tuvo su éxito inicial en el proceso de Haber, en el cual una mezcla de
nitrógeno e hidrógeno se hace circular sobre un catalizador a una temperatura moderadamente elevada,
realizándose la transformación parcial de los gases elementales en amoníaco. Las diversas modificaciones
introducidas para realizar la síntesis del amoníaco a partir de sus elementos han tenido tanto éxito que puede
afirmarse que en la actualidad este proceso es más importante que todos los demás métodos sintéticos
reunidos (Albaladejo, 2013). Las materias primas que se emplean por lo general son aire, agua,
hidrocarburos y energía.
Albaladejo realizó un estudio del proceso de la producción de amoníaco para modificar una planta existente
con el fin de aumentar la producción de 1000 Tn/día a 1500 Tn/día, con una concentración del 93% de
amoníaco (Albaladejo, 2013). Realizó una simulación en Chemcad para estimar los cambios fisicos,
químicos y termodinámicos que se deben tener para llegar a la capacidad deseada.
La producción consta esencialmente de seis fases:
a. Obtención de los gases reactantes.
b. Purificación.
c. Compresión.
d. Reacción catalítica.
e. Recuperación del amoníaco formado.
f. Recirculación.
El proceso propuesto se muestra en la figura 1 y los resultados obtenidos por Albaladejo se presentan en la
Tabla 1.

(7 en hysys)
(9 en hysys)

(8 en hysys)

(6 en hysys)

(5 en hysys)

(10 en hysys)

Figura 1. Proceso de producción de amoníaco simulado por (Albaladejo, 2013).

Corriente en Equivalente Fracción Temperatura Presión Flujo Molar Flujo másico

Chemcad en Hysys vapor (°C) (kPa) (kmol/h) (kg/h)
1 1 1 40 2000 8,247.6 71460
2 2 1 54 10000 8,247.6 71460
3 3 1 105.62 10000 23,778.0 219564
4 4 1 377 10000 20,077.2 219564
7 5 0.815 -28 10000 20,077.2 219564
8 6 1 -29.98 9500 16,351.2 155934
9 7 1 -29.98 9500 15,534.0 148140
10 8 1 -24.42 10000 15,534.0 148140
11 9 1 -29.98 9500 817.6 7796
12 10 0 -29.98 9500 3,724.9 636264
Tabla 1. Resultados obtenidos en la simulación en Chemcad. (Albaladejo, 2013)

La reacción que se lleva a cabo para producir el amoníaco es la siguiente:

𝑁2 + 3𝐻2 ↔ 2𝑁𝐻3
Esta reacción al ser exotérmica emplea una velocidad lenta controlando bajas temperaturas debido a los
requerimientos del sistema de reacción que tiene una alta energía de activación (Barreiro et al, 2015). Lo
anterior se puede comprobar desde la curva de equilibrio que cual muestra que la reacción en un solo reactor
se debería se realizar a una baja temperatura para lograr una alta conversión, sin embargo, la velocidad de
reacción se comportaría a una velocidad demasiado lenta para las abastecer las demandas de producción y
de igual manera se necesitaría un reactor de volumen muy grande. Por este motivo es preferible usar un
sistema de reactores adiabáticos con enfriamiento entre etapas para lograr una alta conversión. El proceso
realizado por Albaladejo esquematizado en la curva de equilibrio de la reacción se presenta en la Figura 2.

Figura 2. Curva de equilibrio de la reacción para producir amoníaco (Albaladejo, 2013).

Después de analizar el proceso de simulación de Chemad se realizó una simulación en Aspen Hysys del
mismo proceso con las mismas especificaciones correspondientes en los corrientes y los equipos
implementados para la producción de amoníaco usando como modelo termodinámico Peng – Robinson, el
diagrama de éste se presenta en la figura 3. La simulación se realizó con el fin de comparar los resultados
obtenidos en Chemcad y Hysys, además de aprender a simular equipos de operaciones en este tipo de
software. En este caso se utilizó un reactor de equilibrio en lugar de varios PBR debido a que, como se
puede observar en la Figura 2, en la simulación realizada por Albaladejo se obtiene una conversión cercana
a la del equilibrio, por lo tanto, con esa base se realizó la suposición de simular más fácilmente en Hysys
como un reactor de equilibrio.

Figura 3. Simulación en Hysys de una planta de producción de amoníaco.

La simulación en Hysys consistió en tomar inicialmente los datos de entrada en la corriente de reactor y así
ir abarcando cada uno de las operaciones en forma de cebolla dados por Albaladejo en (Tabla 1) para la
corriente de vapor formada por gases de nitrógeno, hidrógeno, argón y metano, para posteriormente simular
los siguientes equipos del proceso. También se fijaron igual a la simulación en Chemcad los valores de la
temperatura a la salida del condensador (corriente 4), la relación de recirculación en el mezclador (95% de
recirculación) y las caídas de presión en los compresores.
2. Análisis de resultados
Los resultados obtenidos en la simulación en Aspen Hysys se presentan en la Tabla 2.

Corriente Fracción Temperatura Presión Flujo Molar Flujo másico

vapor (°C) (kPa) (kmol/h) (Kg/h)
1 1 40.00 2000 8247.600 71451.19
2 1 285.23 10000 8247.600 71451.19
3 1 90.35 10000 21860.063 203690.61
4 1 386.16 10000 18114.808 203688.74
5 0.791 -28.15 10000 18114.808 203688.74
6 1 -28.44 9500 14342.115 139394.21
7 1 -28.44 9500 13625.010 132424.50
8 1 -23.57 10000 13612.463 132239.42
9 1 -28.44 9500 717.106 6969.71
10 0 -28.44 9500 3772.693 64294.53
Tabla 2. Resultados obtenidos en la simulación en Aspen Hysys.

Con el propósito de comparar los resultados, se grafican las temperaturas y los flujos másicos obtenidos en
cada corriente para la simulación con Chemcad y con Aspen Hysys. Los resultados se presentan en las
figuras 4 y 5.

Temperatura Flujo másico

500 3.E+05
Flujo másico (kg/h)

400 2.E+05
Temperatura (°C)

300 2.E+05
200 1.E+05

100 5.E+04

0 0.E+00
0 1 2 3 4 5 6 7 8 9 10 11 1 2 3 4 5 6 7 8 9 10
-100
Corriente Corriente

Chemcad Hysys Chemcad Hysys

Figura 4. Temperaturas obtenidas en ambos Figura 5. Flujos másicos obtenidos en ambos

simuladores. simuladores.

En el caso de la temperatura, se puede observar que en general las temperaturas de ambas simulaciones se
aproximan. En la única corriente en donde se presenta una diferencia significativa es en la 2. Sin embargo,
el valor reportado por Albaladejo parece ser incorrecto porque la temperatura a la salida del mezclador
(corriente 3) es mucho mayor a las dos temperaturas de entrada (corriente 2 y 8), lo cual no tiene sentido
ya que debería ser un valor entre estas temperaturas. Por otra parte, en el caso de los flujos másicos, no se
 presenta una diferencia significativa en ninguna de las corrientes. La diferencia en los flujos másicos es
debida al método empleado por Hysys para añadir una corriente de recirculación (corrientes 8 y 8’ de la
figura 3) al proceso, ya que debido a la tolerancia del simulador puede variar en décimas el flujo antes y
después del reciclo (RCY-1).
La Tabla 3 presenta la producción que se obtendría en el proceso de acuerdo con cada simulación. En ambos
casos se cumple con el objetivo de aumentar la producción a 1500 tn/día, sin embargo, en la simulación
realizada en Hysys se obtiene una producción mayor, con 16 toneladas más por día. Esta diferencia puede
ser debida principalmente al tipo de reactor empleado, ya que en esta simulación se utilizó un reactor de
equilibrio el cual puede alcanzar la máxima conversión posible en la reacción, resultando en una mayor
cantidad de producto. Sin embargo, se debe tener en cuenta que la producción obtenida en hysys se aleja
un poco de la realidad porque, según la curva de equilibrio de la reacción, es complicado alcanzar la
conversión de equilibrio. La mejor opción en la realidad es poner varios reactores PBR con
intercambiadores de calor como lo propone Albaladejo.

Producción tn/h tn/día

Hysys 64.2945 1543.0687
Chemcad 63.6264 1527.0336
Diferencia 0.6681 16.0351
Tabla 3. Producción obtenida en ambas simulaciones.

Para determinar la conversión y analizar la reacción y el producto obtenido, las composiciones obtenidas
en cada corriente se presentan en la Tabla 4. Los valores de las composiciones para estas mismas corrientes
obtenidas en Chemcad se muestran en la Tabla 5.

Corriente 1 2 3 4 5 6 7 8 9 10
Nitrógeno 0.247 0.247 0.2326 0.1773 0.1773 0.2237 0.2237 0.2238 0.2237 0.0007
Hidrógeno 0.742 0.742 0.7003 0.535 0.535 0.6748 0.6748 0.6751 0.6748 0.0035
Argón 0.003 0.003 0.0126 0.0152 0.0152 0.0184 0.0184 0.0184 0.0184 0.0027
Metano 0.008 0.008 0.0402 0.0485 0.0485 0.06 0.06 0.0597 0.06 0.0048
Amoníaco 0 0 0.0144 0.2241 0.2241 0.0231 0.0231 0.0231 0.0231 0.9883
Tabla 4. Fracciones molares obtenidas en Aspen Hysys.

Corriente Chemcad 1 2 3 4 7 8 9 10 11 12
Equivalente en hysys 1 2 3 4 5 6 7 8 9 10
Nitrógeno 0.247 0.247 0.228 0.178 0.178 0.218 0.218 0.218 0.218 0.002
Hidrógeno 0.742 0.742 0.7027 0.5556 0.5556 0.6818 0.6818 0.6818 0.6818 0.0019
Argón 0.003 0.003 0.0126 0.015 0.015 0.0178 0.0178 0.0178 0.0178 0.0027
Metano 0.008 0.008 0.0426 0.0504 0.0504 0.061 0.061 0.061 0.061 0.004
Amoníaco 0 0 0.0139 0.2008 0.2008 0.0213 0.0213 0.0213 0.0213 0.9894
Tabla 5. Fracciones molares obtenidas en Chemcad (Albaladejo, 2013).

A partir de los resultados se determinaron las conversiones obtenidas en los reactores de cada simulación,
según los flujos molares en referencia al nitrógeno en las corrientes 3 y 4 que son las correspondientes a la
entrada y salida del reactor. En el caso de Chemcad se obtuvo una conversión del 35.31% con una serie de
reactores PBR y en Hysys se obtuvo una conversión del 36.83% con un reactor de equilibrio.
La conversión otorgada por el hysys presenta una desviación de 4.3% respecto a la de Chemcad. Tiene
sentido que esta conversión sea mayor que la conversión dada por el Chemcad puesto que el reactor
utilizado en la primera es un reactor de equilibrio que ayuda a que el sistema converja rápidamente y a que
el equipo alcance la conversión de equilibrio. Sin embargo, la diferencia entre las conversiones es casi
despreciable y es debida a la serie de intercambiadores puestos justo antes de cada reactor en la simulación
en Chemcad, porque como se puede observar en la figura 2, el complejo de intercambiadores y reactores se
va aproximando hasta casi llegar a la línea de equilibrio, es decir, hasta casi llegar a la conversión de
equilibrio.
Por otra parte, en la Tabla 5 se puede observar que el producto final (corriente 10) tiene una alta
concentración de amoníaco en ambas simulaciones, de 98.83% en hysys y 98.94% en chemcad. Lo que
indica una separación efectiva de los productos de la reacción. Sin embargo, el amoníaco utilizado
industrialmente tiene una pureza mayor al 99.9% (Air products and Chemicals, 2014), por lo tanto, sería
necesario un tratamiento posterior para aumentar la pureza del producto. La diferencia entre ambas
concentraciones se debe principalmente al modelo termodinámico utilizado ya que de éste dependen las
constantes de equilibrio en la separación, sin embargo, la diferencia es pequeña.
Se determinó también el calor retirado en el condensador en la simulación hecha por Albaladejo con el fin
de comparar el resultado. El calor está dado por:
Q=UA(Tsalida-Tentrada)
En el artículo se da un coeficiente global de transferencia de calor (U) y un área de transferencia de calor
(A) de:
U=921.6 W/m2K
A=1376.26 m2
Con una temperatura de entrada y salida de 650K y 308.8K respectivamente, el calor retirado por el
condensador es de:
Q=15.58x108 kJ/h
La magnitud de este mismo calor proporcionado por la simulación en Hysys es de:
Q=3.317x108 kJ/h
Lo cual, comparando estos dos valores, nos da una diferencia que se puede deber a que en el proceso de
Albaladejo se incluyen algunos intercambiadores, como se puede observar en la figura 2, que les ceden
calor a las corrientes de entrada del reactor, con el objetivo de adecuar las temperaturas, mientras que en la
simulación en hysys no hay intercambiadores previos al reactor, como se puede observar en la figura 3. Por
lo tanto, el calor que se tiene que retirar en el artículo de Albaladejo es mayor al que se tiene que retirar en
la simulación en Hysys.
La tabla 6 muestra el calor correspondiente para el compresor K-100 (Chemcad) y K-101 (Hysys).

Simulación Q (kJ/h)
Chemcad 78616123.2
Hysys 59650000
Tabla 6. Calor compresor en ambas simulaciones.

En ambas simulaciones se tienen las mismas condiciones para las corrientes que componen el equipo, a
pesar de esto, el calor requerido por el compresor en la simulación de Chemcad es mayor. Esta diferencia
es debida al cálculo de las respectivas entalpías de las corrientes en cada simulación. En la simulación
realizada en Hysys se utilizó el modelo termodinámico con regla de mezclado de Peng-Robinson mientras
la realizada en Chemcad no define cual modelo termodinámico fue el implementado para realizar los
cálculos correspondientes.
 3. Conclusiones
La simulación del proceso de Haber para la producción de amoníaco en Aspen Hysys a partir de una
corriente de entrada de composición y condiciones conocidas. Se simularon equipos de reacción como un
reactor de equilibrio, de separación, de compresión y de mezclado.
Se compararon los resultados de la simulación con los obtenidos en otra simulación realizada en Chemcad.
Se obtuvo una conversión similar para ambas simulaciones, 35,31% en Chemcad y 36.83% en Hysys. La
diferencia entre estas conversiones es debida a que en Hysys se utilizó un reactor de equilibrio el cual simula
la máxima conversión posible, mientras que en Chemcad fue realizada con una serie de PBR con
intercambiadores de calor el cual se acerca a la conversión de equilibrio, pero es inferior. Sin embargo, esta
última configuración tiene mayor sentido para una planta real.
También se obtuvieron flujos de producción similares que cumplen con el requerido inicialmente (1500
Tn/día). El flujo del producto es de 1543.0687 Tn/día y 1527.0336 Tn/día para Hysys y Chemcad
respectivamente. La diferencia es debida nuevamente al tipo de reactor utilizado, porque la conversión en
el reactor de equilibrio es mayor al sistema simulado en Chemcad. Por otra parte, la pureza de ambos
productos es aproximadamente 98%, un valor inferior al usado industrialmente.
El calor retirado en la simulación realizada en Chemcad es de15.58x108 kJ/h, mientras que este mismo calor
en la simulación hecha en Hysys es de 3.317x108 kJ/h, esta gran diferencia puede deberse a la serie de
intercambiadores de calor que posee la simulación en Chemcad y de la cual carece la simulación en el
hysys. Ésto incide directamente, puesto que, en el Chemcad, a las corrientes de entrada se le suministra
energía y a esto se le suma la energía proporcionada por la reacción exotérmica, por lo tanto, tiene sentido
que en el proceso de condensación se le tenga que retirar más energía.
4. Referencias
Rodríguez, M; Díaz, I; González, E.J y González, M. (2007). Amoníaco y sus principales derivados.
13/02/2019, de Chemical Technology Lab Sitio web:
http://www.diquima.upm.es/old_diquima/docencia/tqindustrial/docs/cap3_amoníaco.pdf.

Albaladejo Carril, A. M. (2013). Revamping de una planta de producción de amoníaco. Universidad

Politécnica de Cartagena.
Barreiro, J. A; Coveña, A. H; Ponce, M. I; Solórzano, R. F; Zambrano, K. E. (2015). Proyecto de simulación
de procesos químicos. Universidad Técnica de Manabí.
Air Products and Chemicals. (2014). Amoníaco, datos técnicos. Formas de suministro. Especificaciones de
producto. 18/02/10, de http://www.refrigerantes.mobi/PDF/Amoniaco-brochure_ES.pdf.