Ejercicio 2
Ejercicio 2
Ejercicio 2
27( R2 T 2 c) RTc
a= b=
64( Pc) 8 Pc
RT a
P= − 2
v−b b
Donde:
v: Volumen específico (m3/kgmol)
R: 8.314 kJ/(kgmol*K)
Esta ecuación de estado se puede reescribir utilizando una expresión
polinómica de grado 3, como la siguiente:
(
v3 − b+
P )
RT 2 a ab
v + v− =0
P P
Pv
Z=
RT
Pc = 4599KPa Tc = 190,6K
Solución:
Evaluaremos el valor del factor de compresibilidad. Para ello, utilizaremos el
método de Newton-Raphson siguiendo los pasos descritos con anterioridad.
Antes de iniciar, debemos tener en cuenta que una ecuación cúbica presenta
tres raíces, una de estas raíces representa al volumen del gas y no del líquido,
y esta es la raíz que nos interesa.
Escribimos nuestra función:
f ( v )=a 3 v 3+ a2 v 2 +a1 v +a 0
(
f ( v )=v 3− b+
P )
RT 2 a
v + v−
P
ab
P
' 2
f ( v )=3 v −2 b+ ( RT
P)v+
a
P
' 2
f ( v )=3 v −2,58034 v +0,11517(2)
Como se requiere calcular el volumen específico del gas, utilizaremos como
estimado inicial la aproximación del gas ideal, el cual es Z=1, sería el inicio más
lógico para la iteración, entonces la ecuación general de los gases quedaría de
la siguiente manera:
RT 3
v= =1,2471m /kmol
P
Ahora, podemos comenzar la iteración, fijamos la tolerancia en 0.001 y
seleccionamos nuestro valor de inicio el cual es:
3
v 0=1,2471 m / kmol
Sustituimos este valor en las ecuaciones (1) y (2), obteniendo los siguientes
resultados:
f ( v 0 )=0,07169
f ' ( v 0 )=1,5630
Seguidamente, evaluamos la solución estimada de la siguiente manera:
f ( v 0)
v1 =v 0− '
=1,20123
f ( v0)