ESCUELA PROFESIONAL DE INGENIERIA QUIMICA

PRODUCCIÓN DE ETILENGLICOL -MÉTODO CONVENCIONAL Y

SIMULACIÓN DEL PROCESO EN ASPEN HYSYS V9.0

CÁTEDRA : ANÁLISIS Y SIMULACIÓN DE PROCESOS

CATEDRÁTICO : Ing. GUEVARA YANQUI, Pascual Víctor

ESTUDIANTES : CAMPOS CARDENAS, Julio César (IQ)

CUYUBAMBA TAPIA, María Teresa (IQ)
SEMESTRE : IX

HUANCAYO-PERÚ
2019
 MATERIA PRIMA TRANSFORMACIÓN PRODUCTOS

ENTRADA PROCESO SALIDA

ESTADO
CONTINUO ESTACIONARIO

INDEPENDIENTE DEL TIEMPO

HIDRÓLISIS DEL ÓXIDO SIN ADICIÓN DE
DE ETILENO CATALIZADOR

OBTENCIÓN DE OBTENCIÓN DE
ETILENGLICOL OTROS GLICOLES
Figura 4. Composición del producto obtenido en la hidrólisis de óxido de etileno como una función de la relación de agua a óxido de etileno
a) monoetilenglicol; b) dietilenglicol; c) Trietilenglicol; d) polietilenglicoles mayores.
Óxido de etileno
Agua Etilenglicol
 Aspecto:
Líquido aceitoso transparente.

 Solubilidad:
Miscible en agua.

 Peso Específico:
20°C/4°C

 Punto de Ebullición:
197.6 °C (388F)
 Punto de Fusión:
-13° C (9F)
 Presión de Vapor (mm Hg):
0.06 a 20°C (68°F)
 Dietilenglicol
ESTADO FISICO Líquido
COLOR Incoloro
OLOR Sin olor
PESO MOLECULAR 106.12
FORMULA MOLECULAR C4H10O3
PUNTO DE EBULLICION 245.8°C
PUNTO DE FUSION -8°C
GRAVEDAD ESPECÍFICA 1.12
SOLUBILIDAD CON AGUA Soluble (agua fría)
pH 7 neutro
VOLATILIDAD -----
VISCOSIDAD Media
SOLUBILIDAD (con solventes) Metanol, dietiléter, agua fría y caliente
 Solvente en las
Anticongelante industrias de
para carros pinturas y
plásticos

Formulaciones
Síntesis de
de tintas para
explosivos
impresoras

Agente
emoliente en el
celofán
SECTOR MANUFACTURERO

Tabla N°5. Producción de las industrias de sustancias y productos químicos

Unidad 2015 2016

Fuente: Ministerio de la Producción-Viceministro de MYPE e Industria

 ANÁLISIS DE DEMANDA
Tabla N°11. Principales empresas importadores de etilenglicol
a) Para el producto: etilenglicol
4.4.1 IMPORTACIÓN DE ETILENGLICOL
- US$: 3,614,222
- Importadores: 36
- Cantidad importada: 57, 953,040 kg
- Precio unitario: 1.257 US$/kg

4.4.2 EXPORTACIONES DE ETILENGLICOL

- US$: 0
- Exportadores: No especificado
- Cantidad importada: 0 kg

Fuente: www.veritradecorp.com
Tabla N°12. Principales Países de origen de importación de etilenglicol

Fuente: www.veritradecorp.com
 CONDICIONES INICIALES

Flujo Temp
Presi
Fracción molar eratu
Corriente ón
molar (kmol ra
(atm)
/h) (°C)

Óxido de
etileno 1 250 25 10
Agua 1 500 25 10
CONSIDERACIONES:

• Las reacciones involucradas son las que se describieron

anteriormente. ocurren en fase líquida y no existe caída
de presión.
• La longitud del reactor es de 10 m y el diámetro de 2,210
m. (Melhem, 1996)
• “A altas temperaturas, hay una significativa formación de
subproducto, mientras que a temperaturas inferiores a 40
°C la reacción no se realiza a velocidad significativa, en
consecuencia, se eligió una temperatura de 55 °C”.
• Operación isotérmica.
Tabla N° 17. Constantes de velocidad con la función de Arrhenius

Fuente: (Benz, 2015)

DESCRIPCIÓN

• Antes del ingreso de la alimentación existe una

válvula en la cual se baja la presión.
• Se requiere la pérdida de carga de presión.
• Previamente se disminuye la presión en 8 atm por
medio de una válvula; para finalmente obtener en el
destilado mayor porcentaje de agua, en los fondos el
etilenglicol y el subproducto no deseado el
dietilenglicol.
Tabla N°14. Especificaciones para el funcionamiento de la columna de destilación

FUENTE: (Ikhyun, 2016)

 1. DISEÑO DEL REACTOR
Ya que el objetivo de la planta es producir 57 953 040 kg de etilenglicol por año, transformando a
kmol por hora se obtiene:

𝑇𝑛 1000 𝑘𝑔 1 𝑘 𝑚𝑜𝑙 1 𝑎ñ𝑜 1 𝑑𝑖𝑎 89.5𝑘𝑚𝑜𝑙

𝑃𝐸𝐺𝐿𝐼𝐶 = 57 953 040 ∗ =
𝑎ñ𝑜 1𝑇𝑛 76 𝑘𝑔 355 𝑑𝑖𝑎𝑠 24 ℎ ℎ

Primero realizaremos el balance en el reactor PFR, para calcular los requerimientos de materia
prima y al mismo tiempo diseñar el reactor.
 FLUJO INICIAL DE ÓXIDO DE ETILENO

89.5 𝑘𝑚𝑜𝑙
𝐹𝐴0 = = 111.88
0.8 ℎ
 CÁLCULO DEL VOLUMEN

𝐹𝐴0 𝑋
𝑉=
−𝑟𝐴

111.88 × 0.8
𝑉= 79376ൗ
1.2𝑒 16 𝑒 − 1.987×55 × 0.8 × 0.2
V=38.30 m3
En base a la bibliografía asumiremos que el radio del reactor es la mitad de su altura, así queda:
V=2πr3=38.30m3
r=1.1 m
h=9.8 m
2. DISEÑO DE LA COLUMNA

500𝑘𝑚𝑜𝑙
𝐹=
ℎ

Para encontrar los moles de agua a la salida del reactor:

𝐹𝐴 = 𝐹 ∗ 𝑋𝐴

330.25𝑘𝑚𝑜𝑙
𝐹𝐴 =
ℎ

Para los moles de etilenglicol haremos lo propio:

𝐹𝐸 = 𝐹 ∗ 𝑋𝐸

89.5𝑘𝑚𝑜𝑙
𝐹𝐸 =
ℎ

Considerando una relación de reflujo de 0.9 (basándonos en la bibliografía) y además observando las
propiedades de los tres compuestos en juego tendremos las características que se mencionan para la
operación:
 Compuesto Óxido de Etileno Agua Etilenglicol

Punto de ebullición (°C) 10.7 96 - 97 195

•Asumiremos que el óxido de etileno por su bajo punto de ebullición no afectará en la separación por lo que consideraremos destilación binaria.
•Producto de cabeza (destilado) 0.01 % de etilenglicol
•Producto de colas 46.82 % de etilenglicol

La composición en la entrada será:

𝑋𝐸 = 0.1790

Realizando el respectivo balance en la torre:

500=D +W W=190.1 kmol/h

0.1790 ×500=0.0001 D +0.4682 W D=309.9 kmol/h

Teniendo en cuenta que:
V =L+ D

𝐷 𝐿
= 1 − = 1 − 0.9 = 0.1
𝑉 𝑉

yn=0.9 xn−1+ 0.002

Para conocer L:

L/V = 0.9 V=L+D

1
𝐿 = 309.9 ∗ = 2789.1
1
0.9 − 1
En consecuencia:

L' =2789.1+500=3289.1 kmol/h

El caudal molar por debajo de la entrada de la alimentación vendrá dado por:

V ' =L' −W

V ' =3289.1−190.1=3099

Con lo que queda:

ym=1.061 xm−1−0.035

Asumiendo por definición que la composición del vapor procedente del piso 1 es la misma que
la del destilado y la del líquido que vuelve como reflujo a la columna, es decir:
y1=xD =0.0001

Los valores en el equilibrio predichos para la presión son:

 Sistema Agua – Etilenglicol
T ºC Fracción mol Fracción mol
líquido vapor
189,23 0 0
170,51 0,05 0,4834
157,16 0,1 0,6963
147,19 0,15 0,8042
139,4 0,2 0,8655
133,09 0,25 0,9034
127,82 0,3 0,9284
123,34 0,35 0,9457
119,44 0,4 0,9582
116 0,45 0,9675
112,94 0,5 0,9745
110,19 0,55 0,98
107,68 0,6 0,9844
105,4 0,65 0,9879
103,29 0,7 0,9907
A partir de estos datos de equilibrio podemos obtener por interpolación:

x1=0.08

Con este valor vamos a la ecuación obtenida para la parte superior de la alimentación en la columna:

y2=0.9 ×0.05+ 0.0001=0.0851

Si continuamos con éste método analítico (método de Sorel y Lewis), obtenemos los siguientes valores:

X10=0.053

Como el valor obtenido en el tercer plato está por debajo de la composición de alimentación, la alimentación
deberá ser en éste plato. Ahora continuaremos con el cálculo total de platos necesarios:

y10=1.061 ×0.053−0.035=0.0212
x11=0.0021
El resultado obtenido está por debajo de la composición requerida por lo que podemos concluir que nuestra torre de
destilación tendrá 10 platos
RESULTADOS DEL MÉTODO CONVENCIONAL
MEZCLADOR

T = 50°C
P = 10 atm
F=750 Kg P =750 Kg
mol/h mol/h

MEZCLADOR

W F, OE=250 Kg mol/h = 0.33

W F, OE=250 Kg mol/h W F, AGUA=500 Kg mol/h = 0.67
W F, AGUA=500 Kg mol/h
 REACTOR

T = 55°C
P = 10 atm

V =500 Kg mol/h

F=750 Kg mol/h
REACTOR PFR

P =250 Kg mol/h
 DESTILADOR

T = 151.758 °C
T = 115.72 °C P = 17 psia
P = 15 psia
D =444.7280 Kg mol/h

W F, OE= 2.1569 Kg mol/h = 0.0005

W F, H2O= 440.41 Kg mol/h = 0.99
R= D/L= 2 W F, ETG= 2.1569 Kg mol/h = 0.0005

P =1000.1 Kg mol/h
T = 98.6212 °C
W F, OE= 0.1 Kg mol/h = 0.00001 P = 10 atm
W F, AGUA= 500.05 Kg mol/h = 0.5
W F, ETG= 499.95 Kg mol/h = 0.4999

F =555.3720 Kg mol/h

W F, H2O= 59.53 Kg mol/h = 0.1072

W F, ETG= 495.83 Kg mol/h = 0.8928
 RESULTADOS

MEZCLADOR MEZCLA PRODUCTO VENTEO DESTILADO FONDO

ESPECIE 1 (kmol/h) 2 3 4 5 6 7
(kmol/h) (kmol/h) (kmol/h) (kmol/h) (kmol/h) (kmol/h)

ÓXIDO DE 250 - 250 0.05 - 440.4141 59.5359

ETILENO

AGUA - 500 500 500.05 - - -

ETILENGLICOL - - - 499.95 - 4.3138 495.8361

DIETILENGLICOL - - - 0.05 - - -

TOTAL - - 750 1000.1 - 444.728 555.372

 SIMULACION

•SELECCIÓN DE LA OPCIÓN NEW

•SELECCIÓN DE COMPONENTES
Clic en el botón Add y elegir componentes
 •SELECCIÓN DEL PAQUETE TERMODINÁMICO
Clic en Add y seleccionar el paquete NRTL

•REACCIONES
Clic en Add, agregar reacciones y el tipo que corresponde a cinética.
 •SIMULACIÓN
Clic en simulación

•AGREGAR CORRIENTES
Mediante la opción Palete, seleccionar:
 •CONDICIONES DE ENTRADA DE MATERIAS PRIMAS
Doble clic en la flecha azul y rellenar:

•AGREGAR MIXER
En la paleta seleccionamos:
 •CONDICIONES INICIALES DEL MEZCLADOR
Doble clic en el mixer y completar:

•AGREGAR EL REACTOR PFR

En la paleta seleccionar:
 •CONDICIONES DE ENTRADA AL REACTOR
Doble clic al reactor y completar:

•VÁLVULA REDUCTORA DE PRESIÓN

En la paleta, agregar:
 •CONDICIONES DE LA VÁLVULA
Doble clic en el equipo y completar la presión:

•COLUMNA DE DESTILACIÓN
En la paleta ubicar en la sección Columns:
 •CONDICIONES PARA LA COLUMNA
Doble clic y completar:

•Para correr la simulación, seleccionamos RUN

 •RESULTADOS OBTENIDOS
Para conocer los flujos de masa y molar del proceso, se hace
clic derecho en cualquier parte del entorno de la simulación,
aparece una ventana en donde se escoge la opción Add
Workbook Table como se muestra en la siguiente figura.

•TABLA DE CORRIENTES DEL

PROCESO
 VIII.COMPARACIÓN Y DISCUSIÓN DE RESULTADOS

Oxido de Etileno Agua Mezcla Vapor Productos Mezcla a destilación Venteo Destilado Etilenglicol
Vapour Fraction 0 0 0 1 0 0 1 0 0
Temperature °C 50 50 50 98.6212 98.6212 98.63 97.5 97.5 140.3
Pressure kPa 1013 1013 1013 1003 1003 810.6 100 100 105
Molar Flow kgmol/h 500
250 1000
500 1500 0 1000.1 1000.1 0 444.728 555.372
Mass Flow kg/h 2.2*10^4 1.8*10^4 4*10^4 0 4*10^4 4*10^4 1.91*10^-4 8054 3*10^4
Liquid Volume Flow m3/h 24.5 18.05 43 0 35 35 1.2*10^7 8.05 29
Heat Flow kJ/h -39000000 -300000000 -310000000 0 -370000000 -370000000 -2.6 -185000000 -240000000

HYSYS EN ESTADO ESTACIONARIO

GRACIAS