Ec De Estado De SRK

ECUACION DE ESTADO DE SOAVE-REDLICH-KWONG:

La ecuación de Soave-Redlich-Kwong fue la primera 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 a la que se le incluye el parámetro ω, el cual es denominado factor acéntrico y es característico de la estructura molecular.

P= RT/(V ̂-b)- (∝a)/(V ̂(V ̂+b))

Dónde:

R=constante de los gases (8.31451 J⁄(K.mol))

a=0.42747 〖〖(RT〗_c)〗^2/Pc

b=0.08664 (RT_c)/P_c

m=0.48508+1.55171ω -0.15613〖ω 〗^2

∝=[1+m(1-√Tr)] 2

T_r=T/T_c

Y donde ω es el factor céntrico del compuesto.

*Para el hidrogeno:

∝ =1.202exp⁡(-0.30288T_r)

En 1972 Soave reemplazó el término a/√(T) de la ecuación de Redlich-kwong por una expresión ∝(T,ω) función de la temperatura y el factor céntrico.

La función α fue concebida para cuadrar con los datos de las presiones de vapor de los hidrocarburos; esta ecuación describe acertadamente el comportamiento de estas sustancias.

EJEMPLO:

Un tanque se carga con 10 libras de gas propano a 25 ° C. El tanque es de 24 pulgadas de alto y 12,5 pulgadas de diámetro. Calcular la presión del tanque utilizando la ecuación de estado de Soave-Redlich-Kwong.

Formulas a usar:

Tenemos que:

P= RT/(V ̂-b)- (∝a)/(V ̂(V ̂+b))

Donde:

R=constante de los gases (8.31451 J⁄(K.mol))

a=0.42747 〖〖(RT〗_c)〗^2/Pc

b=0.08664 (RT_c)/P_c

m=0.48508+1.55171ω -0.15613〖ω 〗^2

∝=[1+m(1-√Tr)] 2

T_r=T/T_c

DATOS

RESOLVIENDO

Pc = 616psia… 41.90 atm… 4.60x106 Pa

Tc = 206.06°F… 369.85K

ω = 0.152

D = 12.5in…

T = 25°C.. 298.15K

R = 6.25in

h = 24in

Masa = 10lb

MM = 44kg/kmol

V ̂ = ?

R= 8.314 (m3.Pa/ mol.K)

Vtan= πr^2 h = π(0.15875m)^2 (0.60960m) = 0.048264m3

η_C3H8= m/MM=4.535kg/(44kg/kgmol)=.1028kgmol

= 102.8383gmol

V ̂= V_gas/n_gas = (0.048264m^3)/102.8383gmol=0.000469 m^3/gmol

T_r=T/T_c = 298.15K/(3690.85 K)=0.806137623

a=0.42747[( 〖((8.314.(m^3.Pa)/(mol.K))(369.85))〗^2)/(〖4.26x10〗^6 Pa)]=〖9.50x10〗^(-01) m^6 Pa/〖mol〗^2

b=0.08664⌈((8.314m^3.pa/mol K)(369 K))/(〖4.26x10〗^6 Pa)⌉=〖6.26x10〗^(-5) m^3/mol

m=0.48508+1.55171(0.152)-0.15613(0.152)^2=0.71733339

α= [1+0.71733339(1-√0.806137623)]2 =1.151917

...