VAPOR- LIQUIDO EQUILIBRIO

¿Ha pasado alguna vez por delante de una refinería y se ha preguntado qué ocurre en esas altas torres?
Algunas de ellas son torres de destilación que se utilizan para separar una mezcla de productos químicos
en dos o más corrientes, cada una relativamente pura de una de las sustancias químicas. El proceso que
rige esta separación es el equilibrio vapor-líquido. Se calcula que que el 10% de la energía utilizada
comercialmente en Estados Unidos se emplea en la destilación. Por lo tanto, es importante que este
proceso sea lo más eficaz posible.

Tome una mezcla de dos o más sustancias químicas en un régimen de temperatura en el que ambas
tengan una presión de vapor importante. La composición de la mezcla en el vapor es diferente a la del
líquido. Aprovechando esta diferencia, se pueden separar dos sustancias químicas, lo que constituye la
base de la destilación. Sin embargo, para calcular este fenómeno es necesario predecir magnitudes
termodinámicas como la fugacidad y, a continuación, realizar balances de masa y energía en el sistema.
En este capítulo se explica cómo predecir las propiedades termodinámicas y, a continuación, cómo
resolver las ecuaciones de una separación de fases. Aunque la separación de fases es sólo una parte del
proceso de destilación, es la base de todo el proceso. En este capítulo aprenderá a resolver problemas
de equilibrio vapor-líquido, y estos principios se emplean en los cálculos de las torres de destilación de
los capítulos 6 y 7. Los problemas de equilibrio vapor-líquido se expresan como ecuaciones algebraicas, y
los métodos utilizados son los mismos que se introducidos en el Capítulo 2

FLASH Y SEPARACIÓN DE FASES

La presión parcial del agua en el aire (en equilibrio) será igual a la presión de vapor del agua a esa
temperatura. Ahora se calienta el recipiente. La presión de vapor aumenta, ya que se incrementa a
medida que aumenta la temperatura y, por tanto, la presión parcial del agua también aumenta. Si la
presión parcial del agua a la temperatura de la olla supera la presión parcial del agua en ambiente
(normalmente fijada por la humedad), el agua se evaporará.

Existe otro escenario que, por desgracia, es puramente imaginario. En ese escenario, a medida que se
incrementa la temperatura, la sustancia química uno se evapora por completo cuando la temperatura
alcanza su punto de ebullición. A medida que la temperatura sigue aumentando, se alcanza el punto de
ebullición de la segunda sustancia química y se evapora toda. ¡Ojalá ocurriera así! A pesar de que en la
vida real la separación de las dos sustancias químicas nunca es completa ni en fase líquida ni en fase
vapor, sigue siendo un fenómeno útil y constituye la base de la destilación.

Este capítulo examina primero las ecuaciones que gobiernan un flash isotérmico, y luego muestra cómo
puede predecir las cantidades termodinámicas que necesita para resolver el problema del destello
isotermo. isotérmico. Los problemas son todos conjuntos de ecuaciones algebraicas, y puede resolver
estos problemas
con Excel y MATLAB. A continuación, el capítulo aborda situaciones más complicadas separaciones
vapor-líquido, pero ahora utilizando Aspen Plus debido a su gran base de datos

FLASH ISOTÉRMICO - DESARROLLO DE ECUACIONES

Considere el diagrama de flujo de la figura 3.1. Supongamos que conoce la temperatura, la presión y la
composición general de la corriente de entrada. Las fracciones molares de los productos químicos en la
entrada se denominan {zi}. En el separador de fases, en cambio, el líquido y el vapor están separados. La
fracción molar de los productos químicos en la fase vapor se denomina {yi} y la de los productos
químicos en la fase líquida se denomina {xi}. fase líquida se denominan {xi}. Cuando el vapor y el líquido
están en equilibrio, se pueden relacionar las fracciones molares de cada sustancia química en el vapor y
el líquido mediante la ecuación

Yi = Kixi

Los miembros del conjunto {Ki} se denominan valores K, y pueden predecirse a partir de la
termodinámica, como se muestra a continuación. Por ahora, sin embargo, supongamos que conocemos
sus valores. Para obtener la ecuación que rige el fenómeno, primero hay que sumar las fracciones
molares de vapor y líquido de todos los componentes. vapor y líquido de todos los componentes.
Entonces, resta esas dos ecuaciones

A la expresión de equilibrio se le sustituye y se reordena para obtener:

A continuación, haz un balance de masas para cada componente sobre el separador de fases. F es el
flujo total, V es el flujo molar del vapor y L es el flujo molar del líquido. del líquido. El balance molar es
entonces

Dividir por F y definir vo como la fracción de la alimentación que es vapor:

Usando de nuevo la expresión de equilibrio, se puede escribir como

Resuelve las fracciones molares en el líquido:

Fracciones molares y valores K para el ejemplo vapor-líquido

Fracción molar Equilibrio de fase Valor K

Introduciendo esta expresión en la Ec. (3.4) se obtiene la ecuación final

Esto se denomina ecuación de Rachford-Rice. Observe que, si los valores K y las composiciones de
entrada {zi} son conocidos, esta es una ecuación no lineal para resolver para v0. Por lo tanto, puede
aplicar aquí los mismos métodos utilizados con Excel y MATLAB en el Capítulo 2. Una vez conocido el
valor de v0 conocido, puede calcular el valor de las composiciones líquidas, {xi}, y de vapor, {yi},
utilizando las Ecs. (3.1) y (3.3). El balance molar está entonces completo.

Ejemplo con Excel

Suponga que tiene una mezcla de hidrocarburos en la corriente de entrada. Quieres encontrar la
fracción de la corriente que es vapor y la fracción molar de cada producto químico en las corrientes de
vapor y líquido. La tabla 3.1 muestra los caudales y los valores K a 1808F y 70 psia, utilizando una base
de un mol por unidad de tiempo. Para resolver este problema, prepare la hoja de cálculo como se
muestra en la Tabla 3.2.

Paso 1 'Término1' es el numerador en la Ec. (3.9) y 'término2' es el denominador, y éstos se calculan en

las columnas D y E.

Paso 2 Se calculan los cocientes de término1 sobre término2 y se ponen en la columna F; esta columna
se suma en la celda F9, que representa la Ec. (3.9).
Ejemplo de hoja de cálculo para Flash

Paso 3 A continuación, puede utilizar Goal Seek para hacer que la ecuación de Rachford-Rice (celda F9)
sea cero

cambiando la fracción de la alimentación que es vapor (celda D1), dando el resultado mostrado. Una vez
que Una vez que se encuentra la fracción de vapor, las fracciones molares en las dos fases son fáciles de
calcular usando (3.8) y (3.1), y éstas se incluyen en las columnas G y H.

¿Cómo puedes comprobar estos resultados?

Paso 1 Por supuesto, el aspecto más importante es tener los valores correctos de K, que sólo se pueden
determinar comparándolos con los valores de K de la ecuación (3.8). con los datos experimentales. Todo
el cálculo depende de esos valores K.

Paso 2 Necesitas que la ecuación de Rachford-Rice sume cero, lo que es evidente en la hoja de calculo.

Paso 3 Puedes comprobar que la suma de las fracciones molares de vapor es igual a uno y que la suma
de las fracciones molares de líquido es igual a uno.

fracciones molares de líquido es igual a uno.

Paso 4 También puede comprobar un término de la ecuación de Rachford-Rice mediante cálculos

detallados, y luego copiar la ecuación. Esto asegura que la fórmula es correcta para todos componentes.

PARAMETROS DE TERMODINAMICA

¿De dónde procede el valor K de equilibrio vapor-líquido? Se define como

donde Yi es el coeficiente de actividad, f0i es la fugacidad de la sustancia química pura øfi es el
coeficiente de fugacidad en fase vapor y p es la presión total. La mayoría de estas cantidades pueden
calcularse utilizando la termodinámica. Por ejemplo, su libro de texto de termodinámica deriva la
siguiente ecuación para la fugacidad de una sustancia pura:

Esto implica que, si conoces el comportamiento p-V-T de una sustancia química, puedes hallar la
fugacidad calculando la integral. En este texto no es necesario hacerlo, porqueAspen Plus lo hará por ti.
En su lugar, puede utilizar los valores K de la bibliografía o realizar los cálculos con un simulador de
procesos.

El otro parámetro importante de la Ec. (3.10) es el coeficiente de actividad. El coeficiente de actividad

depende de las sustancias químicas implicadas, especialmente de si son polares o no (tienen un
momento dipolar). Existen varias opciones:

Opción 1 - si la solución es ideal, el coeficiente de actividad es 1,0.

Opción 2 - la teoría de la solución regular de Scatchard-Hildebrand funciona mejor para compuestos no

polares y compuestos con pesos moleculares similares, y puede funcionar para compuestos débilmente
polares

Opción 3 - la ecuación de Wilson es más robusta, y funciona para compuestos polares y puede manejar
azeótropos, pero no dos fases líquidas. Requiere parámetros de interacción de pares que deben
determinarse a partir de los datos.

Opción 4 - la ecuación no aleatoria de dos líquidos (NRTL) es más robusta y funciona para compuestos
polares y puede tratar azeótropos, pero no dos fases líquidas. compuestos polares y puede manejar
azeótropos así como dos fases líquidas.

Opción 5 - el método de contribución de grupo UNIFAC construye el coeficiente de actividada partir de

parámetros correlacionados con fragmentos moleculares específicos. Funcionaen los mismos casos que
la ecuación NRTL, pero UNIFAC puede predecir coeficientes de actividad en ausencia de datos
específicos, siempre que todos los fragmentos moleculares se hayan correlacionado a partir de otros
datos.se hayan correlacionado a partir de otros datos. en los mismos casos que la ecuación NRTL, pero
UNIFAC puede predecir coeficientes de actividad en ausencia de datos específicos, siempre que todos
los fragmentos moleculares se hayan correlacionado a partir de otros datos. se hayan correlacionado a
partir de otros datos.
El capítulo 6 ilustra otras opciones disponibles en Aspen Plus. El truco está en saber qué opción elegir
para su sistema químico. Desgraciadamente, la única forma de saber con que ha hecho una buena
elección es comparar sus predicciones con los resultados experimentales. experimentales. Para la
mayoría de estas opciones, el coeficiente de actividad depende de las fracciones molares de líquido, y
eso significa que los valores K también dependen de las fracciones molares de líquido. Por lo tanto, no
hay un único valor K que se pueda utilizar en los cálculos. En su lugar, una opción es elegir las fracciones
molares de líquido, determinar los coeficientes de actividad, luego los K y ver si las fracciones molares
de líquido cambian significativamente. Si es así, podría utilizar los nuevos valores y repetir el proceso. Si
este procedimiento no funcionara, habría que para resolver la ecuación de Rachford-Rice utilizando
simuladores de procesos como Aspen Plus. simuladores de procesos como Aspen Plus. Sin embargo, hay
un caso especial que es fácil. Si la solución es ideal, los coeficientes de actividad son 1,0. Si la presión es
lo suficientemente baja (digamos menos de 10 atm), la fase vapor es esencialmente un gas ideal. En ese
caso, el coeficiente de fugacidad del vapor también es 1,0. La fugacidad del líquido en el estado estándar
es la presión de vapor, también, a bajas presiones. Con estos supuestos, el valor K pasa a ser

donde vpi es la presión de vapor a una temperatura dada. En este caso especial, se pueden hallar los
valores K conociendo la presión de vapor en función de la temperatura y la presión total del sistema.
Evidentemente, los componentes ligeros (puntos de ebullición más bajos) tienen un K mayor que los
componentes pesados (puntos de ebullición más altos) y, por tanto, su valor K será mayor. Si coexisten
dos fases en equilibrio, al menos uno de los componentes necesita un valor K mayor que 1 y al menos
uno de los otros componentes necesita un valor K menor que uno. Aun así, no se garantiza la existencia
de dos y el resultado depende tanto de la composición como de los valores K. valores K.

Ejemplo con MATLAB

Puede resolver el mismo problema utilizando MATLAB. En primer lugar, defina una función que
represente el problema que quieres resolver [Ec. (3.9)]. A continuación, comprueba la función para que
es correcta. Por último, utiliza la función 'fzero' para encontrar la solución. La función que define el
problema es:

{%equilibrio vapor-líquido}

función y=vpequil(v)

z=[0.1 0.3 0.4 0.2]

K=[6.8 2.2 0.8 0.052]

suma1=0.;

para i=1:4

num=(K(i)-1)*z(i)

denom=1+(K(i)-1)*v

sum1=sum1+num/denominador

y=suma1}

Paso 1 Programe la función y guárdela en su directorio de trabajo como vpequil.m.

Paso 2 Emita el comando 'vpequil(0.2)' en la ventana de comandos; la salida debería ser:

z=0.1000 0.3000 0.4000 0.2000

K=6.8000 2.2000 0.8000 0.0520

numero=0.5800

denominador=2.1600

suma1=0.2685

....

ans=0.2415

El cálculo con calculadora da las mismas respuestas, por lo que has comprobado que el programa es
correcto.

Paso 3 Añade puntos y comas al final de cada línea del fichero m-fila y guárdalo. Resuelva el problema
emitiendo el comando fzero en el comando

>>fzero(@vpequil,0.2)
ans=0.4258

Como era de esperar, el resultado es el mismo que el obtenido utilizando Excel.

Ejemplo con Aspen Plus

Aspen Plus permite resolver fácilmente este mismo problema utilizando el funcionamiento de la unidad Flash2.
Paso 1 Inicie Aspen y elija Plantilla; Aceptar. Cuando aparezca la ventana, elija General con Unidades inglesas. En
Tipo de ejecución (esquina inferior derecha), elija Hoja de flujo. Haga clic en OK cuando aparezca la ventana del
motor Aspen. (Este último paso es específico de su instalación).

Paso 2 Si la parte inferior de la pantalla no muestra las unidades, utilice el menú desplegable Ver/Biblioteca de
modelos o pulse la tecla F10. En las unidades de la parte inferior, seleccione Separaciones y, a continuación, haga
clic en Flash2.

Paso 3 Haga clic en el diagrama de flujo y aparecerá una unidad flash (separación de fases). Para añadir los flujos de
entrada y salida, haga clic en Material Streams (esquina inferior izquierda), haga clic en el diagrama de flujo y
arrastre un flujo a la flecha roja que es la entrada a la unidad flash. Haga clic en las flechas rojas que salen y arrastre
el flujo fuera. Nota: Para realizar cambios en la ubicación de los flujos o unidades, puede hacer clic en la flecha
situada justo encima del botón Material. Puede alternar entre la flecha y Material Stream en la esquina inferior
izquierda a medida que mejora la presentación de su diagrama de flujo. Cuando se muestre el diagrama de flujo,
haga clic en Corrientes de materiales. Si aparece alguna flecha roja en el diagrama, significa que la unidad no está
bien conectada; hace falta una corriente de entrada o de salida. Corríjalo antes de continuar - puede que haya
colocado el cursor incorrectamente al dibujar los flujos. El resultado debería parecerse a la Figura 3.2.

Paso 4 Haga clic en las gafas. Aparecerá un menú a la izquierda de la pantalla. Las casillas

rojas indican que aún tiene que facilitar información. Empiece por arriba y vaya bajando hasta convertir las casillas
rojas en azules rellenando los formularios.

Paso 5 En la lista de la izquierda, seleccione Componente/Especificaciones e identifique los componentes como hizo
en el capítulo 2 (véase la figura 2.3). Escriba el nombre o la fórmula de los productos químicos. Si Aspen Plus no
reconoce el producto químico, aparecerá una ventana que le permitirá realizar una nueva búsqueda y le sugerirá una
serie de posibilidades. Cuando los componentes estén completamente especificados, deberá tener una entrada para
cada sustancia química en la columna denominada "Nombre del componente". La primera columna es la que da
nombre a las sustancias químicas, pero la tercera es la que utiliza Aspen Plus cuando obtiene las propiedades físicas.
Si esa columna está en blanco, el programa no funcionará.
Paso 6 En la lista de la izquierda, seleccione Propiedades/Especificaciones y elija el modelo termodinámico como
hizo en el capítulo 2 (véase la figura 2.4). En método de propiedades, desplácese hacia abajo hasta obtener RK-
Soave. También debe indicar al ordenador qué base de datos debe utilizar para la interacción

parámetros en RK-Soave. Elija una de las opciones para los parámetros de interacción. Los resultados mostrados a
continuación utilizan RKSKBV-1 de la tabla SDF.

Paso 7 En la lista de la izquierda, elija Streams haciendo doble clic sobre ella. Dentro de esa carpeta hay dos o más
carpetas, una para cada flujo. Elija la corriente de entrada, haga clic sobre ella, e introduzca la temperatura, la
presión y los caudales de cada producto químico, como se muestra en la Figura 3.3. Tenga en cuenta que puede
especificar las unidades que desea utilizar a medida que introduce los núeros. De este modo, puede utilizar los
números tal y como se le presentan, sin tener que hacer usted mismo la conversión de unidades. de unidades.

Paso 8 En la lista de la izquierda, elija Bloques, luego B1 (o como haya llamado a su unidad Flash2), luego Entrada.
Seleccione Especificaciones e introduzca la temperatura y la presión de la unidad como se muestra en la Figura 3.4.
(Para este problema, puede utilizar la misma temperatura y presión que la corriente de entrada, y utilizar la misma
temperatura y presión que la corriente de entrada). presión que la corriente de entrada, y considerar la unidad
simplemente como un separador de fases).

parámetros en RK-Soave. Elija una de las opciones para los parámetros de interacción. Los resultados que se
muestran a continuación utilizan RKSKBV-1 de la tabla SDF.

Paso 7 En la lista de la izquierda, elija Streams haciendo doble clic sobre ella. Dentro de esa carpeta hay dos o más
carpetas, una para cada flujo. Elija la corriente de entrada, haga clic sobre ella, e introduzca la temperatura, la
presión y los caudales de cada producto químico, como se muestra en la Figura 3.3. Tenga en cuenta que puede
especificar las unidades que desea utilizar a medida que introduce los números. De este modo, puede utilizar los
números tal y como se le presentan, sin tener que hacer usted mismo la conversión de unidades.

Paso 8 En la lista de la izquierda, elija Bloques, luego B1 (o como haya llamado a su unidad Flash2), luego "Input".
Seleccione Especificaciones e introduzca la temperatura y la presión de la unidad como se muestra en la Figura 3.4.
(Para este problema, puede utilizar la misma temperatura y presión que la corriente de entrada, y utilizar la misma
temperatura y presión que la corriente de entrada). y considere la unidad simplemente como un separador de fases).

Paso 9 Elija la en la parte superior. En la figura 3.4 puede ver que hay dos opciones. Si la entrada está incompleta,
aparecerá una ventana para avisarle y dirigirle a los datos que faltan. Si la entrada está completa, aparecerá una
ventana para avisarle de ello, también. Haga clic en el botón para que realice el cálculo. Esto hará que se proceda al
cálculo del proceso (aquí una unidad). Una vez finalizados los cálculos (lea los mensajes de error, si los hubiera),
haga clic en la casilla Resultados (la inferior) para volver al menú habitual. A partir de ahí, busque en
Resultados/Corrientes. Los datos del flujo aparecerán en forma de tabla, como se muestra en la Tabla 3.3. Si hace
clic en Stream Table, también se reproducirán en su diagrama de flujo

Paso 10 Para su informe, seleccione el esquema del proceso colocando el cursor en la esquina superior izquierda y
arrastrándolo hasta la esquina inferior derecha. Elija Copiar, cambie a un programa de procesamiento de textos y
elige Pegar. Haga lo mismo con la tabla de flujo haciendo clic sobre ella. Como alternativa, coloque el cursor en la
celda superior izquierda de la tabla Stream copie hasta la celda inferior derecha, copie y pegue en Excel. Su informe
escrito debe incluir algo más que el esquema del proceso y la tabla de flujo. Debe explicar qué problema ha resuelto
y cómo lo ha hecho, centrándose en la información de ingeniería química, no en los pasos detallados de la pantalla
del ordenador. Describe cómo has comprobado los resultados e indícalos. En este caso, has tomado una corriente de
composición, presión y temperatura especificadas que estaba en dos fases, has separado el gas del líquido y has
calculado la composición de cada fase. Dependiendo de tu nivel de conocimientos en ingeniería química, también
deberás especificar y justificar el tipo de modelo termodinámico que has elegido. Puede suponer que Aspen Plus
realizó el cálculo correctamente, pero también debe examinar los resultados para ver si son razonables.

En la tabla 3.3 se puede ver que la fracción de vapor es 0,419 (frente a 0,426 en la tabla 3.2). Dado que Aspen Plus
ha calculado los valores K, es posible que haya utilizado valores diferentes a los suyos. Los resultados son los
siguientes: propano, 5,2; n-butano, 2,1; n-pentano, 0,84; n-octano, 0,067. Aunque estos valores son los mismos que
los de Aspen Plus, el resultado es el siguiente, Aunque estos números difieren ligeramente de los utilizados en la
Tabla 3.2, la fracción molar de cada componente químico es la misma.

en las corrientes de vapor y líquido es próximo. Los valores K también figuran en el informe (véase el Apéndice C).
LÍQUIDOS NO IDEALES - PRUEBA DEL MODELO TERMODINÁMICO

En los ejemplos anteriores, se le han dado los valores K o se le ha dicho qué modelo termodinámico elegir.Por lo
general, la decisión es suya, y sus resultados son tan buenos como su elección. Por lo tanto, es importante
preferiblemente con datos experimentales. Una forma sencilla de hacerlo es modelizar el diagrama T-xy o p-xy para
el vapor-líquido de pares binarios. Puede hacerlo fácilmente utilizando la opción Análisis de propiedades de Aspen
Plus. A continuación se indican los comandos que le permitirán crear estos diagramas. También puede ser necesario
revisar el Apéndice C, que tiene más detalles sobre Aspen Plus.

Paso 1 Inicie Aspen y seleccione Plantilla. Aparecerá una ventana; seleccione General con métricas. En Tipo de
ejecución (esquina inferior derecha), elija Análisis de propiedades. Haga clic en Aceptar cuando aparezca la ventana
del motor Aspen. (Este último paso es específico de su instalación).

Paso 2 En la lista de la izquierda, elija Componente/Especificaciones e introduzca los nombres

o fórmulas de los productos químicos. Si no hay lista a la izquierda, haga clic en las gafas o

elija el menú Datos/Componentes. Si Aspen Plus no reconoce su producto químico, aparece una Si Aspen Plus no
reconoce su producto químico, aparecerá una ventana que le permitirá buscar de nuevo y le sugerirá una serie de
posibilidades. Cuando los componentes estén completamente especificados, deberá tener una entrada para cada en la
columna "Nombre del componente". La primera columna es la que la tercera columna es la que utiliza Aspen Plus
para obtener las propiedades físicas. físicas. Si esa columna está en blanco, el programa no funcionará.

Paso 3 En la lista de la izquierda, seleccione Propiedad/Especificaciones. Establezca el método de la propiedad

como ilustrado en el Capítulo 2 (ver Figura 2.4). Para esta aplicación, desplácese hacia abajo y elija WILS-2. Esta
opción utiliza un gas ideal, la ley de Henry y los parámetros binarios de Wilson.

Paso 4 También debe indicar al ordenador qué base de datos debe utilizar para el parámetro de interacción. Elija una
de las opciones para los parámetros de interacción. Los resultados mostrados a continuación utilizan los parámetros
binarios de Wilson-2. Ahora todos los iconos de la izquierda deben ser de color azul, lo que significa que los menús
están completamente rellenados. Vuelva atrás y complete cualquier menú rojo.

Paso 5 En los menús de la parte superior elija Herramientas/Análisis/Propiedades/Binario. (La opción estará en gris
hasta que realice los pasos dos, tres y cuatro). En la ventana que aparece que aparece, como se ilustra en la Figura
3.5, elija lo siguiente:

Componentes: seleccione dos componentes

Presión: establecer la presión para un diagrama T-xy O

Temperatura: establece la temperatura para un diagrama p-xy

Haga clic en Ir.

Paso 6 Aparecerá un gráfico con el resultado. Puede leer en este gráfico o, si lo mueve ligeramente ligeramente el
gráfico, aparecerá una tabla con los valores tabulados. En
En la figura 3.6 se muestra el diagrama yx para la termodinámica ideal. No es un buen ajuste. Usted Usted sabe que
el etanol y el agua forman un azeótropo, y las predicciones no dan un azeótropo (donde y= x). azeótropo (donde
y=x). Debe concluir que una solución ideal no es un buen modelo termodinámico para este sistema. termodinámico
para este sistema.

Figura 3.6. Diagrama vapor-líquido etanol-agua para termodinámica ideal.

Los resultados de la termodinámica de Wilson se muestran en la Figura 3.7, y los cálculos concuerdan con los datos.
La tabla (detrás del gráfico en Aspen Plus) también enumera valores K y coeficientes de actividad. Cuando se tiene
una buena comparación con los datos, ya sea datos, ya sea en forma de diagrama T-xy, valores K o coeficientes de
actividad, se ha validado un modelo termodinámico para esta ecuación binaria.

modelo termodinámico para este par binario. En un proceso con muchas sustancias tener que hacer esto para todos
los pares binarios antes de decidir cuál es el mejor modelo termodinámico.

RESUMEN DEL CAPÍTULO

En este capítulo se han deducido las ecuaciones que rigen el equilibrio de fases y se ha visto cómo se pueden estimar
los parámetros clave utilizando la termodinámica. Ha resuelto los problemas

problemas resultantes utilizando Excel, MATLAB y Aspen Plus. También ha aprendido a preparar un diagrama T-
xy para comprobar el modelo termodinámico elegido para representar el fenómeno. fenómeno.

PROBLEMAS

3.1. La siguiente corriente está a 100 psia y 1788F. Calcule la fracción que es vapor resolviendo la ecuación de
Rachford-Rice: (1) usando Excel; (2) usando MATLAB.
(Aspen Plus) Un embalse del norte de Luisiana contiene un petróleo volátil. Las condiciones en el momento del
descubrimiento eran de 2468F, 4800 psia. La composición era

Tendrás que elegir cómo modelar los butanos y pentanos: ¿normal? ¿iso-? En También tendrás que elegir qué
sustancia química utilizar para modelizar los heptanos plus: ¿heptano? ¿octano? Utiliza la opción RK-Soave de
Aspen Plus.

(1) El separador gas-líquido en la superficie está a 500 psia, 658F. Hallar la composición de las corrientes gaseosa y
líquida y la fracción de vapor.

(2) El líquido se lleva al tanque de almacenamiento, que está a 14.7 psia, 708F. Encuentre la fracción de vapor

en estas condiciones.

3.3. La siguiente corriente sale de una torre de destilación (descrita detalladamente en el capítulo Está a 138 psia y
197,58F. Si la presión se reduce (adiabáticamente) a 51 psia, ¿cuál será la fracción de vapor y la temperatura?
(Sugerencia: En Aspen Plus, coloca una válvula antes de la unidad Flash2, y reduce la presión con el bloque de
válvulas).

3.4. La siguiente corriente proviene de una torre de destilación (descrita detalladamente en el capítulo 6). Está a 36
psia y 141,58F. Si la presión se reduce (adiabáticamente) a 20 psia, ¿cuál será la fracción de vapor y la temperatura?
(Sugerencia: En Aspen Plus, coloca una válvula antes de la unidad Flash2, y reduce la presión con el bloque de
válvulas).
3.5. Una corriente que sale de un reactor de amoníaco está a 7008F y 3000 psia. Los caudales se indican a
continuación

Si la temperatura se reduce a 808F en una unidad flash, manteniendo la misma presión, determine el caudal de cada
especie en la corriente de vapor y de líquido.

3.6. Los caudales de salida de un reactor de metanol son: CO, 100 lb mol/h; hidrógeno,200 lb mol/h; metanol, 100 lb
mol/h. Esta corriente está a 100 atm y 3008C.Enfríela a 308C y determine los caudales de cada especie en la
corriente de vapor y de líquido.

(1) ¿Cuáles son los valores K de las tres especies en estas condiciones? Una separación flash(Flash2 en Aspen) es
adecuada para modelizar el condensador. Utilizando las fracciones molares yi y xi obtenidas por Aspen Plus para el
vapor y el líquido respectivamente,calcula los valores K, Ki = yi/xi.

(2) Realice la simulación para una solución ideal y también para otro modelo termodinámico de su elección.
termodinámico (a su elección).

(3) El informe debe incluir un diagrama de flujo y los balances de masa. Discuta la diferencia entre los dos casos y
por qué crees que uno puede ser mejor.¿Cómo puedes averiguar cuál es más realista? Nota: Los momentos dipolares
de las tres especies son: CO, 0,1 debyes; hidrógeno, 0 debyes; metanol, 1,7 debyes (Reid et al., 1977).

3.7. Los caudales de salida de un reactor de equilibrio para producir hidrógeno fueron: CO, 5,25 hidrógeno, 1095;
dióxido de carbono, 699; y agua, 364,2 (utilícelos como lb mol/h).

Lo ideal sería enfriar la corriente y eliminar toda el agua. En la práctica agua, y el agua puede contener trazas de los
otros componentes. de los otros componentes. Este problema te reta a averiguar cantidad. La salida del reactor está a
450 K, 1 atm, y se enfría a 808F, todavía a 1 atm.

(1) ¿Cuáles son los valores K de las tres especies en estas condiciones? Una separación

(Flash2 en Aspen) es adecuado para modelar el condensador. Utilizando las fracciones molares yi y xi obtenidas por
Aspen Plus para el vapor y el líquido, respectivamente, calcula los valores K, Ki= yi/xi.

(2) Realice la simulación para una solución ideal y también para otro modelo termodinámico de su elección.
termodinámico (a su elección).

(3) El informe debe incluir un diagrama de flujo y los balances de masa. Discuta la diferencia entre los dos casos y
por qué crees que uno puede ser mejor,¿Cómo puedes averiguar cuál es más realista?

3.8. Utilizando la ecuación de Rachford-Rice, Ec. (3.9), demuestre que para que coexistan dos fases, al menos un
componente necesita un valor K mayor que uno y otro componente necesita un valor K menor que uno.
3.9. Modele los siguientes datos de equilibrio vapor-líquido de acetato de etilo y agua a 760 mmHg utilizando el
modelo WilsonRK en Aspen Plus (Ellis y Garbett, 1960):

3.10. Modele los siguientes datos de equilibrio vapor-líquido de alcohol etílico y acetato de etilo utilizando el
modelo WilsonRK en Aspen Plus (Murti y Van Winkle, 1959):

3.11. Modele los siguientes datos de equilibrio vapor-líquido de acetaldehído y alcohol etílico utilizando el modelo
WilsonRK en Aspen Plus (De Leeuw, 1911):

:D