Universidad San Francisco de Quito

Nombre: Sebastián Hidalgo 00326921

Fecha: 26/09/2023 NRC: 4253

Deber 2

Ejercicio 2.1.3

a) Calcule el volumen molar y el factor de compresibilidad para el amoniaco

gaseoso a una presión P=56 atm, y una temperatura T=450 K, a partir de la
ecuación de estado de van der Waals.

Ecuación de estado de van der Waals

𝑎
(𝑃 + ) (𝑉 − 𝑏) − 𝑅𝑡 = 0 (1)
𝑉2

Donde

27 2 𝑇𝑐 2
𝑎= (𝑅 ∙ ) (2)
64 𝑃𝑐

𝑅 ∙ 𝑇𝑐
𝑏=( ) (3)
8 ∙ 𝑃𝑐

Tenemos como constantes las siguientes expresiones:

𝑇𝑐 : Temperatura crítica (405,5 K para el amoniaco)

𝑃𝑐 : Presión crítica (111,3 atm para el amoniaco)

La presión reducida se define como

𝑃
𝑃𝑟 = (4)
𝑃𝑐

Y el factor de compresibilidad viene dado por

𝑃𝑉
𝑍= (5)
𝑅𝑇
Tabla 1. Constantes planteadas por el ejercicio.

Constantes
P (atm) 56
T (k) 450
R (atm L/g-
mol k) 0,08206
Tc (k) 405,5
Pc (atm) 111,3

Realizamos en Excel cálculos necesarios para hallar el volumen molar para el amoníaco.

Primero en una matriz 4x2 identificamos y colocamos en cada celda las variables de
volumen molar (V), constante (a) y (b) de la ecuación de van der Waals y en la última
columna la ecuación de van der Waals {f(v)}.

Tabla 2. Cálculos para hallar el volumen molar y el factor de compresibilidad del

amoniaco a P=56 atm y T=450 K.

V a b f(v) Z

Segundo paso es colocar un valor al azar para el volumen molar (V) para reemplazarlo en la
ecuación (1), calculamos también la constante (a) y (b) de la ecuación de van der Waals
denotadas por las ecuaciones (2) y (3) respectivamente, por último, reemplazamos el valor
de la constante (a), (b) y (V) en la ecuación (1) perteneciente a la ecuación de estadp de van
der Waals. Para calcular el volumen molar para el amoniaco a la presión y temperatura
establecidas, la ecuación de van der Waals debe ser igual a cero, para hallar este valor
utilizamos la opción ‘Análisis de hipótesis’, ‘buscar objetivo’, definimos la celda {f(v)}
con el valor 0 cambiando la celda (V). Para hallar el factor de compresibilidad se reemplaza
valores en la ecuación (5). De este modo obtenemos los valores resultantes mostrados en la
tabla 3. Para comprobar todo el proceso puede referirse al archivo de Excel denominado
‘Ejercicio 2.1.3’, en el mismo puede verificar que se utilizaron de la ecuación (1) a la (5).
Tabla 3. Resultados del volumen molar para el amoniaco a P=56 atm y T=450 K.

V a b f(v) Z
0,57488437 4,19694587 0,03737122 -0,00033912 0,87181533

b) Repita los cálculos para las siguientes presiones reducidas: 𝑃𝑟 = 1, 2, 4, 10 y 20.

Tabla 4. Constantes proporcionadas por el ejercicio para resolver la parte b).

Constantes
T (k) 450
R (atm L/g-
mol k) 0,08206
Tc (k) 405,5
Pc (atm) 111,3

Para resolver la parte b), el ejercicio nos proporciona la presión reducida como dato, en
Excel calculamos la presión (P) mediante el despeje de la fórmula (4), Procedemos a
realizar una matriz 6x6 en donde vamos a colocar en cada celda las variables de presión
reducida (𝑃𝑟 ), la presión calculada (P), el volumen molar (V), las constantes (a) y (b) y la
ecuación de van der Waals {f(v)}

Tabla 5. Cálculos para hallar el volumen molar para el amoniaco con las presiones
reducidas dadas.

Pr (atm) P (atm) V a b f(v)

Una vez tenemos los valores de presión (P), el segundo paso es colocar un valor al azar
para el volumen molar (V), procedemos hallar el valor para la ecuación (2) y (3) y de este
modo obtener la constante (a) y (b) respectivamente de la ecuación de van der Waals, por
último, reemplazamos los cuatro valores anteriores en la ecuación de van der Waals. Para
calcular el volumen molar para el amoniaco a la presión y temperatura establecidas, la
ecuación de van der Waals debe ser igual a cero, para hallar este valor utilizamos la opción
‘Análisis de hipótesis’, ‘buscar objetivo’, definimos la celda {f(v)} con el valor 0
cambiando la celda (V). De este modo obtenemos los valores resultantes mostrados en la
tabla 6. Para comprobar todo el proceso puede referirse al archivo de Excel denominado
‘Ejercicio 2.1.3’, en el mismo puede verificar que se utilizaron de la ecuación (1) a la (5)
así como verificar los datos obtenidos.

Tabla 6. Resultados del volumen molar para el amoniaco a T=450 K y P dada.

Pr (atm) P (atm) V a b f(v)

1 111,3 0,23352408 4,19694587 0,03737122 0,00091121
2 222,6 0,07727239 4,19694587 0,03737122 0,00096421
4 445,2 0,06065395 4,19694587 0,03737122 -0,00023731
10 1113 0,05087531 4,19694587 0,03737122 7,1365E-05
20 2226 0,04617492 4,19694587 0,03737122 -0,00044801

Ahora utilizando la ecuación (5) hallamos el factor de compresibilidad para cada una de las
diferentes presiones reducidas.

Tabla 7. Factor de compresibilidad para cada una de las presiones reducidas.

Pr (atm) Z
1 0,70385437
2 0,46580641
4 0,73125727
10 1,53340985
20 2,78347448
 c) ¿Cómo varía el factor de compresibilidad en función de 𝑷𝒓 .

Pr vs Z
3
y = 0,1195x + 0,3589
Factor de Compresibilidad

2,5 R² = 0,9776

1,5

0,5

0
0 5 10 15 20 25
Presión Reducida

Figura 1. Factor de compresibilidad en función de la presión reducida.

Ejercicio 2.2.3

Repita el Problema 2.1 utilizando la ecuación de estado de Redlich-Kwong.

a) Calcule el volumen molar y el factor de compresibilidad para el amoniaco

gaseoso a una presión, P=56 atm y una temperatura T=450 K, a partir de la
ecuación de estado de Redlich-Kwong.

Ecuación de estado de Redlich-Kwong

𝑅𝑇 𝑎
𝑃− + =0 (6)
(𝑉 − 𝑏) 𝑉 (𝑉 + 𝑏)√𝑇

Donde

𝑅2 ∙ 𝑇𝑐 5/2
𝑎 = 0,42747 ( ) (7)
𝑃𝑐

Y
 𝑅 ∙ 𝑇𝑐
𝑏 = 0,08664 ( ) (8)
𝑃𝑐

Tenemos como constantes las siguientes expresiones:

𝑇𝑐 : Temperatura crítica (405,5 K para el amoniaco)

𝑃𝑐 : Presión crítica (111,3 atm para el amoniaco)

La presión reducida se define como

𝑃
𝑃𝑟 = (4)
𝑃𝑐

Y el factor de compresibilidad viene dado por

𝑃𝑉
𝑍= (5)
𝑅𝑇

Tabla 8. Constantes planteadas por el ejercicio.

Constantes
P (atm) 56
T (k) 450
R (atm L/g-
mol k) 0,08206
Tc (k) 405,5
Pc (atm) 111,3

Realizamos en Excel cálculos necesarios para hallar el volumen molar para el amoníaco.

Primero en una matriz 4x2 identificamos y colocamos en cada celda las variables de
volumen molar (V), constante (a) y (b) de la ecuación Redlich-Kwong y en la última
columna la ecuación Redlich-Kwong {f(v)}.
Tabla 9. Cálculos para hallar el volumen molar y factor de compresibilidad del amoniaco a
P=56 atm y T=450 K.

V a b f(v) Z

El segundo paso es hallar las constantes (a) y (b) de la ecuación de estado de Redlich-
Kwong, procedemos a colocar un valor al azar para el volumen molar (V), este valor al azar
junto con las contantes (a) y (b) lo reemplazamos en la ecuación (6). Para calcular el
volumen molar para el amoniaco a la presión y temperatura establecidas, la ecuación de
Redlich-Kwong debe ser igual a cero, para hallar el volumen molar utilizamos la opción
‘Análisis de hipótesis’, ‘buscar objetivo’, definimos la celda {f(v)} con el valor 0
cambiando la celda (V). De este modo obtenemos los valores resultantes mostrados en la
tabla 3. Para comprobar todo el proceso puede referirse al archivo de Excel denominado
‘Ejercicio 2.2.3’, en el mismo puede verificar que se utilizaron de la ecuación (7) y (8) para
hallar las contantes (a) y (b) respectivamente y se reemplazó en la ecuación (6) así como
verificar los datos obtenidos.

Tabla 10. Resultados del volumen molar para el amoniaco a P=56 atm y T=450 K.

V a B f(v) Z
0,56981245 85,6348711 0,02590274 -0,00073481 0,864123741

b) Repita los cálculos para las siguientes presiones reducidas: 𝑷𝒓 = 1, 2, 4, 10 y

20.

Tabla 11. Constantes proporcionadas por el ejercicio para resolver la parte b).

Constantes
T (k) 450
R (atm L/g-
mol k) 0,08206
Tc (k) 405,5
 Pc (atm) 111,3

Para resolver la parte b), el ejercicio nos proporciona la presión reducida como dato, en
Excel calculamos la presión (P) mediante el despeje de la fórmula (4), Procedemos a
realizar una matriz 6x6 en donde vamos a colocar en cada celda las variables de presión
reducida (𝑃𝑟 ), la presión calculada (P), el volumen molar (V), las constantes (a) y (b) y la
ecuación de Redlich-Kwong {f(v)}

Tabla 12. Cálculos para hallar el volumen molar para el amoniaco con las presiones
reducidas dadas utilizando la ecuación Redlich-Kwong.

Pr (atm) P (atm) V a B f(v)

Una vez tenemos los valores de presión (P), el segundo paso es colocar un valor al azar
para el volumen molar (V), procedemos hallar el valor para la ecuación (7) y (8) y de este
modo obtener la constante (a) y (b) respectivamente de la ecuación de Redlich.Kwong, por
último, reemplazamos los cuatro valores anteriores en la ecuación de Redlich-Kwong. Para
calcular el volumen molar para el amoniaco a la presión y temperatura establecidas, la
ecuación de Redlich-Kwong debe ser igual a cero, para hallar este valor utilizamos la
opción ‘Análisis de hipótesis’, ‘buscar objetivo’, definimos la celda {f(v)} con el valor 0
cambiando la celda (V). De este modo obtenemos los valores resultantes mostrados en la
tabla 13. Para comprobar todo el proceso puede referirse al archivo de Excel denominado
‘Ejercicio 2.2.3’, en el mismo puede verificar que se utilizaron de la ecuación (6) a la (8)
así como verificar los datos obtenidos.
Tabla 13. Resultados del volumen molar para el amoniaco a T=450 K y P dada.

Pr (atm) P (atm) V a B f(v)

1 111,3 0,23327968 85,6348711 0,025902737 3,3816E-05
2 222,6 0,07371506 85,6348711 0,025902737 -2,1522E-05
4 445,2 0,05087635 85,6348711 0,025902737 1,187E-05
10 1113 0,03980301 85,6348711 0,025902737 0,00018569
20 2226 0,0348354 85,6348711 0,025902737 -0,00019519

Ahora utilizando la ecuación (5) hallamos el factor de compresibilidad para cada una de las
diferentes presiones reducidas.

Tabla 14. Factor de compresibilidad para cada una de las presiones reducidas.

Pr (atm) Z
1
0,70311774
2
0,44436247
4
0,61337646
10
1,19968465
20
2,09991622

c) ¿Cómo varía el factor de compresibilidad en función de 𝑷𝒓 .

 Pr vs Z
2,5

y = 0,0833x + 0,3956
Factor de Compresibilidad

2
R² = 0,9556

1,5

0,5

0
0 5 10 15 20 25
Presión Reducida

Figura 2. Factor de compresibilidad en función de la presión reducida.

Ejercicio 2.3.3

Tenemos la reacción de conversión biológica

𝐶𝐻𝑚 𝑂𝑛 + 𝑎𝑂2 + 𝑏𝑁𝐻3 → 𝑐𝐶𝐻𝛼 𝑂𝛽 𝑁𝛿 + 𝑑𝐻2 𝑂 + 𝑒𝐶𝑂2 (9)

Y tenemos los balances para cada uno de los elementos de la reacción:

Balance de carbono: 1= 𝑐+𝑒 (10)

Balance de hidrógeno: 𝑚 + 3𝑏 = 𝑐𝛼 + 2𝑑 (11)

Balance de oxígeno: 𝑛 + 2𝑎 = 𝑐𝛽 + 𝑑 + 2𝑒 (12)

Balance de nitrógeno: 𝑏 = 𝑐𝛿 (13)

El coeficiente respiratorio, CR, se define como:

𝑒
𝑅𝐶 = (14)
𝑎

a) El sustrato de glucosa, 𝑪𝟔 𝑯𝟏𝟐 𝑶𝟔, reacciona con oxígeno y amoníaco para

formar un material bacteriano, 𝑪𝑯𝟐 𝑶𝟎,𝟐𝟕 𝑵𝟎,𝟐𝟓, agua y dióxido de carbono con
 un coeficiente respiratorio de 1,5. ¿Cuáles son los coeficientes estequiométricos
para esta reacción cuando se escribe en la forma de la ecuación (9)?

El coeficiente respiratorio es 1,5 así que lo podemos reemplazar en la ecuación (14) y

despejar e.

𝑒
1,5 =
𝑎

𝑒 = 1,5𝑎 (15)

Las datos conocidos son: 𝑚 = 2, 𝑛 = 1, 𝛼 = 2, 𝛽 = 0,27 𝑦 𝛿 = 0,25.

Procedemos a reemplazar los datos conocidos en las ecuaciones (10) a (13), agregamos la
ecuación (15) y formamos un sistema de ecuaciones lineales, así obtenemos las nuevas
ecuaciones.

𝑐+𝑒 =1 (16)

3𝑏 − 2𝑐 − 2𝑑 = −2 (17)

2𝑎 − 0,27𝑐 − 𝑑 − 2𝑒 = −1 (18)

𝑏 − 0,25𝑐 = 0 (19)

1,5𝑎 − 𝑒 = 0 (20)

Ingresamos en Matlab y denominamos ‘Alpha’ a nuestra matriz 5x5, a nuestro vector lo

denominamos ‘Beta’ y procedemos a resolver el sistema de ecuaciones lineales sacando la
inversa de la matriz y multiplicando por el vector ‘Beta’.

En la figura 3 podemos observar la creación de la matriz 5x5 denominada ‘Alpha’ que

representa las incógnitas del problema 2.3.3 para los literales a) y b).
Figura 3. Valores de las incógnitas del problema de la fórmula de reacción química de
conversión biológica.

En la figura 4 podemos observar los valores del vector denominado ‘Beta’ el cual indica el
valor al que es igual los sistemas de ecuaciones de la matriz ‘Alpha’.

Figura 4. Valores de la igualdad del sistema de ecuaciones de la matriz ‘Alpha’

En la figura 5 podemos observar el valor de los coeficientes estequiométricos cuando el

coeficiente respiratorio es de 1,5.
Figura 5. Resultados de los valores estequiométricos cuando el coeficiente respiratorio es
de 1,5.

Reemplazando todas las incógnitas de la ecuación (9) tenemos como resultado la reacción
química de conversión biológica con un coeficiente respiratorio de 1,5.

𝐶𝐻2 𝑂1 + 0,2316 𝑂2 + 0,1631 𝑁𝐻3

→ 0,6525 𝐶𝐻2 𝑂0,27 𝑁0,25 + 0,5922 𝐻2 𝑂 + 0,3475 𝐶𝑂2 (21)

b) Repita los cálculos para la parte a) con un coeficiente respiratorio de 2,0.

El coeficiente respiratorio es 2,0 así que lo podemos reemplazar en la ecuación (14) y

despejar e.

𝑒
2,0 =
𝑎

𝑒 = 2𝑎 (22)

Formamos un nuevo sistema de ecuaciones utilizando las ecuaciones (16) a (19), esta vez
agregamos la ecuación (22) y formamos un nuevo sistema de ecuaciones lineales, así
obtenemos las nuevas ecuaciones.

𝑐+𝑒 =1 (16)

3𝑏 − 2𝑐 − 2𝑑 = −2 (17)

2𝑎 − 0,27𝑐 − 𝑑 − 2𝑒 = −1 (18)
 𝑏 − 0,25𝑐 = 0 (19)

2,0𝑎 − 𝑒 = 0 (22)

Ingresamos en Matlab y denominamos ‘Alpha2’ a nuestra matriz 5x5, a nuestro vector lo

denominamos ‘Beta2’ y procedemos a resolver el sistema de ecuaciones lineales sacando la
inversa de la matriz y multiplicando por el vector ‘Beta2’.

En la figura 6 podemos observar la creación de la matriz 5x5 denominada ‘Alpha2’ que

representa las incógnitas del problema 2.3.3.

Figura 6. Valores de las incógnitas del problema de la fórmula de reacción química de

conversión biológica.

En la figura 7 podemos observar los valores del vector denominado ‘Beta2’ el cual indica el
valor al que es igual los sistemas de ecuaciones de la matriz ‘Alpha2’.

Figura 7. Valores de la igualdad del sistema de ecuaciones de la matriz ‘Alpha2’

En la figura 8 podemos observar el valor de los coeficientes estequiométricos cuando el
coeficiente respiratorio es de 2,0.

Figura 8. Resultados de los valores estequiométricos cuando el coeficiente respiratorio es

de 2,0.

Reemplazando todas las incógnitas de la ecuación (9) tenemos como resultado la reacción
química de conversión biológica con un coeficiente respiratorio de 2,0.

𝐶𝐻2 𝑂1 + 0,1310 𝑂2 + 0,1845 𝑁𝐻3

→ 0,7380 𝐶𝐻2 𝑂0,27 𝑁0,25 + 0,5387 𝐻2 𝑂 + 0,2620 𝐶𝑂2 (23)

c) Repita los cálculos de la parte a) cuando el sustrato es ácido benzoico,

𝑪𝟔 𝑯𝟑 𝑪𝑶𝑶𝑯, y forma el mismo material bacteriano en condiciones anaerobias,
es decir, sin la presencia de oxígeno no gaseoso.

Tenemos una nueva reacción de conversión biológica en condiciones anaerobias.

𝐶6 𝐻5 𝐶𝑂𝑂𝐻 + 𝑏𝑁𝐻3 → 𝑐𝐶𝐻𝛼 𝑂𝛽 𝑁𝛿 + 𝑑𝐻2 𝑂 + 𝑒𝐶𝑂2 (24)

Y tenemos los balances para cada uno de los elementos de la reacción:

Balance de carbono: 7= 𝑐+𝑒 (25)

Balance de hidrógeno: 6 + 3𝑏 = 𝑐𝛼 + 2𝑑 (26)

Balance de oxígeno: 2 = 𝑐𝛽 + 𝑑 + 2𝑒 (27)

Balance de nitrógeno: 𝑏 = 𝑐𝛿 (28)

En el problema no necesitamos coeficiente respiratorio ya que son condiciones anaerobias.

Las datos conocidos son: 𝛼 = 2, 𝛽 = 0,27 𝑦 𝛿 = 0,25.

Procedemos a reemplazar los datos conocidos en las ecuaciones (25) a (28) y formamos un
sistema de ecuaciones lineales, obtenemos las nuevas ecuaciones:

𝑐+𝑒 =7 (29)

3𝑏 − 2𝑐 − 2𝑑 = −6 (30)

0,27𝑐 + 𝑑 + 2𝑒 = 2 (31)

𝑏 − 0,25𝑐 = 0 (32)

Ingresamos en Matlab y denominamos ‘Alpha3’ a nuestra matriz 4x4, a nuestro vector lo

denominamos ‘Beta3’ y procedemos a resolver el sistema de ecuaciones lineales sacando la
inversa de la matriz y multiplicando por el vector ‘Beta3’.

En la figura 9 podemos observar la creación de la matriz 4x4 denominada ‘Alpha3’ que

representa las incógnitas del problema 2.3.3 para el literal c).

Figura 9. Valores de las incógnitas del problema de la fórmula de reacción química de

conversión biológica en condiciones anaerobias.

En la figura 10 podemos observar los valores del vector denominado ‘Beta3’ el cual indica
el valor al que es igual los sistemas de ecuaciones de la matriz ‘Alpha3’.
Figura 10. Valores de la igualdad del sistema de ecuaciones de la matriz ‘Alpha3’

En la figura 11 podemos observar el valor de los coeficientes estequiométricos de la nueva

reacción química en condiciones anaerobias.

Figura 11. Resultados de los valores estequiométricos de la nueva reacción química en

condiciones anaerobias.

Reemplazando todas las incógnitas de la ecuación (24) tenemos como resultado la nueva
ecuación (33) de reacción química de conversión biológica en condiciones anaerobias.

𝐶6 𝐻5 𝐶𝑂𝑂𝐻 + 1,5924𝑁𝐻3 → 6,3694𝐶𝐻2 𝑂0,27 𝑁0,25 + 0,9809𝐻2 𝑂 + 0,6306𝐶𝑂2 (33)

Ejercicio 2.4.3

Separación del para-xileno, estireno, tolueno y benceno con un sistema de columnas de

destilación mostrado en la figura 12 extraída del material ‘Cutlip_cap1_cap2’.
Figura 12. Secuencia de separaciones.

a) Calcule las velocidades de flujo molar de 𝑫𝟏 , 𝑫𝟐 , 𝑩𝟏 , 𝑩𝟐 .

Realizamos un balance de materia del rendimiento de cada componente.

Xileno: 0,07𝐷1 + 0,18𝐵1 + 0,15𝐷2 + 0,24𝐵2 = 0,15𝑥70 (34)

Estireno: 0,04𝐷1 + 0,24𝐵1 + 0,10𝐷2 + 0,65𝐵2 = 0,25𝑥70 (35)

Tolueno: 0,54𝐷1 + 0,42𝐵1 + 0,54𝐷2 + 0,10𝐵2 = 0,40𝑥70 (36)

Benceno: 0,35𝐷1 + 0,16𝐵1 + 0,21𝐷2 + 0,01𝐵2 = 0,20𝑥70 (37)

Al realizar un análisis de grados de libertad, observamos que tenemos un sistema de 4
ecuaciones lineales con 4 incógnitas por lo que el problema tiene solución.

Para hallar el valor de las 4 incógnitas utilizamos ‘Matlab’, vamos a denominar ‘Alpha’ a
nuestra matriz 4x4 y denominaremos ‘Beta’ a nuestro vector, procedemos a calcular la
inversa de ‘Alpha’ y la multiplicamos con el vector ‘Beta’.

En la figura 13 podemos observar la matriz ‘Alpha’ de dimensiones 4x4.

Figura 13. Matriz ‘Alpha’ de dimensiones 4x4 realizada en Matlab.

En la figura 14 observamos al vector ‘Beta’ que corresponde al valor de la igualdad de las

ecuaciones de balance de masa de cada componente.

Figura 14. Vector ‘Beta’ realizado en Matlab.

Obtenemos los valores de las velocidades de flujo molar 𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 .

En la figura 15 podemos observar las velocidades de flujo molar 𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 mismas que

𝑘𝑔−𝑚𝑜𝑙
se encuentran en unidades de .
𝑚𝑖𝑛
Figura 15. Velocidades de flujo molar para las corrientes 𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 .

De este modo los valores para las constantes 𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 son:

𝐷1 = 26,25 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐷2 = 8,75 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐵1 = 17,5 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐵2 = 17,5 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

b) Reduzca la velocidad de flujo de alimentación original de la primera columna,

y de forma sucesiva para cada uno de los componentes, primero un 1%, a

continuación un 2%, y calcule las velocidades de flujo correspondiente de

𝑫𝟏 , 𝑫𝟐 , 𝑩𝟏 , 𝑩𝟐 . Explique los resultados.

Primero reducimos un 1% la velocidad de flujo de alimentación original de la primera

columna y realizamos un balance de materia del rendimiento de cada componente.

Xileno: 0,07𝐷1 + 0,18𝐵1 + 0,15𝐷2 + 0,24𝐵2 = (0,15𝑥70) − (70𝑥0,01) (38)

Estireno: 0,04𝐷1 + 0,24𝐵1 + 0,10𝐷2 + 0,65𝐵2 = 0,25𝑥70 − (70𝑥0,01) (39)

Tolueno: 0,54𝐷1 + 0,42𝐵1 + 0,54𝐷2 + 0,10𝐵2 = 0,40𝑥70 − (70𝑥0,01) (40)

Benceno: 0,35𝐷1 + 0,16𝐵1 + 0,21𝐷2 + 0,01𝐵2 = 0,20𝑥70 − (70𝑥0,01) (41)

Para hallar el valor de las 4 incógnitas utilizamos ‘Matlab’, vamos a denominar ‘Alpha2’ a
nuestra matriz 4x4 y denominaremos ‘Beta2’ a nuestro vector, procedemos a calcular la
inversa de ‘Alpha2’ y la multiplicamos con el vector ‘Beta2’.

En la figura 16 podemos observar la matriz ‘Alpha2’ de dimensiones 4x4.

Figura 16. Matriz ‘Alpha2’ de dimensiones 4x4 realizada en Matlab.

En la figura 17 observamos al vector ‘Beta2’ que corresponde al valor de la igualdad de las

ecuaciones de balance de masa de cada componente reducida un 1% del flujo de
alimentación original de la primera columna.

Figura 17. Vector ‘Beta2’ realizado en Matlab.

Figura 18. Velocidades de flujo molar para las corrientes 𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 al reducir un 1% el
flujo de alimentación original de la primera columna.

Realizamos el mismo proceso para obtener los valores de flujo molar para las corrientes
𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 al reducir en un 2% el flujo de alimentación original de la primera columna.

Primero reducimos un 2% la velocidad de flujo de alimentación original de la primera

columna y realizamos un balance de materia del rendimiento de cada componente.

Xileno: 0,07𝐷1 + 0,18𝐵1 + 0,15𝐷2 + 0,24𝐵2 = (0,15𝑥70) − (70𝑥0,02) (42)

Estireno: 0,04𝐷1 + 0,24𝐵1 + 0,10𝐷2 + 0,65𝐵2 = 0,25𝑥70 − (70𝑥0,02) (43)

Tolueno: 0,54𝐷1 + 0,42𝐵1 + 0,54𝐷2 + 0,10𝐵2 = 0,40𝑥70 − (70𝑥0,02) (44)

Benceno: 0,35𝐷1 + 0,16𝐵1 + 0,21𝐷2 + 0,01𝐵2 = 0,20𝑥70 − (70𝑥0,02) (45)

Para hallar el valor de las 4 incógnitas utilizamos ‘Matlab’, vamos a denominar ‘Alpha3’ a
nuestra matriz 4x4 y denominaremos ‘Beta3’ a nuestro vector, procedemos a calcular la
inversa de ‘Alpha3’ y la multiplicamos con el vector ‘Beta3’.

En la figura 19 podemos observar la matriz ‘Alpha3’ de dimensiones 4x4.

Figura 19. Matriz ‘Alpha3’ de dimensiones 4x4 realizada en Matlab.

En la figura 20 observamos al vector ‘Beta’ que corresponde al valor de la igualdad de las

ecuaciones de balance de masa de cada componente reducida un 2% del flujo de
alimentación original de la primera columna.

Figura 20. Vector ‘Beta3’ realizado en Matlab.

En la figura 21 podemos observar las velocidades de flujo molar 𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 cuando se

redujo un 2% el flujo de alimentación original de la primera columna.
Figura 21. Velocidades de flujo molar para las corrientes 𝐷1 , 𝐷2 , 𝐵1 , 𝐵2 al reducir un 2% el
flujo de alimentación original de la primera columna.

Al observar los resultados podemos deducir que al reducir en un 1% el flujo de

alimentación original de la primera columna obtenemos un valor negativo en la corriente de
flujo 𝐷2 lo cual nos indica que ya no habría un flujo de alimentación el cual alimente a la
segunda columna de destilación.

c) Determine las velocidades de flujo molar y las composiciones de B y D para la

parte a).

Para determinar las velocidades de flujo molar para B y D, procedemos a realizar un

balance de masa general para la columna 2 y la columna 3 obteniendo así las siguientes
ecuaciones:

Balance General Columna 2: 𝐷 = 𝐷1 + 𝐵1 (46).

Balance General Columna 3: 𝐵 = 𝐷2 + 𝐵2 (47).

Tenemos como datos conocidos las siguientes constantes 𝐷1 = 26,25 , 𝐷2 = 8,75, 𝐵1 =

17,5, 𝐵2 = 17,5.

Reemplazando las constantes conocidas en las ecuaciones (46) y (47) obtenemos las
velocidades de flujo molar para la corriente B y D

𝐷 = 26,25 + 17,5

𝐷 = 43,75 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛 (48)

 𝐵 = 8,75 + 17,5

𝐵 = 26,25 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛 (49)

Una vez obtenemos la velocidad de flujo molar para la corriente D y B procedemos a hallar
las composiciones de D y B.

Para ello realizamos un balance de cada compuesto en las columnas 2 y 3 respectivamente.

Balance Xileno Columna 2: 𝑊𝐷 = 0,07𝐷1 + 0,18𝐵1 (50)

Balance Estireno Columna 2: 𝑋𝐷 = 0,04𝐷1 + 0,24𝐵1 (51)

Balance Tolueno Columna 2: 𝑌𝐷 = 0,54𝐷1 + 0,42𝐵1 (52)

Balance Benceno Columna 2: 𝑍𝐷 = 0,35𝐷1 + 0,16𝐵1 (53)

Balance Xileno Columna 3: 𝐴𝐵 = 015𝐷2 + 0,24𝐵2 (54)

Balance Estireno Columna 3: 𝐶𝐵 = 0,10𝐷2 + 0,65𝐵2 (55)

Balance Tolueno Columna 3: 𝐸𝐵 = 0,54𝐷2 + 0,10𝐵2 (56)

Balance Benceno Columna 3: 𝐹𝐵 = 0,21𝐷2 + 0,01𝐵2 (57)

Conocemos los valores de las siguientes incógnitas

𝐷1 = 26,25 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐷2 = 8,75 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐵1 = 17,5 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐵2 = 17,5 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐷 = 43,75 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

𝐵 = 26,25 𝐾𝑔 − 𝑚𝑜𝑙/𝑚𝑖𝑛

Procedemos a reemplazar estos valores en las ecuaciones (50) a la (57), despejamos la

incógnita y multiplicamos por 100 para obtener la composición como un porcentaje.

Composición Xileno Corriente D: 𝑊 = 11,4%

Composición Estireno Corriente D: 𝑋 = 12%

Composición Tolueno Corriente D: 𝑌 = 49,2%

Composición Benceno Corriente D: 𝑍 = 27,4%

Composición Xileno Corriente B: 𝐴 = 21%

Composición Estireno Corriente B: 𝐶 = 46,67%

Composición Tolueno Corriente B: 𝐸 = 24,67%

Composición Benceno Corriente B: 𝐹 = 7,67%

Ejercicio 2.9.3

Se propone utilizar un tanque de acero para almacenar dióxido de carbono a 300 K. El

tanque tiene un volumen de 2,5𝑚3 y la presión máxima puede alcanzar, de forma segura,
100atm.

a) Determine el número de moles de 𝑪𝑶𝟐 máximo que se puede almacenar en el

tanque a partir de las ecuaciones de estado mencionadas con anterioridad.

Tenemos como constantes los siguientes valores

𝑇 = 300 𝐾

𝑉𝑡 = 2,5𝑚3

𝑃𝑚𝑎𝑥𝑆 = 100 𝑎𝑡𝑚

𝐿 𝑎𝑡𝑚
𝑅 = 0,08206
𝑔 − 𝑚𝑜𝑙 𝐾

𝑃𝑐 = 79,9 𝑎𝑡𝑚

𝑇𝑐 = 304,2 𝐾

𝜔 = 0,225

Para calcular el número de moles tenemos la siguiente ecuación:

 2,5 ∗ 1000
𝑛𝑚𝑜𝑙𝑒𝑠 = (58)
𝑉

Procedemos a calcular el número de moles máximo que se puede almacenar en el tanque a

partir de las siguientes ecuaciones de estado.

Primero utilizamos la ecuación de Gas Ideal.

𝑃𝑉 = 𝑅𝑇 (59)

Reemplazamos los valores conocidos y de esta manera calculamos el volumen y con ello el
número de moles que se pueden almacenar en el tanque.

Para visualizar todo el proceso de resolución puede referirse al archivo ‘Ejercicio 2.9.3’ en
la parte ‘1.1’, allí se muestra el volumen y el número de moles máxima que se puede
almacenar en el tanque. Obteniendo los siguientes resultados.

Tabla 15. Resultados de volumen y número de moles máxima que se almacena en el tanque
utilizando la ecuación de Gas Ideal.

V f(V) nmoles
0,24618 -1,0303E-13 10155,171

Segundo utilizamos la ecuación de van der Waals.

La ecuación de estado de van der Waals se encuentra al principio de este documento y se la

denominó ecuación (1), la misma tiene dos constantes (a) y (b) cuyas ecuaciones están
escritas como ecuaciones (2) y (3) respectivamente. Primero es necesario hallar el valor de
estas dos constantes, reemplazando los valores dados por el ejercicio, para así hallar el
volumen y número de moles que se puede almacenar en el tanque utilizando la ecuación de
estado de van der Waals. Para visualizar la resolución de literal puede referirse al
documento Excel denominado ‘Ejercicio 2.9.3’ en la parte ‘1.2’ obteniendo, así como
resultado los siguientes datos:
Tabla 16. Resultados de volumen y número de moles máxima que se almacena en el tanque
utilizando la ecuación de estado de van der Waals.

a b V f(V) nmoles
3,29017089 0,03905296 0,07497742 1,36345E-05 33343,3704

Tercero utilizamos la ecuación de Soave-Redlich-Kwong.

Para la ecuación de Soave-Redlich-Kwong tenemos las siguientes ecuaciones

𝑅𝑇 𝛼𝑎
𝑃= −[ ] (60)
𝑉−𝑏 𝑉 (𝑉 + 𝑏 )

Donde

𝑅2 𝑇𝑐 2
𝑎 = 0,42747 ( ) (61)
𝑃𝑐

𝑅𝑇𝑐
𝑏 = 0,088664 ( ) (62)
𝑃𝑐

2
𝑇
𝛼 = [1 + 𝑚 (1 − √ )] (63)
𝑇𝑐

𝑚 = 0,48508 + 1,5517𝜔 − 0,1561𝜔2 (64)

Reemplazamos los valores para hallar las constantes de la ecuación, a continuación,

reemplazamos en la ecuación (60) y hallamos los valores de volumen y número de moles
máximo que se almacena en el tanque, para visualizar la resolución de este literal puede
referirse al archivo de Excel ‘Ejercicio 2.9.3’ en la parte ‘1.3’, de esta manera obtenemos
los siguientes resultados.

Tabla 17. Resultados de volumen y número de moles máxima que se almacena en el tanque
utilizando la ecuación de estado de Soave-Redlich-Kwong.

A b ɑ m V f(V) nmoles
3,33380587 0,02706839 0,0007327 0,826312188 0,27332241 0,000306458 9146,70715
 Cuarto utilizamos la ecuación de Peng-Robinson.

Para la ecuación de Peng-Robinson tenemos las siguientes ecuaciones.

𝑅𝑇 𝑎(𝑇)
𝑃= −[ ] (65)
𝑉−𝑏 𝑉(𝑉 + 𝑏) + 𝑏(𝑉 − 𝑏)

Donde

𝑅𝑇𝑐
𝑏 = 0,07780 ( ) (67)
𝑃𝑐

𝑅2 𝑇𝑐 2
𝑎(𝑇) = 0,45724 ( ) 𝛼(𝑇) (66)
𝑃𝑐

𝛼(𝑇) = [1 + 𝑘 (1 − 𝑇/𝑇𝑐 )0,5 ]2 (68)

𝑘 = 0,37464 + 1,5422𝜔 − 0,26992𝜔2 (69)

Reemplazamos los valores dados por el ejercicio para hallar las constantes de la ecuación, a
continuación, reemplazamos en la ecuación (65) y hallamos los valores de volumen y
número de moles máximo que se almacena en el tanque, para visualizar la resolución de
este literal puede referirse al archivo de Excel ‘Ejercicio 2.9.3’ en la parte ‘1.4’, de esta
manera obtenemos los siguientes resultados.

Tabla 18. Resultados de volumen y número de moles máxima que se almacena en el tanque
utilizando la ecuación de estado de Peng-Robinson.

b a(T) ɑ(T) k V f(V) nmoles

0,02430656 3,60104393 1,00983295 0,7079838 0,3582568 4,50343E-05 6978,23468

Quinto utilizamos la ecuación de Beattle-Bridgeman.

Para la ecuación de Beattle-Bridgeman tenemos las siguientes ecuaciones y constantes.

𝑅𝑇 𝛽 𝛾 𝛿
𝑃= + 2+ 3+ 4 (70)
𝑉 𝑉 𝑉 𝑉

Donde
 𝑅𝑐
𝛽 = 𝑅𝑇𝐵0 − 𝐴0 − (71)
𝑇2

𝑅𝑐𝐵0
𝛾 = 𝑅𝑇𝐵0 𝑏 + 𝐴0 𝑎 − (72)
𝑇2

𝑅𝐵0 𝑏𝑐
𝛿= (73)
𝑇2

Y las constantes proporcionadas por el problema son

𝐴0 = 5,0065

𝑎 = 0,07132

𝐵0 = 0,10476

𝑏 = 0,07235

𝑐 = 66,0𝑥104

Reemplazamos los valores dados por el ejercicio para hallar las constantes de la ecuación, a
continuación, reemplazamos en la ecuación (70) y hallamos los valores de volumen y
número de moles máximo que se almacena en el tanque, para visualizar la resolución de
este literal puede referirse al archivo de Excel ‘Ejercicio 2.9.3’ en la parte ‘1.5’, de esta
manera obtenemos los siguientes resultados.

Tabla 19. Resultados de volumen y número de moles máxima que se almacena en el tanque
utilizando la ecuación de estado de Beattle-Bridgeman.

Beta Gamma Delta V f(V) nmoles

-3,02929165 0,48061113 0,00456107 0,214162128 -1,8014E-05 11673,39914

b) Suponiendo que la ecuación de Beattle-Bridgeman es la más apropiada, ¿qué

porcentaje de error se comete en el cálculo del número de moles si se utilizan
otras correlaciones?

Establecemos que el número de moles obtenido con la ecuación de Beattle-Bridgeman es

nuestro valor teórico, y las ecuaciones de Gas Ideal, van der Waals, Soave-Redlich-Kwong,
Peng-Robinson las asumimos como nuestro valor experimental.
La ecuación para calcular el porcentaje de error es

|𝑉𝐸𝑥𝑝 − 𝑉𝑇𝑒𝑜 |
𝑒% = 𝑥100 (74)
𝑉𝑇𝑒𝑜

Procedemos a calcular el error entre la ecuación de Beattle-Bridgeman y de Gas Ideal.

|10155,171 − 11673,39914|
𝑒% = 𝑥100
11673,39914

𝑒% = 13,01%

Procedemos a calcular el error entre la ecuación de Beattle-Bridgeman y de Gas Ideal.

|10155,171 − 11673,39914|
𝑒% = 𝑥100
11673,39914

𝑒% = 13,01%

Procedemos a calcular el error entre la ecuación de Beattle-Bridgeman y de van der Waals.

|33343,3704 − 11673,39914|
𝑒% = 𝑥100
11673,39914

𝑒% = 185,63%

Procedemos a calcular el error entre la ecuación de Beattle-Bridgeman y de Soave-Redlich-

Kwong.

|9146,70715 − 11673,39914|
𝑒% = 𝑥100
11673,39914

𝑒% = 21,64%

Procedemos a calcular el error entre la ecuación de Beattle-Bridgeman y de Peng-

Robinson.

|6978,23468 − 11673,39914|
𝑒% = 𝑥100
11673,39914

𝑒% = 40,22%
Para revisar los cálculos realizados en este ejercicio puede referirse al archivo Excel
denominado ‘Ejercicio 2.9.3’, en el mismo se comprueba los valores, las fórmulas y el
procedimiento seguido para resolver cada literal.

Ejercicio 2.10.3

a) Calcule la temperatura del punto de burbujeo y la composición del equilibrio

asociado con una mezcla líquida del 10% en moles de n-pentano y el 90% en
moles de n-hexano a 1 atm.
b) Repita los cálculos para mezclas de líquidos que contengan desde un 0% hasta
un 100% en moles de n-pentano.

Para resolver los literales a) y b) necesitamos las siguientes ecuaciones.

1064,63
𝑙𝑜𝑔𝑃𝐴 = 6,85221 − (74)
𝑇 + 232,0

1171,53
𝑙𝑜𝑔𝑃𝐵 = 6,87776 − (75)
𝑇 + 224,336

𝑓 (𝑇𝑏𝑝) = 𝑥𝐴 𝑃𝐴 + 𝑥𝐵 𝑃𝐵 − 760 (76)

𝑥𝐴 𝑃𝐴
𝑦𝐴 = (77)
760

𝑥𝐵 𝑃𝐵
𝑦𝐵 = (78)
760

En el archivo Excel denominado ‘Ejercicio 2.10.3’ se puede observar el procedimiento

realizado para resolver los 3 literales que plantea este ejercicio, para el literal a) primero
debemos hallar los valores para las constantes planteadas por el ejercicio y procedemos a
resolver, Para el literal b) procedemos a realizar el mismo procedimiento del literal a) para
los nuevos valores establecidos en la orden del literal b). Para el literal c) procedemos a
graficar la T vs XA y T vs XB de acuerdo a los valores obtenidos en el literal b).

Para visualizar todos los cálculos realizados y los valores obtenidos puede referirse al
archivo ‘Ejercicio 2.10.3’ en el que se encuentra todo el procedimiento detallado y las
fórmulas detalladas.