>	restart;#"m06_p39"

Se desea ajustar la curva de presión de vapor del n-butano en la forma ln(p/pb)=A*(1-Tb/T). Se pide:

a) Estimar el valor de A a partir del ajuste con el modelo de estados correspondientes ln(pR)=K*(1-1/TR) con K=6 y Tb=0.6Tcr.

b) Estimar el valor de A a partir de la ecuación de Clapeyron integrada en el entorno del punto de ebullición.

Datos:

>	read"../therm_eq.m":read"../therm_proc.m":with(therm_proc):

>	su:="C4H10":dat:=[TRb=0.6,K=6];

Esquema:

>

>

Eqs. const.:

>	eqETg:=subs(eq1_11,eq1_12);eqEE:=eq1_16;gdat:=get_gas_data(su):ldat:=get_liq_data(su):dat:=op(dat),gdat,ldat,Const,SI2,SI1:get_pv_data(su):

a) Estimar el valor de A a partir del ajuste con el modelo de estados correspondientes ln(pR)=K*(1-1/TR) con K=6 y Tb=0.6Tcr,

>	TbMEC:=TRb*T[cr];TbMEC_:=subs(dat,TbMEC):'TbMEC'=evalf(%,2);pRMEC:=exp(K*(1-1/TRb));pRMEC_:=subs(dat,pRMEC):'pRMEC'=evalf(%,2);pbMEC_:=evalf(subs(dat,pRMEC_*p[cr])):'pbMEC'=evalf(%,2);A:=K/TRb;A_:=subs(dat,A):'A'=evalf(%,2);

b) Estimar el valor de A a partir de la ecuación de Clapeyron integrada en el entorno del punto de ebullición.

>	eqClap:=ln(p/p0)=-(h[lv0]/R)*(1/T-1/T[b]);A:=h[lv0]/(R*T[b]);A_:=subs(dat,A):'A'=evalf(%,2);

>	lnp_Clap:=ln(p0)+A*(1-T[b]/T);lnp_MEC:=ln(p[cr])+K*(1-T[cr]/T);plot(subs(dat,SI0,[[[T[b],0],[T[b],ln(p0)]],lnp_MEC,lnp_Clap]),T=100..400);