UNIVERSIDAD NACIONAL EXPERIMENTAL POLITÉCNICA

“ANTONIO JOSÉ DE SUCRE”

VICERECTORADO BARQUISIMETO
DEPARTAMENT
O DE INGENIERÍA QUIMICA
SECCIÓN DE INGENIERIA

ALUMNOS:
​ARIANNA SANCHEZ C.I: 26.561.909
DANIELA SANCHEZ C.I: 25.340.450
MAURICE SEIJAS C.I: 24.354.742
ASIGNATURA: ​Principios de
Ingenieria Química
DOCENTE:​ Eduard Rodriguez
SECCIÓN:​ 01
LAPSO:​ 2018-2

BARQUISIMETO, JUNIO 2019

 ECUACIÓN DE VAN DER WAALS: Johannes D. van der Waals en 1873
propuso una modificación de la ley de gas ideal, para tener en cuenta el tamaño molecular y
las fuerzas de interacción moleculares. Se la refiere normalmente como la ecuación de
estado de van der Waals.

La ecuación de Van der Waals permite describir de forma más satisfactoria el

comportamiento termodinámico de muchos gases, no solo en altas temperaturas y bajas
presiones, pero en un dominio más extenso de temperaturas y presiones nominadamente
cuando el gas está próximo de la condensación

Por encima de la temperatura crítica la ecuación de Van der Waals es una mejora de la
ley del gas ideal, y para temperaturas más bajas la ecuación es también cualitativamente
razonable para el estado líquido y estado gaseoso a baja presión. Sin embargo, el modelo
Van der Waals no es adecuado para los cálculos cuantitativos rigurosos, útil restante sólo
con fines educativos y de calidad.

En la transición de fase de primer orden, el rango de (P, V, T), donde la fase líquida y
la fase gaseosa se encuentran en equilibrio, no lo exhibe el hecho empírico de que p es
constante en función de V a una temperatura dada aunque este comportamiento se pueda
insertar fácilmente en el modelo la Van der Waals, el resultado ya no es un modelo
analítico simple, y otros (como los basados en el principio de estados correspondientes)
logran un mejor ajuste con más o menos el mismo trabajo.

Ecuación

(P + a( Vn )2 ) (V − n * b) = n * R * T

Donde:

P es la presión del fluido,

V es el volumen total del recipiente en que se encuentra el fluido,

A mide la atracción entre las partículas

B es el volumen disponible de un mol de partículas

N es el número de moles

R es la constante universal de los gases ideales,

T es la temperatura, en kelvin.
2 2
27 R *T c
a= 64 *( Pc
)

1 R*T c
b= 8 *( Pc
)

Así mismo está también puede representarse por la siguiente ecuación, empleando el
volumen molar ( V̂ ):

(P + ) (Vˆ − b) = R * T
a
V̂

ECUACIÓN DE PENG- ROBINSON​: Fue desarrollada en la década del 70 por Ding-Yu

Peng y Donald Baker Robinson. Ding-Yu Peng era alumno de doctorado del Prof. Donald
Robinson en la Universidad de Calgary, Canadá. Este trabajo investigativo fue apadrinado
por la Junta Canadiense de Energía (Canadian Energy Board) para solucionar
específicamente sistemas de gas natural. Finalmente en 1976 se publicó este trabajo.

La ecuación Peng-Robinson se diseñó para satisfacer las siguientes metas:

1. Los parámetros se expresen en términos de propiedades críticas y el factor acéntrico.

2. El modelo debe proveer una exactitud razonable cerca del punto crítico, particularmente
para los cálculos del factor de compresibilidad y densidad de líquido.

3. Las reglas que se mezclan no debe emplear más que los parámetros de interacción
binarios, los cuales deben ser independiente de la temperatura, presión y composición.

4. La ecuación debe ser aplicable a todos los cálculos de propiedades de fluidos en procesos
del gas natural.
Gran parte de la ecuación Peng-Robinson muestra un desempeño similar a la ecuación
Soave, aunque es generalmente superior en la predicción de las regiones críticas de fase y
de densidades de líquido de cualquier material, especialmente los no polares por ejemplo en
las moléculas carbonadas por lo que es muy aplicada en la industria petrolera del gas
natural. La ecuación Peng-Robinson se expresa en los siguientes términos:

R*T a*α
P = Vˆ −b
− Vˆ (Vˆ +b)+b(Vˆ −b)

P = Presión (atmósferas)

V̂ = V/n = Volumen molar, el volumen de un mol de gas

T = Temperatura (K)

R = constante de los gases ideales

Tc= temperatura crítica

Pc= Presión Crítica

W= Factor Acéntrico
2 2
a = 0, 45724 * ( R P*Tc c )

b = 0, 07780 * ( RP*Tc c )

k = 0, 37464 + 1, 54226 * (w) − 0, 26992* w2

2
[
α = 1 + k (1 − √T r ) ]

Existen diferentes casos donde se puede aplicar esta ecuación como por ejemplo:

Calculo de la fugacidad y coeficiente de fugacidad mediante la ecuación Peng- Robinson

ECUACIÓN DE REDLICH –KWONG​: Es una ecuación algebraica empírica que
relaciona temperatura, presión, y volumen de los gases. Es generalmente más precisa que la
ecuación de Van der Waals y los Gases Ideales a temperaturas arriba de la temperatura
crítica. Fue formulada por Otto Redlich y Joseph Neng Shun Kwong en 1949.12​ Demostró
que una ecuación de estado cúbica de dos parámetros bien podría reflejar la realidad en
muchas situaciones, junto a un mucho más complicado modelo de Beattie–Bridgeman y la
ecuación de Benedict-Webb-Rubin que se utilizaba en el momento. La ecuación
Redlich-Kwong ha sufrido muchas revisiones y modificaciones, a fin de mejorar su
precisión en términos de predecir las propiedades de más compuestos en fase gas, tanto
como mejorar simulaciones de condiciones a temperaturas bajas, incluyendo el equilibrio
vapor-liquido.

La ecuación de Redlich-Kwong es adecuada para calcular las propiedades de la fase

gaseosa cuando el coeficiente de la presión y la presión crítica es menor que la mitad del
cociente de la temperatura y temperatura reducida

la ecuación de estado de redlich-kwong es muy utilizada para calcular propiedades

termodinámicas para la fase de vapor en combinación con modelos termodinámicos más
complejos.

R*T a
P = Vˆ −b
− √T *Vˆ *(Vˆ +b)

2 2,5
a = 0, 42748 * ( R *PTcc )

b = 0, 08664 * ( RP*Tc c )

ECUACIÓN DE SOAVE-REDLICH –KWONG​: Soave, en 1972 propuso una

modificación de la forma simple de la ecuación de Redlich-Kwong donde el parámetro a
fue hecho dependiente de la temperatura de modo que la curva de presión de vapor pueda
ser reproducida correctamente. La ecuación de estado requiere del ingreso de tres
parámetros por compuesto puro: Tc, Pc y ω
El uso de la ecuación de estado Soave-Redlich-Kwong tiene diversas aplicaciones en el
área de hidrocarburos líquidos gracias al cálculo de su volumen especifico,estas son: Los
parámetros se pueden expresar en términos de propiedades críticas y el factor acéntrico. El
modelo debe proveer una exactitud razonable cerca del punto crítico, particularmente para
los cálculos del factor de compresibilidad y densidad del líquido. La ecuación de
Soave-Redlich-Kwongpermite describir de forma satisfactoria el comportamiento
termodinámico del Estireno, ya que este es un hidrocarburo aromático, liquido e incoloro.
En general sirve para la estimación del factor acéntrico, la temperatura crítica y la presión
crítica usando la ecuación Soave-Redlich-Kwong. También con su ecuación cubica es
posible hallar fácilmente la presión para una sustancia o fluido.

R*T a*α
P = Vˆ −b
− Vˆ (Vˆ +b)

2 2
a = 0, 42748 * ( R P*Tc c )

b = 0, 08664 * ( RP*Tc c )

k = 0, 48508 + 1, 55171 * (w) − 0, 1561* w2

2
[
α = 1 + k (1 − √T r ) ]
FACTOR DE COMPRESIBILIDAD (Z): z es lo que se denomina factor de
compresibilidad, que representa cuán alejado estamos del modelo ideal de los gases. Si z
vale 1, entonces el modelo ideal es aplicable. Dependiendo de la presión los valores de z
pueden ser menores o mayores a 1. Si la presión es baja dominan las fuerzas
intermoleculares atractivas, las cuales reducen el volumen respecto al del gas ideal y z será
menor que 1. A presiones altas dominan las fuerzas repulsivas, las cuales expandirán el gas
respecto al gas ideal y z será mayor que 1. La desviación z se puede calcular con cualquiera
de los modelos matemáticos

P * V̂ = z * R * T
 Ejercicios Aplicando Las ecuaciones de Estado

1) Aplicando la ecuación de Peng-Robinson: Se desea almacenar 30 Kg de SO​2 a 50ºC

en un tanque cuya capacidad es de 250 L se desea saber cuál es la presión dentro del
tanque:

Sol:

Datos

m= 30 Kg SO​2 ​ T​c​= 430,7 K

T= 50 ºC = 323 K P​c​= 77,8 atm

V= 250 L P=? w=0,251

R*T a*α
P = Vˆ −b
− Vˆ (Vˆ +b)+b(Vˆ −b)
(I)

Ahora;

a = 0, 45724 * ( R2 *T c 2
Pc ) (II)
b = 0, 07780 * ( ) (III)
R*T c
Pc

k = 0, 37464 + 1, 54226 * (w) − 0, 26992* w2 (IV)

2
[
α = 1 + k (1 − √T r ) ] (V)

Sustituyendo cada valor, en (II), (III) se tiene que:

a = 0, 45724 * ( 0,0822 *430,72

77,8 ) (II)
atm*L2
a = 0, 20756 mol2

b = 0, 07780 * ( 0,082*430,7
77,8 )
 L
b = 0, 035317 mol

Luego:

T
Tr = Tc

323 K
Tr = 430,7 K = 0, 74994

Resolviendo para α

2
[
α = 1 + (0, 37464 + 1, 54226 * (w) − 0, 26992* w2 )(1 − √T r ) ]
2
[
α = 1 + (0, 37464 + 1, 54226 * (0, 251) − 0, 26992* 0, 2512 )(1 − √0, 74994) ]
α = 1, 2096
V 250 L
V̂ = n
= 468,24 mol
= 0, 53391 L/mol

Sustituyendo todos los valores en (I)

0,082*323 K 0,20756*1,2096
P = 0,53391−0,035317
− 0,53391 (0,53391 +0,035317)+0,035317(0,53391 −0,035317)

P = 52, 340 atm

Por lo que la presión del gas dentro del tanque es de 52,340 atm

2) Se tiene un tanque de almacenamiento el cual contiene gas de O​2​, donde el mismo

se encuentra sometido a una presión de 150 atm y a una temperatura de 4º C, se
sabe también que su factor acéntrico (w) es de 0,021. A) ¿Cuál es su volumen
molar? B) Si el tamaño del tanque es de 1000 m​3​, cuantas toneladas del gas de O​2 se
encuentran dentro? Resolver usando las ecuaciones de Estado.

Sol:

Datos

T​c​= 154,4 K P​c​= 49,7 atm

P= 150 atm T= 277 K

W= 0,021 V̂ =?

V= 1000 m​3​ ton O​2​=?

Primeramente se utilizará el factor de comprensibilidad (Z) para luego comparar con cada
ecuación de Estado, para ello usaremos su ecuación:

P * V̂ = z * R * T

Primero se determinan la Presión reducida y la temperatura reducida

T 277 K
Tr = Tc = 154,4 = 1, 80

P 150 atm
Pr = Pc = 49,7 atm = 3, 01

Usando las graficas para encontrar el factor de comprensibilidad z, se determinó que este
equivalía aproximadamente a 0,92. Ahora bien, despejando el volumen molar de la
ecuación se obtiene que:

z *R*T
V̂ = P

0,92*0,082*277
V̂ = 150

L
V̂ = 0, 1393 mol

a) Utilizando la ecuación de Van der Waals

(P + ) (Vˆ − b) = R * T
a
V̂

Para ello es necesario calcular los valores a y b, para resolver el problema

a= 27
64 * ( 0,0822 *154,42
49,7 ) = 1, 3607
b= 1
8 * ( 0,082*0,03184
49,7 ) = 0,03184
(150 + 1,3607
V̂ ) (Vˆ − 0, 03184) = 0, 082 * 277
Claramente no es sencillo despejar el volumen molar, por lo que se empleo un CAS para
ello, dicho sistema emplea el método de Newton, para resolver la ecuación.

L
V̂ = 0, 1308 mol

b) Utilizando SRK

R*T a*α
P = Vˆ −b
− Vˆ (Vˆ +b)

a = 0, 42748 * ( 0,0822 *154,42

49,7 ) = 1, 3787
b = 0, 08664 * ( 0,082*154,4
49,7 ) = 0, 02202
k = 0, 48508 + 1, 55171 * (0, 021) − 0, 1561* 0, 0212 = 0,5176
2
[
α = 1 + 0, 5176(1 − √1, 8) ] = 0, 6776

0,082*277 1,3787*0,6776
150 = −
Vˆ −0,02202 Vˆ (Vˆ +0,02202)

Nuevamente empleando un CAS

L
V̂ = 0, 1412 mol

c) Utilizando Peng-Robinson
R*T a*α
d) P = Vˆ −b
− Vˆ (Vˆ +b)+b(Vˆ −b)

a = 0, 45724 * ( 0,0822 *154,42

49,7 ) = 1, 4747
b = 0, 07780 * ( 0,082*154,4
49,7 ) = 0, 01982
k = 0, 37464 + 1, 54226 * (0, 021) − 0, 26992* 0, 0212

k = 0, 40690
2
[
α = 1 + 0, 40690(1 − √1, 8) ] = 0, 7413

0,082*277 1,4747*0,7413
150 = −
Vˆ −0,01982 Vˆ (Vˆ +0,01982)+0,01982(Vˆ −0,01982)

Aplicando un CAS para el V̂

 L
V̂ = 0, 1349 mol

e) Utilizando la ecuación de Reldich-Kwong

R*T a
P = Vˆ −b
− √T *Vˆ *(Vˆ +b)

a = 0, 42748 * ( 0,0822 *154,42,5

49,7 ) = 17, 131
b = 0, 08664 * ( 0,082*154,4
49,7 ) = 0,02207
0,082*277 17,131
150 = − √277*Vˆ *(Vˆ +0,02207)
Vˆ −0,02207

Empleando un CAS para resolver la ecuación

L
V̂ = 0, 1373 mol

Apartado B del ejercicio para resolver dicha parte se va a emplear el promedio de los
volúmenes molares determinados anteriormente para indicar la cantidad de toneladas
almacenadas

L
V̂ = 0, 1367 mol

⇨ Ton O​2​= 1*10^6L*1/0,1367(L/mol)*32 ton O2/1*10^6 mol O2 = 234,08 Ton O​2