diagrama de ponchon

Objetivo:
Determinar las propiedades termodinámicas de la mezcla hexano-octano mediante ecuaciones de estado y un modelo de actividad
Introducción:
Las propiedades termodinámicas.
Son características que se pueden observar, medir o cuantificar en las sustancias o en los sistemas.
La cantidad y tipo de propiedades que se puedan establecer para un sistema dependen del tipo de observación quese halla establecido para el análisis del sistema. Por ejemplo si el enfoque usado es el macroscópico se pueden establecer propiedades como temperatura (T), presión (P), energía (e), energía interna (u), y entalpia (h) y otras, que de ningún modo serían establecidas utilizando el enfoque microscópico.
Para apropiar mejor esta definición, a continuación se hace detalle sobre dos aspectosprimordiales de las propiedades.
Las propiedades termodinámicas cuyos valores dependen del tamaño del sistema son llamadas propiedades extensivas (v.g. volumen, entropía, energía interna). Las propiedades denominadas propiedades intensivas son aquellas que no dependen del tamaño del sistema (v. g. presión, temperatura).
Es claro entonces que las propiedades extensivas son aditivas, así, si el sistemase divide en un número de subsistemas, el valor de la propiedad para el sistema entero es igual a la suma de los valores de los subsistemas. El volumen es pues una propiedad extensiva.
Por otro lado las propiedades intensivas no dependen de la cantidad presente de materia y éstas no pueden ser obtenidas como la suma de todos los subsistemas, como es el caso de la temperatura.
Desarrollo:
Parapoder encontrar las propiedades termodinámicas de mi mezcla, primero fue necesario saber qué modelo de actividad y que ecuación tenia que utilizar, para esto utilice un tutorial llamado enciclopedia, que mediante las propiedades que tiene cada compuesto de mi mezcla me dice cual es la ecuación que se ajusta, mi mezcla es HEXANO-OCTANO, y después de analizarlo en la enciclopedia me dijo que debíautilizar PR y WILSON.
Después de acuerdo a la información obtenida anteriormente, sabia que de mi mezcla el compuesto más volátil es el hexano, por lo que mediante otro programa llamado HYSYS obtuve mis propiedades críticas, parámetros de Wilson y factor acéntrico
Propiedades Iniciales:
T=350 C, P=101.325 kPa
PROPIEDAD
HEXANO
OCTANO
Temperatura Crítica (Tc)
234.7 C
295.4 C
PresiónCrítica (Pc)
3032 kPa
2497 kPa
Factor Acéntrico (W)
0.3007
0.4018
Volumen (VL)
0.13115
0.16316

Parámetros de Wilson
12 –λ11=-79.836
21 –λ22=80.236
Kij = q12 =0.00107
Y la R obtenida en el HYSYS es 1.98721 cal/gmolK
Primero en Fortan obtuve las derivadas, para poder obtener las propiedades termodinámicas

USE MSIMSLMS
INTEGER KORDER, NCOUNT, NOUT
REALBGSTEP, DERV, FCN, TOL, X
EXTERNAL FCN
C Get output unit number
X = 350.0
BGSTEP = 2
TOL = 0.01
KORDER = 1
NCOUNT = 1
DERV = DERIV(FCN,KORDER,X,BGSTEP,TOL)
WRITE (NOUT,99999) DERV
99999 FORMAT (/, 1X, 'First derivative of FCN is ', 1PE10.3)
END
C
REAL FUNCTION FCN (X)
REALX,P,Tr1,Tr2,Pc1,Pc2,w1,w2,Tc1,Tc2,R,Y1,Y2,a1,a2,a11,a22
REAL q12,alfa1,alfa2,k1,k2

C FCN = a, ésta a es la de la regla de mezclado. La de parcial de a respecto a T
C datos usados por el profe. En mi caso, lo que cambiaré para cada punto ex Y1 y Y2
P=101.325
R=8.31434
Y1=0.2
Y2=0.8
q12=0.00107
c q12 es la k12

c propiedades criticas
Tc1=234.7+273.15Tc2=295.4+273.15
Pc1=3032
Pc2=2497
w1=0.3007
w2=0.4018
c calculo de parametros de componentes puros
a1=0.45724*(R**2)*(Tc1**2)/Pc1
a2=0.45724*(R**2)*(Tc2**2)/Pc2
k1=0.37964+1.54226*w2-0.26992*w2**2
k2=0.37964+1.54226*w2-0.26992*w2**2
Tr1=X/Tc1
Tr2=X/Tc2
alfa1=(1+k1*(1-Tr1**0.5))**2
alfa2=(1+k2*(1-Tr2**0.5))**2

c parametros cruzados
a11=a1*alfa1
a22=a2*alfa2...