Design og modellering af metanolanlæg til VEnzin-visionen Bilag

methanol.for
d:/DTU/Eksamensprojekt/bilag/
$ sigma,vv)
call sub_alpha_id(alpha_id,tau,omega,f_i,g_i)
s=R*(tau*(dalpha_id_dtau+dalpha_dtau_omega)−alpha−log(omega/
$ omega_a)−alpha_id)
endif
return
end
subroutine calculate_x(x,tau,omega,omega_g,omega_l)
integer IMAX
double precision x,tau,omega,T,rho,rho_l,rho_g,omega_g,omega_l,Z
$ ,P_g,EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C
logical fiters,iterate_rho_l
common/constants/EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C,IMAX
$ ,fiters
common/switches/iterate_rho_l
T=T_ref/tau
rho=omega*rho_ref
if (T.lt.T_c) then
call rho_l__rho_g(tau,omega_l,omega_g)
rho_l=rho_ref*omega_l
rho_g=rho_ref*omega_g
if ((rho.gt.rho_l*.99).and.(rho.lt.rho_l*1.01).and
$ .iterate_rho_l)then
call calculate_Z(Z,tau,omega_g)
P_g=Z*R*T/(1/rho_g*1.d−6)
call iterer_omega_l(P_g,tau,omega_l)
rho_l=omega_l*rho_ref
endif
if ((rho.gt.rho_g).and.(rho.lt.rho_l)) then
x=rho_g*((rho_l/rho)−1)/(rho_l−rho_g)
elseif (rho.lt.rho_g) then
x=1.0+rho_g−rho
elseif (rho.gt.rho_l) then
x=rho_l−rho*10
else
x=rho_g*((rho_l/rho)−1)/(rho_l−rho_g)
endif
if (x.eq.0) x=1.d−12
else
x=10*T/T_c+1
endif
return
end
subroutine calculate_u(u,tau,omega)
integer r_i(44),s_i(44),k_i(44),j_i(44),b_i(44),IMAX
double precision u,tau,omega,T,dalpha_id_dtau,N_i(44),omega_a,
$ c_i(44),f_i(10),g_i(10),sigma,vv,dalpha_dtau_omega,x,omega_g
$ ,omega_l,u_g,u_l,EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C
logical fiters,iterate_rho_l
common/constants/EPSV,R,P_c,T_c,T_ref,rho_ref,rho_star,M,T_0C,IMAX
$ ,fiters
common/switches/iterate_rho_l
common/constants2/sigma,vv,omega_a
iterate_rho_l=.false.
T=T_ref/tau
call IUPAC_constants(N_i,r_i,s_i,k_i,j_i,b_i,c_i,f_i,g_i)
call calculate_x(x,tau,omega,omega_g,omega_l)
if ((x.gt.0).and.(x.lt.1)) then
call sub_dalpha_id_dtau(dalpha_id_dtau,tau,omega_l,f_i,g_i)
call sub_dalpha_dtau_omega(dalpha_dtau_omega,tau,omega_l,N_i
$ ,r_i,s_i,k_i,j_i,b_i,c_i,sigma,vv)
u_l=R*T*tau*(dalpha_id_dtau+dalpha_dtau_omega)
call sub_dalpha_id_dtau(dalpha_id_dtau,tau,omega_g,f_i,g_i)
call sub_dalpha_dtau_omega(dalpha_dtau_omega,tau,omega_g,N_i
$ ,r_i,s_i,k_i,j_i,b_i,c_i,sigma,vv)
u_g=R*T*tau*(dalpha_id_dtau+dalpha_dtau_omega)
u=u_g*x+(1−x)*u_l
else
call sub_dalpha_id_dtau(dalpha_id_dtau,tau,omega,f_i,g_i)
call sub_dalpha_dtau_omega(dalpha_dtau_omega,tau,omega,N_i,r_i
$ ,s_i,k_i,j_i,b_i,c_i,sigma,vv)
13/19
19−03−2007