UNIVERSIDAD NACIONAL AUTÓNOMA DE

HONDURAS
Facultad de Ingeniería
Departamento de Ingeniería Química Industrial

Laboratorio de Operaciones Unitarias

IQ-513

“Simulación 5: Destilador Flash”

Catedrático: Ing. Andrés Abelino Cruz

Presentado por:

Fabiola Alejandra Barahona 20161005197

Julio Raúl Calix 20122004792
Ethan Fabricio Marín 20091010217
Resumen.
El siguiente informe presenta los resultados obtenidos empleando un software para
un proceso de simulación de mezcla a través de un destilador flash, permitiendo así
analizar el comportamiento del destilador y las propiedades de la mezcla que es
introducida en este.

Introducción.
La destilación es una operación unitaria o un proceso de separación física de
componentes más volátiles separados de una mezcla por vaporización parcial de la
mezcla líquida seguida de condensación. La distribución de los componentes en las
dos fases se rige por la relación de equilibrio líquido vapor. La forma más elemental
de destilación simple es la destilación instantánea, en la que una mezcla de
alimentación se vaporiza parcialmente para producir vapor que está enriquecido con
componentes más volátiles y un líquido que carece de componentes volátiles. Las
corrientes de vapor y líquido que salen del tambor están en equilibrio. En
termodinámica e ingeniería química, el equilibrio vapor-líquido (VLE) describe la
distribución de una especie química entre la fase de vapor y una fase líquida.
La concentración de un vapor en contacto con su líquido, especialmente en
equilibrio, se expresa a menudo en términos de presión de vapor, que será una
presión parcial (una parte de la presión total del gas) si hay otros gases presentes
con el vapor. La presión de vapor de equilibrio de un líquido depende en general en
gran medida de la temperatura. En el equilibrio vapor-líquido, un líquido con
componentes individuales en ciertas concentraciones tendrá un equilibrio de vapor
en el cual las concentraciones o presiones parciales de los componentes de vapor
tienen ciertos valores dependiendo de todas las concentraciones del componente
líquido y de la temperatura. La concentración de equilibrio de cada componente en
la fase líquida es a menudo diferente de su concentración (o presión de vapor) en
la fase de vapor, pero existe una relación. Los datos de concentración de VLE pueden
ser determinados experimentalmente, o computados o aproximados con la ayuda
de teorías tales como ley de Raoult, ley de Dalton, ley de Henry o modelos más
complicados.
En las mezclas que contienen dos o más componentes, las concentraciones de cada
componente se expresan a menudo como fracciones molares. Cuando se alcanza
una temperatura tal que la suma de las presiones de vapor de equilibrio de los
componentes líquidos se vuelve igual a la presión total del sistema (de lo contrario
es menor), entonces las burbujas de vapor generadas por el líquido comienzan a
desplazar al gas que mantenía la presión total, y se dice que la mezcla “Hierve”. Esta
temperatura se denomina "punto de burbuja" de la mezcla líquida a la presión dada.
Y tiene una composición de vapor asociada. Lo mismo ocurre cuando se reduce la
temperatura o aumenta la presión, la mezcla empieza a “Condensar” y este punto
se conoce como “punto de roció”, la cual tiene una composición liquida asociada. En
medio de estos dos puntos ocurren mezclas de diversas composiciones en ambas
fases en equilibrio a medida se cambia la presión o temperatura.

Ley de Raoult.
La ley de Raoult es uno de los métodos para describir el equilibrio vapor líquido pro
de solución ideal. Para una solución ideal, la presión parcial de un componente es
igual al producto de su fracción molar en la fase líquida y su presión de vapor en las
condiciones prevalecientes. Se puede escribir como:
𝑃𝑦𝑖 = 𝑥𝑖 𝑃𝑣𝑖
Donde x, es la fracción molar en fase liquida y la fracción molar en fase vapor
representada por y.
Ley de Dalton.

𝑃 = ∑ 𝑃𝑖 = ∑ 𝑦𝑖 𝑃 = ∑ 𝑥𝑖 𝑃𝑣𝑖

Objetivos.
• Analizar el proceso de destilación instantánea.
• Desarrollar habilidades para simular un proceso de destilación.
• Determinar propiedades de mezclas y razonar su comportamiento
Procedimiento.
1. Analizaremos una destilación multicomponente de una mezcla de Benceno –
N-Heptano a distintas proporciones usando la herramienta de simulación
COCO/ChemSep.
2. La mezcla se alimentará a una temperatura 80°C y a 1 atm de presión, con
distintas proporciones.
3. La presión permanecerá constante a 1 atm en la columna de destilación
4. Inicialmente se usarán 20 etapas y la alimentación entrara en la etapa
número 10.
5. Las especificaciones de columna que se usaran será un radio de reflujo de 2
para la parte superior y una fracción combinada de recuperación para ambos
flujos de 0.5 para la parte inferior de la columna.
6. Comprobar el punto de burbuja para la mezcla y el número de etapas
requerido para obtener un destilado con un alto porcentaje de recuperación.

Asignación
1. Una mezcla de 45mol% de benceno y 55mol% de tolueno se está destilando de
forma flash a una velocidad de 15 kmol/hrs a 1 atm de presión total. Calcular las
cantidades y la composición de los productos superiores e inferiores si el
producto líquido no debe contener más de 35 mol% de benceno. La volatilidad
relativa del benceno en esta mezcla es de 2.6. El componente más volátil es el
benceno, por tanto, es el componente 1 y el tolueno el componente 2.
Datos
• Z1=0.45
• Z2=0.55
• F=15 kmol/hora
• Volatilidad relativa=2.6
• P=1 atm
Balance Global
𝐹 =𝐿+𝑉
15 = 𝑉 + 𝐿
Balance por componente (Benceno)
𝑍1 𝐹 = 𝑥1 𝐿 + 𝑦1 𝑉
0.45 ∗ 15 = 𝑥1 𝐿 + 𝑦1 𝑉
6.75 = 𝑥1 𝐿 + 𝑦1 𝑉
Balance de Operación
𝑍1 = 𝑥1 𝐿
0.45 ∗ 15 ∗ 0.35 = 𝑥1 𝐿
2.362 = 𝑥1 𝐿
𝐿 = 2.362/𝑥1
Sustituyendo L en Balance por componente
2.362
6.75 = 𝑥1 + 𝑦1 𝑉
𝑥1
𝑦1 𝑉 = 6.75 − 2.3625
4.3875
𝑉=
𝑦1
Volatilidad Relativa
𝑦1
𝑥
𝛼 = 𝑦1
2
𝑥2
Sustituyendo
𝑥1 + 𝑥2 = 1
𝑥2 = 1 − 𝑥1
𝑦1 + 𝑦2 = 1
𝑦2 = 1 − 𝑦1

𝑦1
𝑥1
𝛼=
1 − 𝑦1
1 − 𝑥1
1 − 𝑦1 𝑦1
∗𝛼 =
1 − 𝑥1 𝑥1

(1 − 𝑦1 )𝑥1 𝛼 = 𝑦1 (1 − 𝑥1 )

𝛼𝑥1 − 𝛼𝑥1𝑦1 = 𝑦1 − 𝑦1𝑥1

𝛼𝑥1 = 𝛼𝑥1𝑦1 + 𝑦1 − 𝑦1𝑥1
 𝛼𝑥1 = 𝑦1(𝛼𝑥1 + 1 − 𝑥1)

𝛼𝑥1
𝑦1 =
(1 + 𝑥1(𝛼 − 1))

Ley de Raoult
𝑃 = 𝑃1 𝑥1 + 𝑃2 𝑥2 = 𝑃1 𝑥1 + 𝑃2 (1 − 𝑥1 )
𝑃 = 𝑃1 𝑥1 + 𝑃2 −𝑃2 𝑥1
𝑃 − 𝑃2 = 𝑥1 (𝑃1 − 𝑃2 )
𝑃 − 𝑃2
𝑥1 =
𝑃1 − 𝑃2

L
T(k) P1(Pa) P2(Pa) x1 x2 (Kmol/h) V(Kmol/h) L+V
370.65 168029 68578 0.329 0.56 7.174 7.824 15

Polinomio para presiones de Vapor Antoine

Sustancia C1 C2 C3 C4 C5
Benceno 83.107 -6486.2 -9.2194 6.98E-06 2
Tolueno 76.945 -6729.8 -8.179 5.30E-06 2

2. Utilizando la ecuación de Antoine, calcule el punto de burbuja de una solución

ideal que contenga 40 %mol de benceno y 60 %mol de tolueno a 1 atm de
presión total. A, B y C para el Benceno son respectivamente 15.9037, 2789.01 y
220.79 y A, B, C para el Tolueno son respectivamente 16.00531, 3090.78 y
C=219.4
Datos
• xb=0.4
• xt=0.6
• PT=1 atm= 760 mmHg
Ecuación de Antoine

Sustancia A B C
Benceno 15.9037 2789.1 220.79
Tolueno 16.00531 3090.78 219.4

𝑦𝑖 = ∑ 𝐾𝑖𝑥𝑖 = 1

𝑃𝑖
𝐾𝑖 =
𝑃𝑡
𝑃𝑏 𝑃𝑡
𝑥𝑏 + 𝑥𝑡 =1
𝑃𝑇 𝑃𝑇
𝐵
Sustituyendo Pb y Pt por 𝑒 (𝐴− 𝐶+𝑇) respectivamente obtenemos que T= 95.01 °C

3. Utilizando la ecuación de Antoine, calcule el punto de roció de un vapor que

contiene 65% en moles de benceno y 35% en moles de tolueno a 1.5 bar de
presión total, donde:

Sustancia A B C
Benceno 15.9037 2789.1 220.79
Tolueno 16.00531 3090.78 219.4

Empleando el Simulador COFE de COCO Simulator se realizó la siguiente resolución:

a. Con la herramienta TEA se creó un paquete de propiedades con los
compuestos y el modelo para la resolución.
b. Se selecciono el benceno y el tolueno, para el modelo matemático se escogió
SOAVE Redlich Kwong.

c. Insertar una línea de corriente con el benceno y el tolueno.

d. Se colocan los datos proporcionados por el ejercicio, como la presión, la

fracción molar de cada compuesto, asumiendo 1 mol/s de flujo y una
temperatura de 25 °C.
e. Se generan automáticamente los datos de corriente, se seleccionó “Overall
Properties” donde puede encontrarse el punto de roció. 𝑇=107.93 ℃ 𝑦
𝑃=0.06926 bar.

Cálculos en Excel

Datos
P(bar) 1.5
P(Kpa) 150
P(atm) 1.4804
T(°C) 107.83

Empleando la herramienta de análisis de hipótesis en Excel.

Ecuación de Antoine
 Punto de Roció

∑ 𝑦𝑖/𝑘𝑖 = 1

𝑃𝑠𝑎𝑡
𝐾𝑖 =
𝑃𝑇
Cálculos Benceno Tolueno
Psat (Kpa) 221.93 93.64
Psat (atm) 2.19 0.92
Ki 1.479 0.624
yi 0.65 0.35
Punto de
Rocío 0.439 0.56
Σ 1

Grafica T (°C) vs Fracción Molar de Benceno.

4. Una mezcla de 35 mol% n-hexano, 30 mol% n-heptano y 35 mol% n-octano a

350 K y 1.5 atm de presión total. ¿Se separa la mezcla en dos fases? Si es así,
calcular la cantidad de los productos en fase de vapor. Asumir que la solución es
ideal.
Se realizó una simulación de la operación de destilación flash de una mezcla
multicomponentes a temperatura y presión especificada mediante la herramienta
de chemsep del Coco Simulator.

Configurar los detalles de la operación necesarios para la simulación.

1. Ingresar los componentes de la alimentación.

2. Seleccionar los detalles de la operación.
3. Propiedades:
• Utilizar la ley de Raoult con la suposición de una mezcla ideal ya
que son sustancias semejantes.
• Utilizar la ecuación de Antoine para modelar las presiones de vapor.

4. Especificar los detalles de alimentación.

 Resultados

5. Una mezcla equimolar de benceno y tolueno se somete a destilación flash a una

presión de 1 bar en el separador. Determine las composiciones (en fracción molar
de benceno) del líquido y del vapor que salen del separador cuando la alimentación
está vaporizada en un 30%. Estimar la temperatura en el separador. Los datos de
equilibrio para el sistema benceno-tolueno a 1 bar se indican en la siguiente tabla.

Los datos de equilibrio para el sistema benceno-tolueno a 1 bar se indican en la

siguiente tabla.

Fracción Fracción molar

molar de de benceno en
T(°C)
benceno en vapor
líquido
80 1.0 1.0
85 0.780 0.9
90 0.58 0.78
95 0.41 0.64
100 0.26 0.45
105 0.13 0.26
110 0.0 0.0
Cuáles son las concentraciones en el vapor y en el líquido, y la temperatura de
separación de la alimentación:

• 0% vaporizado
• 50% vaporizado
• 75% vaporizado
• 100%vaporizado

Graficar las líneas de operación para cada uno de estos casos y analizar cómo
cambian las líneas de operación al cambiar f.

Datos:

Componente T. Ebullición
(ºC)
Benceno (a) 80.1
Tolueno (b) 110.6

P=1 bar= 100.007 Kpa

Balance Global:
F = L +V
Balance Por Componente:
Fa= Xa*L + Ya*V

Asumiendo F=1 mol de alimentación: V= 1- L

Aplicando Ley de Raoult ideal: Ya ∗ P = Xa ∗ Pa(T)

P = Xa ∗ Pa + Xb ∗ Pb

Regla de la Palanca: L/Zal. = (y − z) / (y − x) considerando Zal. =1

Donde la Presión de saturación de cada componente la calculamos con: La
ecuación de Antoine:

Fracción mol: Xa +Xb=1

𝑋𝑎 = (P-Pb) /(Pa-Pb)

Sustancia A B C
Benceno
137,819 2726,81 217,572
(a)
Tolueno
13,932 3056,96 217,625
(b)

Resultados

Evaporación
T(°C) xa xb
(%)
0 91.403 0.5 0.5
30 93.58 0.4334 0.5666
50 94.98 0.3888 0.6112
75 96.7 0.3363 0.6637
100 98.2785 0.2904 0.7096
 115 EQUILIBRIO BENCENO -TOLUENO

110

105
Xa(Benceno)
Ya(Benceno)
100
Lineal (Xa(Benceno))
Lineal (Ya(Benceno))
95

90
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Xa, Ya

Análisis de Resultados.

El trabajo desarrollado en el primer ejercicio permitió determinar las propiedades de

mezcla en el destilador flash. Empleando la ecuación de Antoine para determinar las
presiones de vapor y haciendo uso de la regla de la palanca para determinar
porcentaje de condensado y vapor teniendo así un resultado que nos proporciona
un mayor porcentaje de Tolueno en el condensado comprobando así que el Tolueno
es de los él menos volátil.

L
T(k) P1(Pa) P2(Pa) x1 x2 (Kmol/h) V(Kmol/h) L+V
370.65 168029 68578 0.329 0.56 7.174 7.824 15
También se empleo la ley de Raoult idealizada para encontrar el punto de burbuja
de la mezcla benceno/tolueno dando resultados aceptables puesto que la estructura
de ambas sustancias es muy parecida formando una solución ideal por lo que es
aplicable, si estas sustancias fueran diferentes estructuralmente los resultados
empleando esta ley serian poco confiables y descartables.

Empleando el simulador COFE de COCO Simulator y Excel para realizar

comparaciones obtuvimos los siguientes resultados.

En este programa seleccionamos que como método de resolución se utiliza la

ecuación de estado de Soave Redlich Kwong como ecuación de gases para la
resolución de la práctica asumimos 1 mol/s de flujo y una temperatura de 25̊C lo
cual el programa nos dio los puntos de roció. 𝑇=107.93 ℃ 𝑦 𝑃=0.06926 𝑏𝑎𝑟.

Analizando el grafico nos arrojó el programa para esta mezcla se pudo notar lo que
nos comprueba la teoría de que al ser estructuras similares entonces muestra un
comportamiento de una solución ideal este grafico este hecho a una presión de
𝑃=0.06926 𝑏𝑎𝑟. Para la siguiente parte de nuestro trabajo se analizó una mezcla de
tres compuestos n-hexano, n-heptano y n-octano a una respectiva composición. Se
realizó una simulación de la operación de destilación flash de una mezcla
multicomponentes a temperatura y presión especificada mediante la herramienta de
chemsep del Coco Simulator. Este ejercicio se basó en la ley de Raoult con la
suposición de una mezcla ideal ya que son sustancias semejantes, Se utilizó la
ecuación de Antoine para modelar las presiones de Vapor. Al hacer los cálculos en
el programa se pudo determinar que toda la solución se encontraba en fase liquida
comprimida por lo tanto las fracciones molares y su concentración debería de estar
con la concentración de alimentación.

En el último punto de nuestro trabajo se nos pide que determinemos las

composiciones (en fracción molar de benceno) del líquido y el vapor que salen del
separador cuando la alimentación está vaporizada a diferentes porcentajes que
serían nuestros sub- incisos de una mezcla benceno tolueno este lo realizamos con
la formula. Con la cual se obtuvieron resultados bastante aceptables para los
distintos porcentajes de vaporizado.

Conclusiones.

• Se logro alcanzar los objetivos de destilación que coinciden con las

propiedades físicas de las sustancias de la mezcla al ser el benceno el
componente mas rico en el vapor pues es el que posee la presión de vapor
mas alta (a cierta presión) lo que lo convierte en el componente más volátil
entre el benceno y el tolueno.

• Con los resultados encontrados podemos deducir que la destilación flash es

un proceso muy eficaz por lo que es sumamente utilizado sobre todo cuando
este es empleado para una sola fase donde se puede tener un mejor control
en la producción as como la disminución de gastos y perdidas que pueda
presentar este proceso.
Bibliografía
• McCabe, W. L., Harriott, P., & Smith, J. C. (2007). Operaciones unitarias en
ingeniería química. México: McGraw-Hill.
• MHRD. (n.d.). Recuperado de
http://vlabs.iitkgp.ernet.in/cpd/exp3/index.html#
• Smith, J. M., C., V. N. H., & Abbott, M. M. (2015). Introducción a la
termodinámica en ingeniería química. México, D.F.: McGraw-Hill
Interamericana.

Anexos.

Ejercicio Excel 1

Ejercicio Excel 2
Ejercicio Excel 3

Ejercicio Excel 4

Ejercicio Excel 5

30% Vaporizado
Tasum.
Pa[Kpa] Pb[Kpa] Xa Ya L V
(ºC)
92 1,445,431 579,103 0,4859 0,7023 0,9349 0,0651
94 1,529,512 616,986 0,4198 0,6421 0,6392 0,3608
93,5 1,508,146 607,335 0,4360 0,6575 0,7110 0,2890
93,7 1,516,665 611,181 0,4295 0,6513 0,6821 0,3179
93,58 1,511,549 608,871 0,4334 0,6550 0,6994 0,3006
0% Vaporizado
Tasum. (ºC) Pa[Kpa] Pb[Kpa] Xa Ya L V
92 1,445,431 579,103 0,4859 0,7023 0,9349 0,0651
91,403 1,421,039 579,103 0,5000 0,7105 10,000 0,0000

50% Vaporizado

Tasum.
Pa[Kpa] Pb [Kpa] Xa Ya L V
(ºC)
94,98 1,572,064 636,255 0,3888 0,6111 0,4998 0,5002

75% Vaporizado

Tasum.
Pa [Kpa] Pb [Kpa] Xa Ya L V
(ºC)
96,7 1,648,948 671,229 0,3363 0,5546 0,2500 0,7500

100% Vaporizado

Tasum.
Pa [Kpa] Pb [Kpa] Xa Ya L V
(ºC)
982785 1,722,023 704,653 0,2904 0,5000 0,0000 10000